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FOREWORD 


Volume  1  of  this  study  on  the  effect  of  long  wav  ts  in  the  lower 
Chesapeake  Bay  dealt  particularly  with  the  practical  and  engineering 
aspects  of  this  problem  and  was  based  in  part  on  new  theoretical  in¬ 
vestigations.  These  theoretical  investigations  are  presented  in  this 
Volume  II  as  various  appendices.  Most  of  them  are  original  contri¬ 
butions  to  the  science  of  hydrodynamics  of  long  waves.  They  may 
have  other  applications  than  those  presented  in  Volume  1  of  this  study. 
Volume  II  is  made  up  of  the  following  studies: 

Appendix  1:  Long  Waves  Generated  by  Nuclear  Explo¬ 

sions.  Types  of  Waves  and  Initial  Decay  in  Deep  Water.  (Dr. 
Basil  Wilson) 

Appendix  II:  Theoretical  Considerations  and  Computations 

for  Water  Waves  produced  by  Explosion.  (Dr.  Larry  Armijo 
and  Miss  Mary  Ann  Noser) 

Appendix  III:  Surface  Waves  Generated  by  Disturbance  on 

Sea  Bed  in  Constant  Depth  Open  Sea.  (Dr.  J.  A.  Hendrickson) 

Appendix  IV;  The  Principle  of  Superposition  and  Theory 
of  Cauchy- Poisson.(Or.  B.  Le  Mehaute) 

Appendix  V:  The  Shoaling,  Damping,  Breaking  and  Run¬ 

up  of  Long  Waves  over  the  Continental  Shelf.  On  Saturated 
and  Nonsaturated  Breakers.  (Or.  B.  Le  Mehaute) 


iii 


Appendix  VI:  The  Wave  Run-* up  by  the  Method  of  Charac¬ 

teristics.  (Dr.  B.  Le  Mehaute.  Dr.  J.  Freeman  was 
scientific  advisor  and  Mr.  R.  Grewal  did  the  numerical  com 
putations . ) 

Appendix  VII:  Two-Dimensional  Nonlinear  Wave  Motion 

in  an  Estuary.  (Drs.  J.  C.  Freeman  and  Larry  Armijo  and 
Miss  Mary  Ann  Noser) 
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ABSTRACT 


This  volume  aesembles  a  number  of  theoretical  studies  rele¬ 
vant  to  the  field  of  hydrodynamics  of  long  waves.  They  are: 

1.  Literature  survey  on  theoretical  and  experimental  informa¬ 
tion  on  the  problem  of  cylindrical  waves  generated  by  a  local 
disturbance,  with  application. 

2.  Application  of  the  theory  of  Kranzer  and  Keller  on  cylindrical 
waves  generated  by  explosion  to  three  magnitudes  of  power. 

3.  Theoretical  study  of  cylindrical  waves  generated  by  a  cylin¬ 
drical  upthrust  on  the  sea  bottono.. 

4.  The  principle  of  superposition  is  applied  to  the  Cauchy- 
Poisson  solution  for  determining  the  cylindrical  wave  motion 
due  to  a  finite  sea  surface  disturbance. 

5.  Wave  deformation  on  a  very  gentle  slope,  wave  damping  by 
bottom  friction.  Saturated  and  nonsaturated  breakers.  A 
survey  on  experimental  data  on  the  wave  run-up. 

6.  A  method  of  characteristics  is  presented  for  analyzing  the 
wave  deformation  over  a  gentle  slope,  wave  breaking,  spilling 
breakers,  bore,  and  run-up  on  a  dry  bed. 

7.  A  numerical  procedure  for  calculating  the  penetration  of  a  bore 
and  nonlinear  long  wave  into  estuaries  is  given  and  applied  in 
the  case  of  a  bore. 


XV 


INDEX  TO  APPENDICES 


APPENDIX  I . . 

LONG  WAVES  GENERATED  BY  NUCLEAR  EXPLOSIONS. 
TYPE  OF  WAVES  AND  INITIAL  DECAY  IN  DEEP  WATER. 
By  Dr.  Basil  Wilson 


APPENDIX  II . 

THEORETICAL  CONSIDERATIONS  AND  COMPUTATIONS 
FOR  WATER  WAVES  PRODUCED  BY  EXPLOSIONS.  By 
Dr.  Larry  Armijo  and  Miss  Mary  Ann  Noser 


APPENDIX  III.  . .  .  . . .  , 

SURFACE  WAVES  GENERATED  BY  DISTURBANCE  ON 
SEA  BED  IN  CONSTANT  DEPTH  OPEN  SEA.  By  Dr. 

J.  A.  Hendrickson 


APPENDIX  IV  . . . 

THE  PRINaPLE  OF  SUPERPOSITION  AND  THE  THEORY 
OF  CAUCHY -POISSON.  By  Dr.  B.  Le  Mehaute 


APPENDIX  V . 

THE  SHOAUNG,  DAMPING,  BREAKING  AND  RUN-UP  OF 
LONG  WAVES  OVER  THE  CONTINENTAL  SHELF.  ON 
SATURATED  AND  NONSATURATED  BREAKERS.  By 
Dr.  B.  Le  Mehaute 


APPENDIX  VI .  . 

THE  WAVE  RUN-UP  BY  THE  METHOD  OF  CHARACTERISTICS 
By  Dr.  B.  Le  Mehuate 


APPENDIX  VII . 

TWO-DIMENSIONAL  NONLINEAR  WAVE  MOTION  IN  AN 
ESTUARY.  By  Drs.  J.  C.  Freeman  .md  Larry  Armijo  and 
Miss  Mary  Ann  Noser 


APPENDIX  I 


LONG  WAVES  GENERATED  BY  NUCLEAR  EXPLOSIONS 


TYPE  OF  WAVES  AND  INITIAL  DECAY 
IN  DEEP  WATER 


By 


Basil  W.  Wilson 


1.  INTRODUCTION 


The  literature  on  the  effect*  of  surface  and  sub- surface  distur¬ 
bances  of  water  both  theoretical  and  experimental  is  quite  extensive 
and  this  tentative  survey  of  available  results  makes  no  claim  to  being 
exhaustive.  Satisfactory  treatment  of  the  initial  conditions  prevailing 
in  a  thermo- nuclear  underwater  explosion  has  not  yet  been  achieved, 
to  the  writer's  knowledge,  in  any  theoretical  analysis  or  laboratory 
experiment.  Both  mathematical  and  experimental  models  remain 
relatively  crude,  though  they  can  nonetheless  provide  useful  guide¬ 
lines  to  the  natural  behavior  of  the  archetype  event.  The  hydrody- 
namical  difficulties  of  achieving  an  effective  theoretical  simulation 
of  the  explosion  are  considerable,  though  progress  in  this  direction 
undoubtedly  will  come.  In  the  end,  however,  the  most  reliable  infor¬ 
mation  on  the  effects  of  nuclear  explosions  in  water  will  be  that 
secured  from  actual  prototype  experiments.  In  this  report  an  attempt 
will  be  made  to  examine  the  problem  of  wave  propagation  from  a 
nuclear  explosion  in  the  deep  ocean  (water  depth  circa  16,  000  ft.)  up 
to  the  point  that  the  waves  encounter  the  submerged  continental 
slope  of  the  nearest  land  mass.  In  accomplishing  this,  recourse  is 
made  to  such  theoretical,  laboratory  and  field  studies  of  impulsive 
water  waves  as  it  has  been  possible  to  analyze  within  the  time  avail- 


1. 


THEORIES  OF  RADIALLY  SYMMETRIC  DISPERSIVE  WAVES 


The  great  pioneering  theoretical  treatment  of  disturbances 
in  water  resulting  from  initial  elevations  or  depressions  of  lo¬ 
calized  extent  came  from  Cauchy  (1815)  and  Poisson  (1816).  Their 
combined  achievement,  now  made  classic  through  the  elegant  pre¬ 
sentation  and  extension  of  it  by  Lamb  (1904),  examined  the  effects 
in  water  of  infinite  depth  of  certain  shapes  of  paraboloidal  and 
ellipsoidal  depressions  or  elevations.  Lamb  generalized  the  solu¬ 
tion  in  the  two-dimensional  case  for  both  initial  elevation  and 
initial  impulse  by  making  use  of  the  Fourier  Integral  Theorem 
and  in  the  three-dimensional  case  by  use  of  Neumann's  (1862) 
Theorem.  As  our  present  interest  lies  really  in  the  three-dimen¬ 
sional  problem,  having  cylindrical  symmetry,  we  shall  give  only 
passing  consideration  to  two-dimensional  solutions. 

Within  the  limits  of  the  assumption  of  a  concentrated  point 
impulse  or  elevation  applied  to  the  surface,  Lamb  derives  an  exact 
hydrodynamical  solution.  The  solution  of  the  free  surface,  in  the 
form  of  an  infinite  series  (Lamb,  1932,  fz55),  has,  however, 

rather  limited  physical  significance  because  the  input  energy  is 
implanted  on  an  area  of  infinitely  small  extent.  Kelvin's  method 
of  st.ationary  phase  (1887),  applied  by  Lamb,  nevertheless  serves 
to  show  that  at  large  distances  from  the  source  the  wave  forms 
(initial  elevation)  and  (initial  impulse)  will  be  given  by 
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cotkt' 
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wherein  and  are  respectively  the  concentrated  elevation  and 
impulsive  pressure  at  the  point  source  per  unit  area, .  r  is  radial 
distance  from  the  source,  t  is  variable  time  and  p  and  g  have 
their  usual  meanings  of  fluid  density  and  gravitational  acceleration. 

The  period  T  of  these  waves,  if  t  »  T,  is 


T  » 


2tr 

<r 


4wr 


(1-2) 


and  their  wave  length  X  .  provided  r  »  X  ,  is 


X  = 


Birr' 


(1-3) 


making  the  wave  velocity  c, 


^  A 

c «  —  *  2 


=  2V 


(1-4) 


or  twice  the  group  velocity,  V. 
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Since  the  wave  length  and  period  df  the  waves  at  considerable 
distance  from  the  source  tend  to  change  very  slowly  with  result  that 
r/t  (the  group  velocity)  varies  only  gradually,  Eqs,  (I-l)  tend  to  show 
that  wave  amplitude  will  decline  approximately  as  r 

Terazawa  (1915)  applied  Lamb's  methods  to  the  case  of  an 

initial  disturbance  spread  over  a  finite  extent  of  the  free  surface  in 

infinitely  deep  water,  and  likewise  found  the  amplitude  decay  to  be 

proportional  to  r~^.  Terazawa  also  investigated  the  effect  of  an 

impulsive  explosion  at  a  finite  depth  h  below  the  surface  and  found 

the  initial  amplitude  of  the  wave  disturbance  at  surface  zero  to  be 
-3  /2 

proportional  to  h  .  The  same  result  is  deduced  by  Lamb  (1922) 

for  a  rather  different  model  of  explosion  in  which  an  abrupt  pressure 
rise  is  followed  immediately  by  a  gradual  fall.  If  the  initial  pressure 
rise  is  more  gradual  the  amplitude  of  surface  elevation  tends  to  vary 
as  h"^. 

The  high  frequency  of  occurrence  of  tsunami-generating  earth¬ 
quakes  off  the  Japanese  islands  has  led  Japanese  scientists  to  pay  con¬ 
siderable  attention  to  the  problems  of  waves  generated  by  impulsive 
movements  of  the  sea  bed  and  ocean  surface.  Prominent  among 
authors  who  have  contributed  to  analyses  of  these  problems  have  been 
Sano  and  Hasegawa,  Syono,  Homma,  Nakamura,  Sezawa  and  Kanai, 
Takahasi,  Ichiye,  Matzawa,  Miyoshi  and  Unoki  and  Nakano. 

We  shall  commence  by  considering  briefly  the  remarkable 
series  of  papers  of  Unoki  and  Nakano  (1953,  (i),  (ii),  (iii)  ),  which 
extend  the  work  of  Lamb  and  Terazawa  for  a  surface  disturbance  or 
impulse  of  finite  amount  spread  over  a  finite  area  in  deep  water  and 
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compare  the  results  with  observed  waves  from  a  volcanic  explosion. 
For  the  case  where  the  initial  elevation  is  a  uniform  piston-like  rise 
Q  of  the  surface  over  a  circle  of  radius  R  at  the  origin,  the  wave 
disturbance  is  described  by 

Jj  (KR)  C08(kr)  (1-5) 

provided  r  3i>  R,  this  being  the  asymptotic  solution  of  the  problem 
using  Kelvin's  powerful  method  of  stationary  phase.  In  Eq.  (1-5)  fi 
is  a  coefficient  of  "virtual  viscosity"  or  friction  coefficient  in  the  ex¬ 
ponential  time  decay  arising  from  the  assumption  that  friction  from 
eddy  viscosity  is  proportional  to  fluid  velocity;  is  a  Bessel 
function  of  the  first  order  which  modulates  the  last  cosine  term. 

The  equivalent  result  for  a  piston- like  impulse  of  uniform 
amount  I  imposed  on  the  surface  over  a  circle  of  radius  R  at  the 
origin  is 

“■  J.  (Kr)  (1-6) 

again  under  the  condition  r  ^  R. 

In  both  Eqs.  (1-5)  and  (1-6),  since  the  group  velocity  V  (=r/t) 
changes  only  slowly  for  large  values  of  r,  wave  amplitude  decay, 
discounting  the  exponential  decay  with  time,  is  proportional  to  r~^. 
The  system  of  waves  that  arise  in  these  two  instances  have  beats 
whose  modes  are  determined  by  zero  values  of  tho  Bessel  functions 
Jj  (kR). 
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Unoki  and  Nakano  successfully  applied  these  results  to  the  case 
of  the  Myojinsho  reef  submarine  volcanic  explosions  and  concluded 
that  the  wave  system  described  by  Eq.  (1-6)  best  fitted  the  observa¬ 
tional  data.  The  volcanic  explosions  (which  have  been  described  as 
not  dissimilar  to  small  nuclear  explosions)  thus  accorded  reasonably 
well  with  the  mathematical  model  of  a  uniform  cylindrical  surface  im¬ 
pulse.  The  authors  estimated  the  energy  of  the  explosions  as  being 
19 

from  3  to  8  X  10  ergs  (about  the  equal  of  a  1  kiloton  nuclear  blast), 
and  concluded  that  most  of  the  energy  went  into  wave  formation. 

From  a  comparison  of  observations  with  theory,  they  concluded  that 
R  had  the  value  2.2  km. 

In  some  cases  the  volcanic  explosions  gave  rise  to  wave  trains 
which  showed  no  beat  effects.  Unoki  and  Nakano  explain  this  on  the 
basis  that  the  initiating  impulse  was  sometimes,  probably  of  Gaussian 
form.  Thus  by  assuming  the  instantaneous  impulse  to  have  the  form 

I  (r)  *  I,e  7^  (l-O 

These  authors  shew  that  the  resulting  asymptotic  approximation  to 
the  wave  form,  i^j  ,  for  r  »  R  is 

In  this  case  there  is  an  absence  of  beating  and  wave  amplitude  merely 
decays  monotonicly. 
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We  may  note  that  Unoki  and  Nakano  also  elaborated  the  two- 

dimensional  case  of  wave  disturbances  originating  from  finite  source 

disturbances  and  showed  that  wave  amplitude  decay  in  deep  water  in 

these  circumstances,  remote  from  the  source,  is  proportional  to 
-1/2 

X  ,  where  x  in  this  case  defines  the  horizontal  distance.  Both 
Jeffreys  and  Jeffreys  (1956  Edn. ,  p.  517)  and  Eckart  (1948,  p,  409), 
however,  independently  show  that  for  this  two-dimensional  case  the 
wave  amplitude  near  the  front  of  the  train,  where  group  velocity 
approaches  the  value  d  being  the  water  depth,  declines  as 

r  so  that  the  front  of  the  wave  train  becomes  increasingly  more 

prominent  with  lapse  of  time  and  distance  in  the  dispersion  of  the 
system.  Eckart's  solution  shows  that  this  wave  front  is  an  amplitude 
modulation  of  a  sinusoidal  carrier  system  of  waves  whose  wave  length 
is  infinite  (at  least  for  the  usual  assumed  inconfipressible  water 
medium).  Effective  wave  length  near  the  front  of  the  train  is  thus 
dictated  by  the  modes  of  the  modulating  Airy  integral,  while  towards 
the  rear  the  sinusoidal  carrier  waves  increasingly  assume  dominwce 
in  defining  the  wave  length. 

Reverting  again  to  the  three-dimensional  problem,  we  find  that 
Kranzer  and  Keller  (1959),  again  confining  attention  to  initial  surface 
elevation  or  impulse,  but  introducing  the  influences  of  water  depth, 
d,  and  finite  areal  disturbance,  derive  for  the  case  of  initial  elevation 

£s (kd)  co»  (kr-<rt)  (1-9) 
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where  ?3{k)  is  the  Hankel  transform  of  the  function  Q(r)  describing 
the  initial  elevation  as  a  function  of  r  and  ^^(kd)  is  a  continuous 
function  of  kd  which  varies  in  value  from  about  1  when  kd  is  less 
than  tt/IO  (shallow  water)  to  an  asymptotic  value  (  ^^kd)  when 
kd  >  TT  (deep  water). 

For  the  comparable  case  of  initial  impulse  Kranzer  and 
Keller's  result  is 


17.  (kd)  sin  (kr-at)  (I-IO) 

in  which  r(k)  is  the  Hankel  transform  of  the  function  l(r)  describing 
the  initial  impulse  as  a  function  of  r  and  ^'^(kd)  is  a  continuous 
function  of  kd  which  approaches  the  value  unity  for  kd  <  tt/IO 
(shallow  water)  and  becomes  asymptotic  to  the  value  2 kd  when 

kd  >  TT  (deep  water). 

Kranzer  and  Keller's  derivations  are  the  asymptotic  solutions 
of  the  surface  disturbance  problems,  applicable  only  for  large  values 
of  kr  which  justify  the  use  of  Kelvin's  method  of  stationary  phase. 
They  follow  the  transform  techniques  applied  first,  apparently,  by 
Sneddon  (1951),  and  elaborated  by  Stoker  (1957).  The  Hankel  trans¬ 
forms  ?3(k)  and  r(k)  are  defined  as 


(i)  0  (k)  » /o  (r)r  Jo(kr)dr 

(ii)  I(k)»  J  I  (r)r  Jq  (kr)dr 


(l-li) 
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In  which  J  is  a  Bessel  function  of  sero  order, 
o 

In  I!qs.  (1-9)  and  (1-10)  a  distinction  has  now  to  be  recognized 
in  the  values  of  k  and  O’  over  those  given  by  Eqs.  (I-l  iii)  and 
(1-1  iv).  Here  k  is  defined  by  the  root  of  an  equation: 

^  (kd)  «  r/(t>/gd)  (1-12) 

while  O'  is  given  by 

O’*  «  gk  tonh  (kd)  (1-13) 

The  function  ^  (kd)  varies  between  asymptotic  values  of  I  for 
kd<  ir/10  (shallow  water)  and  1/(2  ^^d)  for  kd  >  ir  (deep  water). 
Thus  for  the  deep  water  case  with  (kd)  =  1/(2  Jkd),  k  and  O' 
assume  the  values  specified  in  £q.  (I-l). 

If  in  (1-9)  we  introduce  the  deep  water  value  at  ^kd,  and 
take  Q(r)  &  Q  for  0  <  r  <  R  with  Q(r)  s  0  for  r  R,  then  we 
obtain  exactly  the  same  amplitude  result  as  Unoki  and  Nakano  in 
Eq.  (1-5),  from  the  special  property  (in  this  case)  that 

j^Q(r)kr  Jo(kr)dr  a  0RJ,(kR),  (1-14) 

the  only  difference  residing  in  the  absence  of  the  friction  term  e~ 

In  like  manner  the  amplitude  terms  of  Eq.  (1-10)  reduce  identically  to 
those  of  Eq.  (1-6)  when  the  deep  water  value  =  ^/zkd  is  taken  along 
with  the  special  case  I(r)  s  I  for  0  <  r  <  R  and  I(r)  &  0  for  r  ^  R. 

Again  for  the  special  case  treated  by  Unoki  and  Nakano  in  Eq. 

(1-7)  we  find  that  Kranzer  and  Keller's  generalised  solution  applied  to 


1-9 


deep  water  (  =  v2  kd)  gives  identically  the  same  amplitude  re¬ 

sult  as  Eq.  (1-8)  (e~^^  excepted)  in  virtue  of  the  fact  that  the  Hankel 
transform  of  (1-7)  is 


I(k)«  2  loR* 


(1-15) 


If  now  we  apply  Eqs.  (1-9)  and  (I- 10)  to  shallow  water,  for  which 
kd  <  ir/10  and  si  1  si  the  wave  disturbances  at  a  long 
distance  from  the  source  for  the  cases  of  piston-like  surface  eleva¬ 
tion  Q  and  impulse  I  applied  at  the  origin  over  a  radius  R,  become, 
respectively: 


(i)  Vt  ^  Q 

(ii)  Vi  “  »in(kr-^<rt)  . 


(1-16) 


-  1  ■  -  1 /2  ' 

Dependence  of  upon  depth  is  thus  as  d  and  as  d  ; 

dependence  on  distance  in  both  cases  is  as  r~^. 

All  this  discussion  of  the  three-dimensional  form  of  the  Cauchy- 
Po'isson- Lamb  problem  suggests  that  wave  amplitude  decay  at  distances 
remote  from  the  source  is  proportional  to  r*^.  We  might  infer,  how¬ 
ever,  from  analogy  to  the  two-dimensional  problem  of  Jeffreys  and 
Eckart,  that  wave  amplitude  near  the  head  of  the  wave  train  will 
follow  another  law.  As  pointed  out  by  Munk  (1952)  this  in  fact  con¬ 
forms  to  r  The  reason  for  this  is  given  by  Takahasi  (1961) 
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who  points  out  that  the  method  of  stationary  phase  is  no  longer 

valid  in  the  neighborhood  of  the  wave  front.  For  there  the  third 

term  in  the  Taylor  expansion  of  (kr  -  c  t)  (cf.  Lamb,  jT  241)  is 

indeterminate  because  d^<r  /dk^  is  zero,  in  virtue  of  the  fact  that 

d  or/dk{=  r/t  =  is  a  constant.  Takahasi  shows  that  wave 
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height  at  the  front  is  proportional  to  r  for  the  case  of  waves 

resulting  from  a  piston-like  upthrust  of  the  sea-bed  over  a  radius 
R.  Since  the  behavior  of  dispersive  waves  at  a  large  distance  from 
the  source  is  practically  independent  of  the  nature  of  the  source  dis¬ 
turbance,  the  r  law  may  be  considered  to  prevail  in  general  at 

the  front  of  the  wave  train  and  the  r  ^  law  in  the  main  body  of  the 

waves.  Takahasi  (1961)  has  demonstrated  experimentally  that  the 
-5/6 

r  amplitude  decay  law  does  in  fact  prevail  at  the  leading  crest 

and  trough  of  wave  trains  generated  by  the  sudden  upthrust  of  a 
circular  portion  of  sea-bed  in  shallow  water. 

In  all  the  theoretical  results  presented  above,  the  nature  of  the 
wave  disturbance  near  the  source  is  undefined  because  of  the  mathe¬ 
matical  difficulties  of  describing  the  fluid  motions  near  the  moving 
boundaries.  Theoretically  complete  solutions  have,  however,  been 
obtained  by  Lamb  (1932,  jf  238,  239)  in  series  form  for  the  rather 
hypothetical  cases  of  concentrated  elevation  and  impulsive  pressure 
within  the  limits  of  linear  theory.  Recently,  the  near-source  dis¬ 
turbances  generated  by  various  forms  of  sea-bed  movement  which 
could  simulate  earthquake  displacements,  have  been  examined  in  a 
number  of  mathematical  treatments  by  Takahasi,  Ichiye,  Honda, 
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Nakamnura,  Keller  and  others.  Many  of  these  are  two-dimensional 
in  scope  and  are  therefore  not  of  direct  interest  or  application  to  this 
study.  In  other  respects,  however,  it  is  doubtful  whether  any  of 
their  models  can  be  considered  to  simulate  an  underwater  nuclear 
explosion  to  a  degree  which  could  justify  reliance  on  the  mathemati¬ 
cal  prediction  of  the  wave  forms  generated  near  the  source.  More 
satisfactory  mathematical  models  of  underwater  explosions  have  been 
developed  by  Penny  (1950),  Kirkwood  and  Seeger  (1950)  and  Fuchs 
(1952),  but  in  one  way  or  another  these  are  rather  poor  approxima¬ 
tions  to  the  nuclear  underwater  burst  near  its  source,  some  of  the 
features  of  which  have  been  illustrated  by  Snay  (1957)  and  discussed 
by  Cole  (1948)  and  by  Lane  and  Green  (1956). 

We  note  in  passing  that  the  theoretical  result  obtained  by 
Hendrickson  (see  Appendix  III)  for  a  piston- like  upthrust  of  the 
bottom  gives  a  wave- amplitude  decay  law  proportional  to  r"  ^  .  The 
result  is  quite  similar  to  that  of  Takahasi  (1961)  for  large  values  of 
kr  in  which  the  method  of  stationary  phase  is  applicable  for  deriving 
an  asymptotic  solution. 

In  regard  to  the  input  data  supplied  to  this  project  from  the 
analysis  of  Kaplan,  Wallace  and  Goodale  (1962)  (Fig.  1),  a  question 
really  needing  investigation  is  whether  the  application  of  the  asymp¬ 
totic  solution  of  Kranzer  and  Keller  (1959),  as  given  in  Eq.  (1-9), 
is  valid  at  the  relatively  small  values  of  kr  pertaining  when  r  is 
only  20  miles  from  surface  zero,  and  whether  full  reliance  can 
therefore  be  placed  on  the  results  of  such  an  application. 
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3. 


EXPERIMENTS  ON  DISPERSIVE  WAVE  BEHAVIOR 


In  Fig.  I-I  are  aeeembled  such  reeulte  of  obeervatlone  on  ex¬ 
plosion  waves  as  are  readily  available  in  "Effects  of  Nuclear  Weapons" 
(Glasstone,  1957«  1962)  from  the  underwater  nuclear  e3q>losion  tests 
in  the  Pacific. 

Shot  BAKER  of  a  20  Kiloton  nuclear  charge  at  Bikini  in  1946, 
in  a  water  depth  approximating  200  feet,  set  up  trains  of  waves  whose 
maximum  height,  H,  at  various  distances  shows  a  decay  propor¬ 
tional  to  r  ^  to  a  distance  of  about  1  mile  from  surface  zero  (Fig.I-1) 
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and  a  decay  thereafter  proportional  to  r  . 

In  Fig.  1-1  we  also  plot  the  standard  results  for  a  1  kiloton 
nuclear  imderwater  burst  as  presented  in  "Effects  of  Nuclear  Weapons." 
The  earlier  version  of  this  publication  (Glasstone,  1957)  suggests  a 
decay  law  of  wave  height  proportional  to  r~^  in  a  water  depth  of  85  ft. 
The  1962  version  on  the  other  hand  gives  the  decay  law  as 
over  a  range  of  distance  r  from  1  to  100  miles.  For  explosions  in 
deep  water  (  >  400  ft.)  the  decay  law  conforms  to  r’^. 

An  attempt  has  been  made  in  Fig.  1-2  to  represent  the  Bikini 
field  data  in  a  dimensionless  plot  for  comparision  with  certain  labora¬ 
tory  experiments  of  Johnson  and  Bermel  (1949).  Crest  elevation, 
above  still  water,  of  the  maximum  waves,  as  a  ratio  of  the  diameter 
of  the  crater  or  cylinder  of  the  burst  at  the  surface,  has  been 
plotted  against  dimensionless  distance  r/d,  in  which  d  is  the 
water  depth.  In  Fig.  1-2  the  first  part  of  the  BAKER  test  results 
(for  r/d  <  35),  no  longer  subscribes  to 
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FIGURE  1-2 

MAXIMUM  CREST  ELEVATION  OF  DISPERSIVE  IMPULSIVELY  GENERATED  WAVES 
AS  A  FUNCTION  OF  RADIAL  DISTANCE  FROM  THE  ORIGIN 


17^,/Dc  0C(r/d) 

because  in  the  initial  stages  of  the  wave  propagation  in  shal¬ 

low  water  was  in  general  greater  than  half  the  wave  height  H.  It 
was  necessary  to  infer  the  appropriate  value  of  from  Fig. 

1-3  which  was  originally  compiled  by  Reid  and  Bretschneider  (1953) 
and  subsequently  modified  slightly  by  Rechtin,  Steele  and  Scales 
(1957),  by  Bretschneider  (1958)  and  currently  by  the  present  writer. 
For  this  application  it  was  necessary  to  know  the  wave  period  T, 
which  was  calculated  from  Eq.  (1-13)  in  terms  of  the  known  dis¬ 
tance  r  and  elapsed  time  t.  We  may  point  out  here  that  Unoki  and 
Nakano  (1953,  i,  ii,  iii)  had  demonstrated  the  accuracy  of  this  for¬ 
mula  (Eq.  (1-2),  for  deep  water)  in  their  analysis  of  the  volcanic 
explosions  of  the  Myojinsho  reef. 

The  slope  of  the  latter  part  of  the  curve  representing  the 
BAKER  test  in  Fig.  1-2  accords  with  the  decay  law  because 

the  waves  for  r/d  >35  are  largely  oscillatory  with  iji^/H£S£  0.5. 
The  field  experiments  of  Van  Dorn  (1961),  related  to  nuclear  under¬ 
water  explosions  in  the  Pacific  in  1956,  show  that  relative  wave 
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amplitude  declines  as  (r/d)  over  a  range  of  values  of  the  latter 

parameter  from  100  to  1000,  Although  the  standard  of  reference  of 
relative  amplitude  is  not  given  by  Van  Dorn,  it  seems  reasonable 
to  infer  from  his  results  that  the  BAKER  test  results  of  Fig.  1-2 
which  extend  to  about  r/d  =  100  could  be  projected  at  the  same  slope 
(  OC  as  far  as  r/d  =  1000.  The  validity  of  this  step,  however. 
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may  be  open  to  queation,  on  grounds  that  different  values  of  d/D^ 
(alternatively  d/  X  )  are  involved. 

The  experiments  of  Johnson  and  Bermel  (1949)  measured  the 
characteristics  of  waves  generated  impulsively  by  horizontal  cir¬ 
cular  discs  falling  vertically  from  different  heights  into  still  water. 
The  data  of  their  Fig.  1-3  has  been  recalculated  and  adapted  to  the 
requirements  of  our  Fig.  1-2.  Although  individual  test  results  show 
considerable  variability  there  is  a  remarkable  overall  parallelism 
with  the  Bikini  data  e'.  en  to  the  extent  of  an  indicated  slope  change 
at  r/d  ^  35.  For  values  of  r/d  <  35  particular  test  results  show 
higher  than  average  and  some  lower  than  average  rates.  The  fact 
that  initial  decay  rates  greater  than  r~  ^  were  found  in  particular 
Instances  lends  emphasis  to  the  need  for  knowing  more  about  the 
special  conditions  that  may  lead  to  energy  dissipation  which  could 
detract  from  the  energy  available  for  wave  formation. 

Generally  speaking,  the  experimental  results  on  wave  ampli¬ 
tude  decay  satisfactorily  confirm  the  theoretical  predictions  of 
Section  2.  The  experiments  tend  to  show  that  in  a  dispersive  system 
of  linear  waves,  with  beat  characteristics,  the  decay  law  r  ^ 
prevails  until  the  distance  becomes  sufficiently  great  to  give  greater 
prominence  to  the  front  of  the  wave  train  and  render  invalid  the  theo¬ 
retical  approximation  of  stationary  phase,  thus  causing  the  decay  law 
to  change  to  one  proportional  to  The  transition  is  probably 

gradual  even  though  the  experiments  suggest  quite  a  sharp  change 
at  r/d  Si  35. 
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RATIO  OF  ELEVATION  OF  WAVE  CREST  ABOVE  STILL 
WATER  LEVEL  TO  TOTAL  WAVE  HEIGHT,  «/H 


FIGURE  1-3 

RATIO  OF  ELEVATION  OF  WAVE  CREST  ABOVE  STILL  WATER 
LEVEL  TO  TOTAL  WAVE  HEIGHT,  '7o/H,  AS  A  FUNCTION  OF 
RELATIVE  DEPTH  d/T^,  AND  RELATIVE  HEIGHT  H/T^ 
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4. 


TYPES  OF  WAVES  IN  INPULSIVELY  GENERATED  SYSTEMS 


Th«  data  of  Figs.  1>1  and  1-2  all  pertain  to  long  waves  in 
shallow  water.  In  Eqs.  (1-9)  and  (I- 10)  we  have  already  noted  that 
the  Fourier  transforms  and  T  and  the  functions  ^  and  Uf. 
are  dependent  on  kd.  or  alternatively,  d/K  ,  the  relative  depth 
referred  to  wave  length.  In  consequence,  the  type  of  wave  that  will 
result  from  any  given  initiating  disturbance  will  be  a  function  of  the 
relative  depth  d/X  .  Eqs.  (1-9)  and  (I- 10)  show  however  that  the 
carrier  waves  at  a  great  distance  from  the  source  (large  kr)  are 
siniply  sinusoidal.  This  obviously  is  not  necessarily  true  in  general 
of  the  Waves  near. the  source.  To  investigate  this  aspect  we  have 
recourse  to  Fig.  1-4  (adapted  from  Wilson  (1962)  ).  This  portrays 
the  area  of  existence  of  waves  of  different  types  in  a  plot  of  H/  X^ 
versus  d/ X  ^,  in  which  X^  is  the  equivalent  deep-water  wave 
length  of  any  given  wave  type.  Waves  are  theoretically  unable  to 
exist  outside  of  the  shaded  areas,  which  are  bounded  by  the  break¬ 
ing  wave  criteria  of  McGowan  (long  waves)  and  Michell  (short  waves), 
and  an  intermediate  limit  largely  defined  by  experiment  and  inter¬ 
polation;  Short  waves  are  normally  defined  as  those  for  which  the 
relative  depth  d/X  (=d/  X^)  >  0.5.  Long  waves  are  usually 
interpreted  as  those  for  which  d/ X  <  0.05  (d/X  ^  <  0.015). 

In  Fig.  1-4  a  set  of  isolines  of  the  dimensionless  parameter 

2  3 

(  X  /d  ),  in  which  <n  is  wave  crest  elevation  above  still  water 
o  'o 

level,  is  shown  crossing  the  diagram  in  diagonal  fashion.  The  signi¬ 
ficance  of  this  parameter  has  been  stressed  by  Ursell  (1953)  who 
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FIGURE  1-4 

CHARACTERISTICS  OF  IMPULSIVELY  GENERATED  WAVES  AS 
DETERMINED  EXPERIMENTALLY,  IN  RELATION  TO  THEIR 
THEORETICAL  PREDICTION 
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EQUIVALENT  VALUES  OF  H/T  -  (ff/«,e 


points  out  that  linear  wave  theory,  such  as  that  of  Airy,  is  normally 

2  3 

applicable  only  if  (  X  /d  )  «  1  and  H/  X^  is  also  small. 

Ursell's  analysis  showed  that  long  waves  are  inevitably  non- 
2  3 

permanent  when  (  X  /d  )  >>  1.  The  crests  of  such  waves  ad¬ 
vance  more  rapidly  than  the  troughs  and  finally  reach  a  condition  of 
instability  which  manifests  itSelf  in  breaking  or  bore-formation. 

Ursell  was  disposed  to  consider  that  the  solitary  wave  of  permanent 
type  could  exist  at  or  near  (  X  /d  )  =  1.  The  present  writer, 
however,  has  compiled  evidence  (Wilson,  1962)  to  show  that  the 
most  favored  range  for  the  existence  of  solitary  waves  is  in  the  belt 

shown  hatched,  or  more  broadly,  between  isolines  of  (  7)^  X  /d  ) 

2  3 

from  10  to  35,  In  the  zone  of  (  X  /d  )  from  1  to  10  cnoidal 
waves  of  permanent  type  are  possible,  which  overlap  permanent  waves 
derived  from  Stokes  theory  of  second  or  higher  orders.  Beyond 
(  7^^  X^/d^)  >  35,  waves  are  likely  to  be  non- permanent. 

Since  dispersive  waves  are  composite  of  the  interference  effect 
of  a  broad  spectrum  of  frequency  components,  the  relevance  of 
Fig.  1-4  may  at  first  seem  questionable,  as  Fig.  1-4  really  has  ap¬ 
plication  to  uniform  wave  trains.  However,  the  theory  of  Cauchy- 
Poisson-Lamb,  as  elaborated  by  Unoki  and  Nakano  and  by  Kranzer 
and  Keller,  clearly  shows  that  the  resultant  type  of  wave  in  a  disper¬ 
sive  wave  train  has  all  the  characteristics  of  an  equivalent  wave  in 
a  non- dispersive  system  at  any  particular  moment  of  its  existence, 
before  influences  of  distances  and  time  bring  about  change.  Conse¬ 
quently  it  is  of  importance  to  recognize  the  type  of  wave  likely  to 
exist  in  a  dispersive  system,  particularly  in  regions  close  to  the  source 
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where  available  theories  are  difficult  to  evaluate. 

Here  we  resort  to  a  consideration  of  the  experiments  of  Prins 
(1956,  1958)  who  generated  waves  in  a  laboratory  flume  from  sud¬ 
denly  released  initial  elevations  or  depressions  of  finite  extent  in 
water  of  constant  depth.  The  initial  waves,  near  the  source,  gen¬ 
erated  under  different  conditions  of  water  depth  d,  height  of 
elevation  and  length  of  initial  elevation,  have  been  plotted  in  Fig. 
1-4,  each  plotted  point  referring  to  a  different  test.  Prins  gave 
valuable  information  on  the  type  of  wave  initiated,  shown  in  the 
legend  of  Fig.  1-4.  It  is  therefore  of  considerable  interest  to  find 
that  his  types  lie  within  specific  belts  of  the  parameter  (  7) ^  X  /d  ) 
and  that  the  occurrence  of  single  solitary  waves  (followed  usually 
by  a  trailing  train  of  small  oscillatory  waves)  agrees  remarkably 
well,  with  the  zone  of  existence  for  solitary  waves  that  we  have 
already  specified.  For  (  7^^  X^/d^)  >  40  Prins  found  waves  of 
"complex- solitary"  type,  and  under  conditions  which  brought  the 
leading  wave  characteristics  close  to  the  McGowan  limiting  criter  ia 
the  wave  was  inevitably  a  type  of  bore. 

Also  shown  plotted  in  Fig.  1-4  are  some  experimental  results 
of  Wiegel  (1955),  who  generated  waves  under  laboratory  conditions 
by  causing  submerged  block  masses  to  slide  down  slopes  of  varying 
steepness,  thereby  simulating  submarine  landsliding  in  earthquakes 
Most  of  these  waves  were  of  oscillatory  character  and  plot  in  a 
region  of  Fig.  1-4  which  would  identify  them  as  Stokes  waves. 

The  experimental  results  of  Johnson  and  Bermel  (1949)  have 
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also  been  plotted  in  Fig.  1-4  on  the  assumption  that  the  effective 
initial  wavelength  X  of  the  primary  wave  near  the  source  was 
four  times  the  radius  of  the  circular  plate  dropped  on  the  water  sur¬ 
face.  The  data  are  found  to  fall  in  a  zone  which  would  classify  the 
waves  as  breaking  complex  solitary  or  single  solitary  waves. 

Finally,  the  initial  waves  of  the  nuclear  underwater  explo¬ 
sion  test  BAKER  has  been  located  in  Fig.  1-4  by  determining  the 
applicable  values  of  H/T^  and  d/T^  and  this  plotted  point  is 

found  to  lie  in  the  unstable  (breaking)  complex- solitary  or  single- 

2  3 

solitary  zone  for  which,  (17^^  /d  )  is  from  about  30  to  40.  The 
BAKER  test  result  accords  remarkably  well  with  Johnson  and  Bermel's 
laboratory  data  which  were  scaled  to  produce  results  approximately 
equivalent  to  the  field  explosion  of  1946.  That  the  initial  waves  from 
the  BAKER  test  were  actually  of  solitary  type  is  confirmed  by  the 
following  quotation  from  Glasstone  (1962): 

"Observations  of  the  properties  of  the  waves  indicated  that 
.  the  first  wave  behaved  differently  from  succeeding  ones  in 
that  it  was  apparently  a  long  solitary  wave  generated  directly 
by  the  explosion,  receiving  its  initial  energy  from  the  high 
velocity  outward  motion  of  the  water  accompanying  the  ex¬ 
pansion  of  the  gas  bubble.  " 

Assuming  that  the  primary  wave  form  near  the  explosion  source 

in  a  nuclear  underwater  burst  is  solitary  or  complex- solitary,  as 

suggested  by  Fig.  1-4,  the  extent  to  which  the  wave  form  will  change 

2  3 

will  obviously  be  dictated  by  the  prevailing  value  of  (  17^  X  /d  )  in 
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which  each  of  the  variables  of  this  parameter  is  a  function  of  r,  the 
radial  distance  of  progression.  If  the  depth  is  constant  the  parameter 
will  depend  only  on  the  product  ^X(r)  j  Since  “f]  is 

a  decaying  function,  and  X(r)  an  increasing  function,  of  the  distance 
r,  the  value  of  the  product  will  be  dictated  by  whichever  of  these 
variables  prevails.  If  the  product  should  remain  constant.  Fig.  1-4 
shows  that  wave  form  will  remain  unchanged.  Thus  the  wave  could 
propagate  as  a  distended  and  flattened  solitary  wave.  If  the  product 
declined  in  value,  the  wave  form  would  become  cnoidal  and  even¬ 
tually  oscillatory.  On  the  other  hand,  if  it  increased  in  value,  the 
wave  form  would  become  complex  and  unstable,  finally  leading  to 
bore  formation.  Since  depth  is  seldom  constant  in  wave  propagation, 
this  variable  (d)  obviously  will  have  a  powerful  influence  on  the 
value  of  the  parameter  (  X^/d^)  and  hence  on  the  ultimate  wave 

form.  As  the  depth  decreases  towards  the  coast  the  parameter  must 
inevitably  increase  and  cause  the  wave  form  to  assume  the  unstable 
shape  that  finally  leads  to  bore  or  breaker. 

The  nature  of  the  decay  of  with  distance  r  has  already 

received  considerable  attention  in  Sections  2  and  3  of  this  appendix. 
There  remains  to  consider  how  the  wave  length  X  and  period  T 
of  the  waves  are  modified  by  the  dispersion. 

Before  leaving  this  section  it  may  be  noted  that  the  input  data 

originally  supplied  to  this  project  gave  the  maximum  wave  height  at 

r  =  20  miles  as  n  =  H/2  =  388  ft.  (calculated  from  Kranzer-Keller 
'  o 

theory  for  d  =  18,  000  ft,  ),  with  a  corresponding  period  T  =  1. 2 
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minute*.  This  give*  ua  H/T^  =  0.  15  and  d/T^  =  3.47  (ft.  /*ec.^ 
unit*).  When  this  point  i*  located  in  Fig.  1-4  (not  actually  shown) 
it  fall*  in  a  region  which  would  not  seem  to  justify  the  use  of 
linear  theory  (the  premise  upon  which  Kranzer- Keller,  Unoki- 
Nakano,  and  Cauchy- Poisson- Lamb  theories  are  all  founded). 
About  the  highest  value  that  H/T^  could  have,  in  order  not  to 
violate  the  necessary  conditions  of  linear  theory,  would  be  0.042 
(cf.  Wilson,  1962,  Fig.  28).  This  is  an  additional  reason,  per¬ 
haps,  beyond  that  expressed  earlier,  for  questioning  the  use  of 
Kranzer- Keller  theory  in  derivation  of  the  input  data. 


5.  INCREASE  OF  WAVE  PERIOD  WITH  DISPERSION 

From  Eq.  (1-9)  it  is  clear  that  a  wave  crest  will  be  encoun¬ 
tered  whenever  cos(kr  -  cr  t)  s  1.  This  will  happen  whenever 


(i) 

kr-<rt  =  ±  Zirm 

(ii) 

m-  0,l,2,3,.. . 

Eq.  (1-17  i)  may  be  written  in  the  form 

kr(l-  2  TT  m 


(1-18 


and  since  by  definition 


(i) 

c *  X/T  *  c/k 

(ii) 

(1-19) 
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V  being  the  group  velocity,  or  rate  of  progression  of  the  wave 
groups,  Eq.  (1-18)  is  adaptable  to  the  form 


kd(l-^)=  ± 


(r/d) 


{I--J0) 


We  have  already  noted  that  waves  generated  from  an  initial  distur¬ 
bance,  although  composite  of  a  spectrum  of  frequencies,  tend  through 
interference  effects  to  assume  momentarily  the  form  of  waves  of 
non- dispersive  type,  subject  to  the  same  conditions  under  which  the 
latter  can  exist.  Thus,  as  the  linear  theory  of  Kranzer- Keller 
predicts 


(ii) 


I  + 


2kd 


sin  h  2kd 


] 


(l-cil) 


Since  V  in  general  is  less  than  c,  we  must  take  the  negative  sign 
with  the  right  hand  side  of  (l-<iO).  Eq.  (1-20)  thus  transforms  to 


sinh  2kd-2kd1  27rm 
sinh2kd+2kdJ  (r/d) 


Van  Dorn  (1961)  has  evaluated  the  function  of  (kd)  forming  the  left 
hand  side  of  Eq.  1-22),  It  is  sufficient  to  note  here  the  two  extremes 
of  its  values. 

First,  for  deep  water  conditions  (kd  >  V  ),  (1-22)  becomes 
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a-^3) 


and  since  for  the  same  conditions  c^  =  g/k  or  tr  =  ^^k  from 
Eq.  (1-13),  elimination  of  k  in  (I-Z^)  in  favor  of  cr  (  =  2  tt/T) 
yields 


Eq.  (1-24)  shows  that  the  period  of  those  waves  for  which  d/  X  >  -j, 
which  would  inevitably  include  most  of  the  trailing  waves  in  the  dis¬ 
persive  system  and  possibly  the  bulk  of  the  waves  if  the  depth  is 
great  and  the  explosion  small,  will  increase  as  the  square  root  of 
the  distance  (r^^^). 

For  the  opposite  extreme  the  function  (kd)  of  Eq.  (1-22) 
assumes  the  asymptotic  value  of  f  (kd)^/3]  /[  1  +  (kd)^/3]  when 
kd  <  tt/IO,  applicable  to  shallow  water  conditions.  In  this  case 
then  Eq.  (1-22)  reduces  to 


(kd)^~  2rrm 
3  (r/d) 


(1-25) 


For  this  condition  Eq.  (1-13)  yields  O’  ^  =  gdk*^  so  that  (1-25)  re 
duces  to  the  form 


(1-26) 
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and  the  inference  is  that  wave  period  increases  as  the  cube  root  of 
the  distance  (r^^^)  when  depth  conditions  are  such  that  d/X  <  1/20. 
This  result  agrees  with  Van  Dorn  (1961)  and  accords  with  the  two- 
dimensional  theory  of  Eckart  (1948,  p.  409). 

In  Eqs.  (1-24)  and  (1-26)  an  anomaly  arises  for  m  =  0  at  the 
front  of  the  wave  train  for  which  the  wave  period  and  wave  length 
are  theoretically  infinite,  for  here  from  Eq.  (1-17  i),  (kr  -  a  t)  =  0 
or 


(i) 


C»  o-/k  =  r/t=  V 


(ii) 


c=^/5d 


(1-27) 


The  anomaly  is  implicit  in  the  Kranzer- Keller  theory  since  Eq.  (1-27) 
is  in  agreement  with  (1-12),  provided  (kd)  =  1  with  kd  =  0.  It 
brings  to  focus  the  statement  made  earlier  in  reference  to  Eckart's 
(1948)  solution  of  a  two-dimensional  dispersive  wave  system,  that 
the  modulating  factor  apparently  becomes  the  criterion  determining 
the  wave  length.  It  is  of  interest  to  quote  Eckart  in  this  regard; 

"The  situation  near  this  front  it  very  peculiar.  It  may  be 
described  as  a  carrier  wave  which  is  amplitude  modulated. 
However,  contrary  to  the  customary  case,  the  carrier  has  a 
longer  (infinitely  longer!)  wave  length  than  the  modulation. 
Consequently,  the  empirically  determined  "wave  length" 
bears  no  relation  whatever  to  the  spectrum  of  the  disturbance. 
Such  wave  fronts  are  of  considerable  importance  as  they  con¬ 
stitute  the  seismic  sea  waves,  commonly  called  'tidal  waves'.  " 
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To  resolve  this  problem  we  shall  attempt  an  approximate 
analysis  based  on  Eq.  (I-16i).  which  gives  the  wave  form  for 
shallow  water  conditions  (applicable  to  the  wave  front)  for  a 
piston-like  surface  elevation  Q  applied  over  a  radius  R  at  the 
origin.  If,  in  conformity  with  Eckart,  the  modulator,  Jj(kR)/kd 
in  this  case,  becomes  conditional  in  prescribing  the  effective 
wave  length  and  period  at  the  front  of  the  dispersive  system, 
then  it  can  be  seen  that  for  small  values  of  kR  the  above  term 
approximates  to 


Jl(kR) 

kd“ 


(kR)^  .  (kR)'* 

s—  +  "W 


>1' 


(T-^8) 


so  that,  approximately 


kR  =  4  [  'I'  ]  (1-29) 

if  represents  the  function  J^(kR)/kd. 

Crests  of  the  carrier  waves  near  the  front,  however,  occur 
whenever  cos(kr  -  <r  t)  =  1,  so  that  the  same  condition  (1-25) 
prevails  as  before,  giving 


kd  = 


Z  tr  m 

WTSf 


1/3 


(1-30) 


On  multiplying  Eqs.  (1-29)  and  (1-30)  and  noting  that 


Igd 
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for  small  values  of  kd,  as  pertain  to  the  wave  front,  we  find 


(1-31) 

The  series  approximation,  Eq.  (I-<!8),  is  valid  so  long  as 
kR  is  less  than  about  2.0  and  since  the  Bessel  function  Jj(kR) 
reaches  its  first  maximum  value  at  kR  =  1 . 84,  it  may  be  con¬ 
cluded  that  Eq.  (1-31)  is  pertinent  for  values  of  m  (^  1,  2,  3.  , .  .  ) 
up  to  the  maximum  waves  in  the  first  group  so  long  as  kd  <  1/20. 
For  kR  =  2.0,  Eq.  (1-29)  gives 

(1-32) 

resulting  in  a  simplification  of  Eq.  (1-31)  to 

If  the  p-th  wave  is  the  highest  in  the  first  beat  then  its  period 
will  be  given  approximately  by  Eq.  (1-32)  for  m  =  p. 

The  general  picture  we  gain  of  increr.se  of  wave  period 
with  distance  on  dispersion  from  the  source  is  now  as  follows: 
first,  the  size  of  the  initial  waves  set  up  by  an  underwater  explo¬ 
sion  or  other  disturbance  is  a  function  of  the  magnitude  of  the 
disturbance,  R.  Wave  length  and  period  would  appear  to  increase 
very  slowly  with  distance  according  to  r^^^  with  long  waves  some 
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distance  behind  the  front  the  period  increase  would  conform  to 

1/3 

r  and  with  short  waves  near  the  rear  of  the  train,  or  in  other¬ 
wise  exceptionally  deep  water,  the  period  increase  would  follow 
the  law  r  ^ 

6.  EXPERIMENTAL  EVIDENCE  FOR  WAVE  PERIOD  INCREASE 
In  a  classic  paper  on  the  subject  of  wave  period  increase, 

Munk  (1947)  had  indicated  that  the  wave  period  of  tsunamis  must 
bear  some  relationship  to  the  original  disturbance,  because  of  an 
indeterminateness  of  the  function  in  the  general  solution  of  the 
governing  differential  equation: 

^  +  V  =  0  (1-34) 

O  t  OX 

This  is  patent  from  the  data  drawn  from  his  paper  and  plotted  in 
Fig.  1-5,  which  relates  the  dimensionless  quantities  T  and 

r/d.  Clearly  the  tsunamis  of  April  1,  1946  (Aleutian  Trench), 

April  1.1,  19<i3  (Kamchatka),  and  November  10,  1922  (Chile),  obey 
different  laws  of  period  increase  with  distance,  as  also  the  tsunami 
originating  from  the  submarine  volcanic  eruption  at  Myojinsho  Reef, 
on  March  11,  1953  (Unoki  and  Nakano,  1953). 

For  the  tsunami  of  April  1,  1946,  we  have  extracted  from  the 
data  given  by  Munk,  the  periods  of  the  waves  20  mins,  behind  the 
wave  front,  as  recorded  at  Valparaiso,  Matarani,  La  Jolla  and 
Honolulu.  Also  assembled  in  Fig.  1-5  are  period  data  for  the  leading 
waves  obtained  from  Takahasi's  (1961)  model  experiments  on  a 
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DIMENSIONLESS  WAVE  PERIOD  -  T-/g7d 


FIGURE  1-5 

PERIOD  OF  THE  LEADING  WAVE  IN  A  DISPERSIVE,  IMPULSIVELY 
GENERATED  SYSTEM  AS  A  FUNCTION  OF  RADIAL  DISTANCE  FROM 
THE  ORIGIN 
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piston-like  upheaval  of  sea  bed  and  from  Van  Dorn's  (1961)  mea¬ 
surements  of  underwater  nuclear  explosion  waves  from  Bikini  Atoll, 
as  measured  at  Ailinginae,  Eniwetok,  Wake  and  Johnston  Islands 
in  the  Pacific.  It  is  quite  fortuitous  that  all  these  last  mentioned 
data  are  in  virtual  alignment  in  Fig.  1-5,  suggesting  that  from  a 
dimensional  standpoint  the  source  disturbances  in  each  case  must 
have  been  quite  similar.  What  is  of  special  interest  is  that  for 
small  values  of  r/d  Takahasi's  data  show  T  Of  r^^^,  in  agree¬ 
ment  with  Eq.  (1-33)  while  Van  Dorn's  and  Munk's  data  conform 

1/4 

almost  exactly  to  a  law  T  Of  r  What  is  further  of  extraordinary 

interest  is  that  the  change  from  to  r^^^  occurs  at  r/d  Ci  35, 

suggesting  some  parallel  with  the  height- change  law  of  Fig.  1-2. 

We  may  recall  too  that  the  experimental  height-change  law  (  cC  r”^) 

for  small  r/d  was  in  accordance  with  the  prediction  of  the 

theoretical  statiqnary- phase  asymptotic  approximation  and  note  the 

similarity  here  that  the  experimental  period- change  law  (  Otr^^^) 

for  small  r/d  is  also  in  accord  with  the  theoretical  stationary- 

phase  asymptotic  approximation  (Eq.  (1-33)  ).  Beyond  r/d  —  35, 

the  height  change  law  (  OCr^^^)  exhibited  the  invalidity  of  the  approx 

imation.  It  seems  evident  therefore  that  the  same  invalidity  is  res- 

1/4 

ponsible  for  the  period-change  law  following  a  new  pattern  (r  ). 
Data  from  the  two-dimensional  experiments  of  Prins  (1956,  1958) 
have  not  been  invoked  in  Fig.  1-5,  as  they  would  be  inappropriate. 

Eq.  (1-33)  can  be  expressed  alternatively  as 


f 


(i) 


(ii) 


K^{ 


(1-35) 


It  is  of  interest  to  compare  the  theoretical  value  of  K  given  by 
(1-35  ii)  with  the  value  K  (=  8.2),  found  from  Fig.  1-5  as  appli¬ 
cable  to  Takahasi's  experiments.  For  this  purpose  we  take  p  =  1 
for  the  first  wave  and  insert  R  =  6  ins,  d  =  1.875  ins,  the  dimen¬ 
sions  of  piston  radius  and  water  depth  respectively  in  the  model 
experiments.  Eq.  (1-35  ii)  then  yields  K  =  7.03  which  is  of  the 
correct  order  of  magnitude.  Since  R  should  really  be  the  radius 
of  the  initial  surface  disturbance  which  probably  exceeded  6  ins.  , 
the  congruency  could  probably  be-improved.  Eqs.  (1-33)  dr  (1-35) 
thus  appear  to  be  a  reliable  prediction  of  initial  wave  period 
increase. 

From  the  alignment  of  the  April  1,  1946.  tsunami  and  the  1956 
REDWING  nuclear  explosion  data  with  that  of  Takahasi's  (Fig.  1-5), 
we  conclude  that  the  same  R/d  value  of  about  o.2  prevailed.  For 
the  tsunami  this  would  imply  an  initial  radius  of  surface  disturbance 
R  Si  52,  200  ft,  (d  S!  16,  320  ft. )  or  8.  60  n.  mi.  Interpreting  this 
as  related  to  an  earthquake  fault- length  of  19.  2  n.  mi.  or  35.6  km, 
the  equivalent  earthquake  magnitude  would  be  about  7.40  (cf.  Wilson, 
1962  (Fig.  8)  ),  whereas  the  actual  earthquake  magnitude  was  in 
fact  7.4.  For  the  REDWING  explosion  series,  the  implication  is 
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that  if  the  waves  had  emanated  from  an  explosion,  unobstructed 

by  atolls  or  reefs,  the  bubble  crater  diameter  2R)  would 

have  been  6.4d,  or  about  6  times  the  water  depth.  For  d  ~  200  ft, 

D  would  thus  have  been  of  the  order  of  1280  ft. 
c 

7.  CALCULATION  OF  DESIGN  WAVE  HEIGHT  AND  PERIOD 
AT  THE  CONTINENTAL  SLOPE 

Despite  any  reservations  made  regarding  the  input  design 
data,  we  proceed  to  use  them  in  conjunction  with  the  laws  of  height 
and  period  change  with  radial  distance  of  dispersion,  justified  in 
this  appendix,  in  determining  the  wave  height  and  period  that  would 
prevail  at  the  continental  slope  off  Chesapeake  Bay. 

The  law  of  maximum  wave  height  change,  following  the  dis¬ 
cussions  of  Sections  2  and  3,  may  be  generalized  as 

^  '  On  (1-36) 

in  which  n  is  a  numerical  exponent  and  the  corresponding  pro¬ 

portionality  constant.  It  is  convenient  to  write  this  in  the  form 

Y=a„X-"  (1-37) 

where 


(i) 


(1-38) 


(ii) 


X  »  r/d 
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In  order  to  apply  Eq.  (1-36)  derived  for  constant  depth,  to  conditions 
under  which  the  depth  is  a  function  of  r,  namely  d  =  d(r),  it  is 
necessary  to  integrate  di),/dr  as  a  function  from  which  X  has 
been  eliminated,  (cf.  Wilson,  1961).  Thus  by  differentiating  Eq. 
(1-37),  regarding  d  as  constant,  and  eliminating  X  in  favor  of  Y  by 
use  of  (1-37),  we  arrive  at  the  expressions 


(i) 

dr  d(r)  ' 

(ii) 

L  At' 

(1-39) 

The  law  of  period  increase,  following  the  discussions  of  Sections  5 

and  6,  may  be  generalized  in  a  similar  way 

as 

(1-40) 

with  m  a  numerical  exponent  and  the 

corresponding  constant 

of  proportionality. 

Writing  this  as 

(1-41) 

in  which 

(i) 

Z=TVf 

(1-42) 

(ii) 

X*  r/d 
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we  may  follow  the  same  general  procedure  as  for  wave  height  in 
deriving 


(i) 

if.  g-4  [d(r)]- 

(ii) 

Z'  *  m  5*  Z 

(1-43) 


In  case  III,  the  given  data  prescribes  a  half-wave  amplitude 

=  388  ft.  at  r  =  20  miles  resulting  from  a  bubble  crater  dia¬ 
meter  s  20,  800  ft.  For  an  assumed  point  of  detonation  at 
latitude  35°  30%  longitude  65°  27',  off  the  east  coast  the  water 
depth  profiles  along  great  circle  paths  directed  towards  Chesapeake 
and  New  York  Bays  are  as  shown  in  Fig.  1-6,  with  depth  at  the  de¬ 
tonation  point  d  15,  500  ft.  The  leading  wave  period  is  given  as 
Tj  =  2.75  mins,  and  that  of  the  highest  wave  (crest  No.  4  or  trough 
No.  5)  as  =  1.2  mins. 

With  this  information  we  find  =  1.87  x  10  ^  and 

Tj  =  3.  87  at  r/d  =  6.77.  These  points  are  located  on  Figs. 

1-2  and  1-5  and  design  curves  have  been  constructed  through  them 
parallel  to  the  theoretical  empirical  relationships  found  to  correspond 
with  observations.  The  same  indicated  change  of  slope  at  r/d  =  35 
has  been  adopted  in  both  cases. 


*  This  location  was  the  first  originally  prescribed  for  the  project  but 
was  later  amended  to  one  closer  to  the  shore  (see  Chapter  V  -  Vol.  1 
Part  B).  The  calculation  here  given  is  left  unchanged,  however, 
on  grounds  that  general  deductions  are  still  pertinent. 
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DEPTH  BELOW  MEAN  LOW  WATER,  tf  -  (faat) 


DISTANCE  FROM  BAY- MOUTH  -  ( navlloal  milM I 
It  B  10  I  B  10*  t  B  lo’ 


FIGURE  1-6 

PROFILES  OF  SEA-BED  ALONG  GREAT  CIRCLE  PATHS 
FROM  MOUTHS  OF  CHESAPEAKE  AND  NEW  YORK  BAYS 
TO  OCEAN  STATION  AT  LAT.  35°  30',  LONG.  65°  27' 
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Applicable  values  of  n,  0^,  m,  now  readily 

evaluated  from  the  design  curves  of  Figs.  1-2  and  1-5  are  recorded 
in  Table  I- 1 . 


TABLE  I-  1 

DESIGN  CONSTANTS  GOVERNING  WAVE  HEIGHT  AND  PERIOD 


Relative  Distance 

Wave  Height  Decay 

Wave  Period  Increase 

r/d 

n 

«n 

m 

<35 

1 

0.  127 

1/6 

2.90 

>35 

5/6 

0. 0734 

1/4 

2.  12 

Eqs.  (1-39)  and  (1-43)  have  been  used  in  conjunction  with  the 
parameters  of  Table  I- 1  in  a  stepwise  numerical  integration  pro¬ 
cedure  for  computing  and  T  over  the  (full-line)  depth  profile 

shown  in  Fig.  1-6,  carried  as  far  as  the  foot  of  the  continental 
slope.  The  results  of  the  calculations  are  portrayed  in  Fig.  1-7. 

The  highest  wave,  without  regard  to  its  position  in  the  wave 
groups,  will  decline  from  an  initial  elevation  7f^  =  388  ft.  at 
r  =  19  n.  mi.  to  17^  =  17.3  ft.  at  r  =  428  n.  mi.  from  the  origin. 
Over  this  distance  the  water  depth  decreases  from  15,  420  ft.  to 
7800  ft.  Maximum  water  depth  of  17,  760  ft.  is  encountered  at  the 
intermediate  distance  r  -  155.5  n.  mi.  The  further  history  of 
wave  height  has  not  been  pursued  beyond  r  s  428  n.  mi.  because 
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of  the  complications  of  reflection  at  the  continental  slope,  which 
require  separate  study  (as  given  in  Chapter  IV,  Part  B,  Vol.  I). 

Over  the  same  wave  path  the  period  of  the  leading  wave, 

=  165  secs,  at  r  =  19  n.  mi.  is  found  to  increase  to 
Tj  =  192  secs,  at  r  =  428  n.  mi.  On  the  assumption  that  the 
fourth  wave  in  the  leading  group  of  waves  follows  the  same  laws 
of  period  increase  as  the  actual  wave  front,  its  period  of 

T.  =  72  secs,  at  r  =  19  n.  mi.  is  found  to  increase  to  81  secs. 

4 

at  r  =  428  n.  mi.  The  legitimacy  of  this  calculation  may  be  open 
to  question  as  the  governing  law  may  be  more  in  accord  with 
Eq.  (1-26).  For  the  initial  values  of  d  =  15,500  ft.  at  r  =  19  n.  mi. 
the  latter  equation  gives  =  70.  5  secs,  for  m  =  3  (the  fourth 
wave)  which  is  in  surprising  agreement  with  =  72  secs,  of  the 
given  data.  Since  Eq.  (1-26)  is  valid  only  for  constant  d,  we 
must  take  a  mean  value  of  depth  over  the  distance  to  estimate 
at  r  =  428  n.  mi.  For  d  =  15,000  ft.  then,  is  found  to  be 

197  secs,  at  r  =  428  n.  mi.  on  this  basis.  This  value  exceeds  the 
computed  value  of  T^  =  192  secs,  and  is  therefore  not  likely  to  be 
correct.  The  inference  must  be  that  Eq.  (1-26)  is  only  applicable 
much  further  back  in  the  body  of  the  waves  for  larger  values  of  m. 

8.  CONCLUSIONS 

On  the  basis  of  the  arguments  evolved  in  this  appendix,  and 
on  the  assumption  that  a  nuclear  explosion  with  a  bubble  crater 
diameter  =  20,  800  ft.  produces  maximum  waves  of  a  height 
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7)  =  388  ft.  above  still  water  at  a  distance  of  r  =  20  miles  from 
'o 

surface  zero  in  a  water  depth  of  15,  420  ft. ,  it  is  found  that  over 

the  shelving  depth  to  the  base  of  the  continental  shelf  (d  fit  7800  ft.) 

the  maximum  waves  will  have  decreased  to  n  =  17.3  ft.  at 

'  o 

r  s  428  n.  mi.  The  period  of  the  leading  wave  (which  will  be  of 
much  lower  height)  will  approximate  192  secs.  (3.2  mins.),  but 
the  period  of  the  highest  wave  will  very  likely  be  in  the  neighbor¬ 
hood  of  80  secs,  at  this  distance. 

At  a  radius  of  300  n.  mi.  from  surface  zero  (d  ^  10,  000  ft. ) 

Fig.  1-7  shows  that  the  probable  maximum  height  above  still  water 
of  the  highest  wave  will  be  about  =  25  ft.  and  its  corresponding 

period  about  T  s  78  secs. ,  which  is  in  very  good  agreement  with 
what  has  been  found  in  Appendix  11. 

A  question  as  yet  unanswered  is  whether  the  Kranzer- Keller 
theory  and  the  asymptotic  solution  which  it  invokes  is  legitimately 
applicable  to  a  system  of  high  waves  which,  in  the  initial  stages  at 
least,  are  likely  to  be  extremely  non-linear.  However,  an  encouraging 
aspect  of  this  query  is  that  experimental  results  of  wave  height  and 
period- change,  even  for  comparatively  small  values  of  r/d,  do  con¬ 
firm  the  theoretical  predictions  based  on  the  asymptotic  solutions  of 
the  linearized  theory.  Since  the  Kranzer-Keller  theory  is  merely 
an  adaptation  of  the  more  fundamental  Cauchy- Poisson- Lamb  linear 
theory  which  conceives  the  emergent  waves  as  the  sum  of  an  infinite 
number  of  small  amplitude  waves  distributed  over  a  continuous  spec¬ 
trum  of  frequencies,  there  is  the  possibility  that  the  aggregate 
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outcome  of  its  legitimate  linear  processes,  applied  to  the  spectral 
elements,  is  a  very  close  approximation  (if  not  an  exact  solution) 
to  the  characteristics  of  large  waves,  which  as  entities  must  other¬ 
wise  be  considered  as  obeying  non-linear  laws. 
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LIST  OF  SYMBOLS 
FCF  APPENDIX  I 


c  phase  velocity  of  water  waves 

d  water  depth 

diameter  of  bubble  crater  at  the  water  surface  in 
underwater  explosions 

e  universal  constant  (2.  718.  .  . ) 

g  acceleration  due  to  gravity 

h  depth  below  the  surface  of  an  underwater  explosion 

H  wave  height,  crest  to  trough 

I  uniform  impulse  (per  unit  area)  over  source  area 

I  (r)  function  describing  the  initial  impulse  as  a  function 

of  r 

T  (k)  Hankel  transform  of  I  (r)  (Eqs.  (I- 11)  ) 

concentrated  impulse  at  the  point  source  (per  unit  area) 

J  Bessel  function  of  zero  order 

o 

Bessel  function  of  first  order 

k  wave  number  (  =  ^  ^  ) 

A 

K  constant  of  proportionality,  function  of  R/d,  (Eqs. 

(1-36  ii) 

m  (1)  numerical  exponent;  (2)  corresponding  subscript 

n  (1)  numerical  exponent;  (2)  corresponding  subscript 

p  integer  number  representing  the  wave  with  greatest 

height  counted  from  the  front  of  the  train 

Q  vertical  surface  elevation  over  a  circle  of  radius  R 

Q  (r)  initial  surface  elevation  as  a  function  of  r 
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Hankel  transform  of  Q  (r),  (Eqs.  I- 11)  ) 

initial  surface  elevation  of  point  source  disturbance 
at  the  origin  (per  unit  area) 

variable  horizontal  radial  distance  in  cylindrical 
symmetry 

radius  of  a  cylindrically  symmetrical  surface  or 
bottom  disturbance 

variable  time 

wave  period 

period  of  the  initial  or  leading  wave  in  a  dispersive 
system 

group  velocity  of  water  waves 

charge-weight  of  nuclear  explosion 

that  part  of  the  total  energy  W  which  is  responsible 
for  the  generation  of  waves 

variable  horizontal  distance 

dimensionless  variable,  (Eq.  (1-43  ii)  ) 

dimensionless  variable  (Eq.  1-39  i)  ) 

differential  coefficient  of  Y  with  respect  to  X 
(Eq,  1-40  ii)  ) 

dimensionless  variable  (Eq.  (1-43  i)  ) 

differential  coefficient  of  Z  with  respect  to  X, 

(Eq.  (1-44  ii)  ) 

constant  of  proportionality,  (Eq.  (1-38)  ) 

constant  of  proportionality,  (Eq.  (1-42)  ) 

elevation  of  wave  surface  above  still  water  level 

maximum  value  of  at  the  wave  crest 

value  of  7j  deriving  from  an  initial  elevation  over 
a  limited  source  region 


value  of  rj  deriving  from  an  initial  impulse  applied 
over  a  limited  source  region 

wave  length  of  a  wave  in  water  of  any  depth 

deep-water  wave  length  of  any  given  wave  type 

wave  length  of  the  leading  wave  in  a  dispersive  system 

coefficient  of  virtual  viscosity  (sec.  units),  (Eq, 
(1-5)  ) 

universal  constant  (3. 14159. . . ) 
mass  density  of  (sea)  water 
angular  frequency  (  =  2  ir/T) 

function  of  (kd)  in  Kranzer- Keller  theory,  common  to 
cases  of  initial  elevation  and  initial  impulse 

function  of  kd  in  Kranzer- Keller  theory,  case  of 
initial  elevation 

function  of  kd  in  Kranzer- Keller  theory,  case  of 
initial  impulse 

function  of  kR  and  kd  (£q.  (1-28)  ) 
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THEORETICAL  CONSIDERATIONS  AND  COMPUTATIONS 
FOR  WATER  WAVES  PRODUCED  BY  EXPLOSIONS 


By 

Larry  Armijo 
and 


Mary  Ann  Noser 


1 .  INTRODUCTION 


A  critical  literature  survey  on  the  various  theories  for 
studying  gravity  waves  generated  by  underwater  explosion  has 
been  given  in  Appendix  I.  It  is  recalled  that  a  description  of  the 
motion  of  the  free  surface  of  a  body  of  water  of  constant  depth 
following  an  arbitrary  initial  distribution  of  impulse  applied  to  the 
surface,  or  an  initial  elevation  or  depression  of  the  surface  of 
arbitrary  shape,  is  given  by  the  Kranzer- Keller  theory  (1955)  for 
the  radially  symmetric  case.  Kaplan,  Wallace  and  Goodale  (1962) 
have  carried  out  computations  of  the  Kranzer- Keller  equations  for 
three  cases  of  initial  paraboloidal  depressions  of  the  surface  for 
a  point  on  the  surface  twenty  miles  from  the  center  of  the  initial 
disturbance.  These  results  have  been  presented  in  Fig.  1  (Vol.  I) 
of  this  report.  The  primary  object  of  this  appendix  is  to  present 
a  similar  set  of  computations  for  the  same  three  cases  considered 
by  Kaplan,  Wallace  and  Goodale,  but  at  a  point  on  the  surface  300 
miles  from  the  center  of  the  initial  disturbance.  For  the  cases 
under  consideration,  certain  theoretical  results  have  been  obtained 
by  the  authors  and  these  results  are  also  presented. 

2.  PRINCIPAL  EQUATIONS 

The  displacement  t)  (r,  t)  of  the  free  surface  is  given  by 

rf{T,\)  »  ^  B  COS  [zir  (II- 1) 


II- 1 


where 


B  » 


*\tr) 


i 


and  cr 

X 

T 

h 

R 

r 

Vo 

t 


_  2.tr  h 
X 

=  wave  length 
=  wave  period 
s  fluid  depth 

=  effective  radius  of  disturbance 
=  distance  from  center  of  disturbance 
=  depth  of  disturbance  at  r  =  0 
=  time 


ai-2) 


The  function  is  the  Hankel  transform  of  the  function  £(r), 

the  initial  paraboloidal  depression,  which  is  assumed  to  be  of  the 
form 


where 


E(r) 


Ar*  -  170  ,  0  <  r  <  y, 

0,  Y  -  ^  ^ 


A 


J2_ 

yz 


See  Figure  II- 2. 

It  follows  that 

E(8) 


2^0 

s< 


Jg  (ys) 


(11-3) 


(II- 4) 


(II- 5) 


II- 3 


11-4 


NOTATION  FOR  THE  ORICaN 


where  Bessel  function  order  2. 

The  function-  ^(o*)  is  given  by 


^(<r) 


X 


2 


I 

2(cosh  a- 


<T 

sinh  O’ 


O’  >  0 


(II- 6) 


The  function  ^(<r)  is  a  positive  monotone  decreasing  function  of 

O’  and  lim  ^(O’)  =  1.  i.e.  ^o)  =  1. 

O’—^o 

The  parameter  o’  is  related  to  the  quantities  r  and  t  by 
the  equation 


VgTT  I  ■  '  ■  V5h 
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so  that 


t).  ’  - 


(II- 8) 


where  ^  ^  is  the  inverse  of  ^  and  g  is  the  gravitational  ac¬ 
celeration  constant. 

The  period  T  is  given  by  the  equation 


2»r 


O’  tonh  cr 


(II- 9) 


The  effective  radius  R  is  related  to  the  quantity  y  by  the  equation 


R  » 


yr 


(II- 10) 


It  follows  from  (II-2),  (II-5)  and  (11-10)  that 
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2yr  h 


y<r* 


(■^) 


f  (0-) 


(11-11) 


where 


e{o-) 


s  f-o- 

L 


(11-12) 


For  convenience,  we  have  defined 


B*(r,t)  =  B  =  -  ^^0;  Jg  (11-13) 

'  r<r*  \  n  / 


and 


«*(M)S  t>^ 


ysT 


(11-14) 


so  that  by  (II- 1),  (11-13)  and  (11-14) 


17  (r,t)  =  B*{r,t)  cose*(r,t) 


(11-15) 


The  quantity  B  (r,  t)  therefore  represents  the  displacement 
of  the  amplitude  of  the  wave  envelope.  The  quantity  0  (r,  t)  given 
by  (11-14)  may  be  shown  to  be  a  monotone  increasing  function  of 


t  for  t  > 


Moreover,  lim  Q  (r,t)  ®  0  and  Mm  9  (r,t)  *  00 


t/x) 


11-6 


The  times  of  arrival  of  crests  and  troughs  are  given  by  the  solu¬ 
tions  t  of  the  equations 
xn 

(r.tn)  -  nir,  n  -  1,2...  (11-16) 


3.  APPROXIMATIONS 

Following  the  example  set  by  Kaplan,  Wallace  and  Goodale, 

we  have  replaced  the  functions  ^(a)  and  ((c)  in  the  above 

equations  by  their  asymptotic  expressions  - and  ^/2(r  ’ 

2  va 

respectively,  for  values  o*  >  4  .  The  parameter  o-  may  be 
eliminated  from  the  preceding  equations  and  the  following  simplifi¬ 
cations  occur: 


B*  (r.O  « 


8^/2  rjo  r 


(11-17) 


and 


9  (r.t) 


4Trr 
Q  t  ’ 


X 


Str  fg 
gt* 


t  > 


(11-18) 
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The  times  of  arrival  t^  of  crests  and  troughs  may  be  computed 
from  the  equations _ _ _ _ _  '  " 


t 


n 


m 


Trrn 

g 


for  positive  integers 


1.  >  4k  <"■>’> 


The  following  theoretical  results  have  been  noted  by  the 

authors  of  this  report.  The  aeros  of  B  (r,  t)  occur  at  the  times 

e 

t  given  by 


U-7 


(11-20) 


where  x^* 

If  -^  > 
the  times  t 


js  the  positive  zero  of  J2(x). 

1.739.  the  maximum  values  of  |  B*(r,  t)  | 
given  by 


occur  at 


t 


» 

k 


(11-21) 


where  x^*  is  the  positive  zero  of  J2(x)  -  x  Jj(x). 

To  compute  the  first  two  packets  of  waves,  we  have  made  use 
of  the  values 


k 

- 5 - 

*k 

’‘k 

1 

5.  136 

2.  300 

2 

8.417 

6.541 

The  restriction 


>  1.739,  is  satisfied  for  Cases  I  and 


II  considered  by  Kaplan,  Wallace  and  Goodale,  For  their  Case  III, 
y  =1.731  and  the  restriction  is  very  nearly  satisfied.  Eq.  (11-20) 
follows  from  (11-13),  and  Eq.  (11-21)  follows  from  (11-17),  which 
follows  the  asymptotic  expressions  for  ^(cr)  and  ((o') 


It  should  be  emphasized  that  the  value  t 


«»< 

1 


'4r^x 


2w*« 


is  the  time  of  arrival  of  the  maximum  value  of  the 


r  9 


absolute  value  of  the  amplitude  of  the  wave  envelope.  This  maximum 
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B  (r,  }|  .  The  value  will  be  very 

nearly  equal  to  the  time  of  arrival  of  the  highest  crest  and  the 
value  iB^ir.  tj'*"'')|  will  be  very  nearly  equal  to  the  displace¬ 
ment  of  the  highest  crest. 

4.  DESCRIPTION  OF  COMPUTATIONS 

In  calculating  the  displacement  1)  (r,  t),  it  is  first  necessary 

to  compute  the  functions  and  ((cr)  A  table  of  these 

functions  has  been  prepared,  using  a  Bendix  G-15  computer,  with 

O'  incremented  by  0.01  in  the  interval  from  0  to  4,  in  order  to 

simplify  any  future  hand  computations.  These  results  are  found  in 

Table  11-4.  Also,  the  following  numerical  values  were  used  to 

a 

compute  ^  (r,  t)  and  B  (r,  t): 

300  miles  »  1,  584,  000  feet 
18,  000  feet 

1.  1592  X  10^  ft.  /min.  ^ 

These  computations  were  carried  out  for  the  following  three  cases 


r  ■ 

h  s 

g 


Dimensions  of  Initial  Displacement 

Depth 

Radius 

(feet) 

(feet) 

Case  I 

3,  380 

4,  440 

Case  II 

6,  010 

7,900 

Case  III 

7,920 

10,  400 

TABLE  n-  1 
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The  time  at  which  disturbances  begin  at  a  distance  r 
from  the  center  of  the  initial  depression  is  the  same  for  all  three 
cases  and  is  given  by 


t 


0 


34.67  min. 


■S'  tjt 

The  times  of  arrival  t^  and  t^  of  the  first  two  zeros  of  the  dis- 

»{* 

placement  of  the  wave  envelope  B  (r,  t)  were  computed  from  Eq. 
(11-20)  and  are  listed  in  Table  II-2  ■ 


Times  of  Arrival  of  the  First  Two  Zeros  of 

the  Displacement  of  the  Wave  Envelope 

k 

tj^  (min.  ) 

Case  I 

1 

2 

Case  II 

1 

237.2 

2 

303.  6 

Case  III 

1 

206.7 

2 

264.  6 

TABLE  II- 2 

The  times  of  arrival  t^  and  t^  of  the  first  two  maxima 
of  the  absolute  value  of  the  displacement  of  the  wave  envelope 
I  B  (r,  t)  I  as  well  as  the  maximum  values  I  B  (r,  tj  )|  , 

B  (r,  t^  )l  .  and  the  periods  and  wave  lengths  at  times  tj 
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and  t2  were  computed  from  Eqe.  (11-21)  and  (U- 17),  and  are 
listed  below. 


Times  of  Arrival  of  the  First  Two  Maxima  of  the 

Absolute  Value  of  the  Displacement  of  the  Wave 

Envelope  and  the  Corresponding  Maximum  Values, 

Periods,  and  Wave  Lengths 

Case 

k 

B*(r,  t^*) 

T 

X 

(min) 

(feet) 

(min) 

(feet) 

1 

1 

0.811 

2 

HEEBli 

0.481 

11 

1 

15.268 

1.082 

2 

4.018 

0.  642 

U1 

1 

138.3(7) 

26,  488(7) 

1,242 

28, 440 

2 

233.3 

6.971 

0.736 

10,010 

TABLE  11-3 


The  values  followed  by  ( ?)  are  doubtful  because  the  restriction 
>  1 . 739  is  not  satisfied. 

In  the  range  34.67  <  t  <  138.01,  which  corresponds  to 

the  range  0  <  O'  <  4,  the  values  of  the  displacement  of  the  wave 

envelope  B*(r,  t)  were  computed  from  Eq.  (11-13)  using  the  Bendix 

G-15  computer.  For  values  138.01  <  t  <  t^  ,  the  values  of 

B  (r,  t)  were  computed  by  hand  using  Eq.  (11-17). 

*  The  tables  of  the  Bessel  Functions  of  the  First  Kind  of  Orders  2 
and  3  by  the  Staff  of  the  Computation  Laboratory,  published  in  Cam¬ 
bridge  by  the  Harvard  Press,  were  used  in  these  computations. 
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TABLE  n>4 


THE  FUNCTIONS  ^(<r)AND  (  («r) 


or 

€  («■) 

0,01 

.9999500 

1.0000277 

0.02 

.9998000 

1. 0001111 

0.03 

.9995502 

1.0002500 

0.04 

.9992006 

1.0004444 

0.05 

.9907516 

1.0006944 

0.06 

.9982034 

1.0010000 

0.07 

.9975563 

1.0013612 

0.08 

.9968107 

1.0017779 

0.09 

.9959672 

1.0022503 

O.IO 

•9950262 

1.0027782 

0.11 

.9939884 

1.0033617 

0.12 

•9928543 

1.0040009 

0.13 

•9916247 

1.0046957 

0.14 

.9903004 

1.0054462 

0.15 

.9888^21 

1.0062523 

0.16 

.9873708 

1.0071141 

0.17 

.9857673 

1.0080317 

0.18 

.9840727 

1.0090049 

0.19 

•9822680 

1.0100339 

0.20 

.9804142 

1.0111186 

0.21 

.9784525 

1.0122592 

0.22 

.9764040 

1.0134556 

0.23 

.9742700 

1.0147078 

0.24 

.9720518 

1.0160158 

0.25 

.9697506 

1.0173798 

0.26 

.9673677 

1.0187997 

0.27 

.9649046 

1.0202756 

0.28 

.9623627 

1.0218075 

0.29 

.9597435 

1.0233954 

0.30 

.9570483 

1.0250395 

0.31 

.9542788 

1.0267396 

0.32 

.9514365 

1.0284959 

0.33 

.9485229 

1.0303084 

0.34 

.9455397 

1.0321771 

0.35 

.9424884 

1.0341022 

0.36 

.9393708 

1.0360836 

0.37 

.9361884 

1.0381214 

0.38 

.9329430 

1.0402157 

0.39 

.9296362 

1.0423665 

0.40 

.9262699 

1.0445739 
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TABLE  II-4  (continued) 

9 

^  (o) 

C(») 

0.41 

.9228496 

1.0468379 

0.4t 

.9193691 

1.0491986 

0.43 

.9198302 

1.0919361 

0.44 

.9122427 

1.0539704 

0.4S 

.9086042 

1.0564617 

0.46 

.9049169 

1.0590099 

0.47 

.9011819 

1.0616191 

0.40 

.8974008 

1.0642779 

0.49 

.8939762 

1.0669972 

0.90 

.8897099 

1.0697741 

0.91 

.8898000 

1.0726000 

0.92 

.8818600 

1.0759000 

0.93 

.8778800 

1.0784000 

0.94 

.8738600 

1.0815000 

0.99 

.8698100 

1.0845000 

0.96 

.8657200 

1.0876000 

0.97 

.8616100 

1.0908000 

0.90 

.8974700 

1.0941000 

0.99 

.8533000 

1.0974000 

0.60 

.8491100 

1.1007000 

0.61 

.8448900 

1.1041000 

0.62 

.8406500 

1.1076000 

0.63 

.8363900 

1.1111000 

0.64 

.8321200 

1.1147000 

0.69 

.8278200 

1.1184000 

0  •  66 

.8235100 

1.1221000 

0.67 

.8191900 

1.1258000 

0.68 

.8148500 

1.1297000 

0.69 

.8105100 

1.1336000 

0.70 

.8061500 

1.1375000 

0.71 

.8017900 

1.1415000 

0.72 

.7974100 

1.1456000 

0.73 

.7930300 

1.1497000 

0.74 

.7886500 

1.1539000 

0.79 

•7842700 

1.1581000 

0.76 

.7798800 

1.1624000 

0.77 

.7755000 

1.1668000 

0.78 

.7711100 

1.1712000 

0.79 

.7667300 

1.1757000 

0.80 

.7623900 

1.1803000 

U-13 


TABLE  11-4  (continued) 
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i  (<r) 

0.81 

.7579700 

1.1849000 

0.82 

•7936000 

1.1899000 

0.8) 

•7492400 

1.1943000 

0.84 

.7448800 

1.1991000 

0.85 

•7409300 

1.2039000 

0.86 

.7362000 

1.2088000 

0.87 

•7318700 

1.2138000 

0.88 

.7275500 

1.2189000 

0.89 

.7232900 

1.2240000 

0.90 

•7189600 

1.2291000 

0.91 

•7146900 

1.2344000 

0.92 

.7104200 

1 .2397000 

0.93 

•7061600 

1.2450000 

0.94 

•7019900 

1.2504000 

0.99 

•6977400 

1.2559000 

0.96 

•6935500 

1.2615000 

0.97 

.6893700 

1.2671000 

0.98 

.6852200 

1.2728000 

0.99 

.6810800 

1.2785000 

1.00 

•6769700 

1.2643000 

1.01 

.6728700 

1.2902000 

1.02 

.6688000 

1.2961000 

1.03 

•6647500 

1.3021000 

1.04 

•6607300 

1.3082000 

1.09 

•6967200 

1.3144000 

1.06 

•6527400 

1.3206000 

1.07 

.6467900 

1.3268000 

1.08 

.6448600 

1.3332000 

1.09 

.6409500 

1.3396000 

1.10 

.6370700 

1.3460000 

1.11 

.6332100 

1.3526000 

1.12 

.6293900 

1.3592000 

1.13 

.6255800 

1.3658000 

1.14 

•6218100 

1.3726000 

1.19 

•6160600 

1.3794000 

1.16 

•6143400 

1.3863000 

1.17 

•6106400 

1.3932000 

1.18 

.6069700 

1.4002000 

1.19 

•6033300 

1.4073000 

1.20 

•9997200 

1.4144000 

TABLE  U-4  (continued) 


9 

^  {9) 

f  (ff) 

1.21 

.5961400 

1.4217000 

1.22 

.5925900 

1.4290000 

1.29 

.5890600 

1.4363000 

1.24 

.5855600 

1.4497000 

1.2S 

.5820900 

1.4512000 

1.26 

.5786500 

1.4588000 

1.27 

.5752400 

1.4664000 

1.28 

.5718600 

1.4741000 

1.29 

.5685000 

1.4819000 

1.30 

.5651800 

1.4898000 

1.31 

.5618800 

1.4977000 

1.32 

.5586200 

1.5057000 

1.33 

.5553800 

1.5138000 

1.34 

.5521700 

1.5219000 

1.35 

.5489900 

1.5301000 

1.36 

.5458400 

1.5384000 

1.37 

.5427200 

1.5467000 

1.38 

.5396200 

1.5552000 

1.39 

.5365600 

1.5637000 

1.40 

.5335200 

1.5722000 

1.41 

.5305100 

1.5809000 

1.42 

.5275400 

1.5896000 

1.43 

.5245900 

1.5984000 

1.44 

.5216600 

1.6072000 

1.45 

.5187700 

1.6162000 

1.46 

.5159000 

1.6252000 

1.47 

.5190600 

1.6343000 

1.48 

.5102600 

1.6434000 

1.49 

.5074700 

1.6526000 

1.50 

.5047200 

1.6619000 

1.51 

.5019900 

1.6713000 

1.52 

.4992900 

1.6808000 

1.53 

.4966200 

1.6903000 

1.54 

.4939700 

1.6999000 

1.55 

.4913500 

1.7096000 

1.56 

.4887600 

1.7193000 

1.57 

.4861900 

1.7291000 

1.58 

.4836500 

1.7390000 

1.59 

.4811400 

1.7490000 

1.60 

.4786500 

1.7591000 
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TABLE  11-4  (continued) 


er 

^  (©•) 

i  (tr) 

1.61 

.4761900 

1.7692000 

1.62 

.4737500 

1.7794000 

1.63 

.4713400 

1.7896000 

1.04 

.4689600 

1.8COOOOO 

1.65 

.4665900 

1.8104000 

1.66 

.4642600 

1.8209000 

1.67 

.4619500 

1.8315000 

1.68 

.4596600 

1.8421000 

1.69 

.4573900 

1.8528000 

1.70 

.4551500 

1.8636000 

1.71 

.4529400 

1.8745000 

1.72 

.4507400 

1.8854000 

1.73 

.4485600 

1.8964000 

1.74 

.4464300 

1.9075000 

1.75 

.4443100 

1.9187000 

1.76 

.4422000 

1.9299000 

1.77 

.4401300 

1.9412000 

1.78 

.4380700 

1.9526000 

1.79 

.4360400 

1.9641000 

1.80 

,4340200 

1.9756000 

1.81 

,4320300 

1.9872000 

1.82 

.43C0600 

1.9989000 

1.83 

.4281100 

2.0106000 

1.84 

.4261900 

2.0225000 

1.85 

.4242600 

2.0344000 

1.86 

.4223900 

2.0463000 

1.87 

.4205300 

2.0584000 

1.88 

.4186800 

2.0705000 

1.89 

.4168600 

2.0827000 

1.90 

.4150500 

2.0949000 

1.91 

.4132600 

2.1073000 

1.92 

.4114900 

2.1197000 

1.93 

.4097400 

2.1322000 

1.94 

.4080100 

2.1447000 

1.95 

.4063000 

2.1573000 

1.96 

.4046100 

2.1700000 

1.97 

.4029400 

2.1828000 

1.98 

.4012800 

2.1956000 

1.99 

.3996400 

2.2085000 

2.00 

.3980200 

2.2214000 
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TABLE  U-4  (continued) 


9 

4  (o) 

€  (9) 

2.01 

•3964100 

2.2345000 

2.02 

.3948300 

2.2476000 

2.03 

.3932500 

2.2608000 

2.04 

.3917000 

2.2740000 

2.05 

•3901600 

2.2873000 

2.06 

.3886400 

2.3007000 

2.07 

•3871400 

2.3141000 

2.06 

.3856500 

2.3276000 

2.09 

•3641800 

2.3412000 

2.10 

.3827200 

2.3548000 

2.11 

.3812800 

2.3685000 

2.12 

.3798500 

2.3823000 

2.13 

.3784400 

2.3961000 

2.14 

.3770400 

2.4100000 

2.15 

.3756600 

2.4240000 

2.16 

.3742900 

2.4380000 

2.17 

.3729400 

2.4521000 

2.18 

•3716000 

2.4663000 

2.19 

•3702800 

2.4805000 

2.20 

•3689700 

2.4947000 

2.21 

.3676700 

2.5091000 

2.22 

•3663900 

2.5235000 

2.23 

•3651100 

2.5379000 

2.24 

•3638600 

2.5524000 

2.25 

.3626100 

2.5670000 

2.26 

•3613B00 

2.5816000 

2.27 

•3601600 

2.5963000 

2.28 

.3589600 

2.6111000 

2.29 

,3577600 

2.6259000 

2.30 

.3565800 

2.6407000 

2.31 

•3554100 

2.6556000 

2.32 

.3542500 

2.6706000 

2.33 

•3531000 

2.6856000 

2.34 

.3519700 

2.7007000 

2.35 

•3508500 

2.7158000 

2.36 

.3497300 

2.7310000 

2.37 

.3486300 

2.7462000 

2.36 

•3475400 

2.7615000 

2.39 

•3464600 

2.7768000 

2.40 

.3453900 

2.7922000 
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TABLE  11-4  (continued) 


<r 

^  (o’) 

{  (o’) 

2.41 

.3443400 

2.8076000 

2.42 

.3432900 

2.8231000 

2.43 

.3422500 

2.8386000 

2.44 

.3412200 

2.8542000 

2.45 

.3402100 

2.8698000 

2.46 

.3392000 

2.8855000 

2.47 

.3382000 

2.9012000 

2.48 

.3372200 

2.9170000 

2.49 

.3362400 

2.9328000 

2.50 

.3352700 

2.9486000 

2.51 

.3343100 

2.9645000 

2.52 

.3333600 

2.9804000 

2.53 

.3324200 

2.9964000 

2.54 

.3314900 

3.0124000 

2.55 

.3305600 

3.0285000 

2.56 

.3296500 

3.0446000 

2.57 

.3287400 

3.0607000 

2.58 

.3278400 

3.0769000 

2.59 

.3269500 

3.0931000 

2.60 

.3260700 

3.1093000 

2.61 

.3252000 

3.1256000 

2.62 

.3243400 

3.1419000 

2.63 

.3234800 

3.1583000 

2.64 

.3226300 

3.1746000 

2.65 

.3217900 

3.1910000 

2.66 

.3209500 

3.2075000 

2.67 

.3201300 

3.2240000 

2.68 

.3193100 

3.2405000 

2.69 

.3185000 

3.2570000 

2.70 

.3177000 

3.2736000 

2.71 

.3169000 

3.2902000 

2.72 

.3161100 

3.3068000 

2.73 

.3153300 

3.3235000 

2.74 

•3145500 

3.3401000 

2.75 

•3137800 

3.3568000 

2.76 

•3130200 

3.3736000 

2.77 

•3122600 

3.3903000 

2.78 

.3115100 

3.4071000 

2.79 

.3107700 

3.4239000 

2.80 

.3100300 

3.4407000 
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TABLE  II-4  (continued) 


cr 

^  (v) 

(  (cr) 

2.81 

.3093000 

3.4576000 

2.82 

•3085600 

3.4745000 

2.83 

•3078600 

3.4913000 

2.84 

.3071500 

3.5082000 

2.85 

.3064400 

3.5252000 

2.86 

.3057400 

3.5421000 

2.87 

.3050500 

3.5591000 

2.88 

.3043600 

3.5761000 

2.89 

.3036800 

3.5931000 

2.90 

.3030000 

3.6101000 

2.91 

.3023300 

3.6271000 

2.92 

.3016700 

3.6442000 

2.93 

.3010000 

3.6612000 

2.94 

.3003500 

3.6783000 

2.95 

.2997000 

3.6954000 

2.96 

.2990500 

3.7125000 

2.97 

.2984100 

3.7296000 

2.98 

.2977800 

3. 7467000 

2.99 

.2971500 

3.7638000 

3.00 

•2965300 

3.7810000 

3.01 

•2959100 

3.7981000 

3.02 

•2952900 

3.8153000 

3.03 

.2946600 

3.8324000 

3.04 

.2940700 

3.6496000 

3.05 

•2934700 

3.8668000 

3.06 

•2928600 

3. 6640000 

3.07 

•2922800 

3.9012000 

3.08 

•2917000 

3.9184000 

3.09 

.2911200 

3.9356000 

3.10 

•2905400 

3.9526000 

3.11 

•2899600 

3.9700000 

3.12 

.2893900 

3.9872000 

3.13 

.2888300 

4.0044000 

3.14 

.2882600 

4.0216000 

3.15 

.2877100 

4.0386000 

3.16 

.2871500 

4.0560000 

3.17 

•2866000 

4.0733000 

3.18 

•2860600 

4.0905000 

3.19 

•2855200 

4.1077000 

3.20 

.2849800 

4.1249000 
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TABLE  II-4  (continued) 


9 

^  (o’) 

^  (o’) 

3.21 

.2844400 

4.1421000 

3.22 

.2839100 

4.1594000 

3.23 

.2833900 

4.1766000 

3.24 

.2828700 

4.193B000 

3.25 

.2823500 

4.2110000 

3.26 

.2818300 

4.2282000 

3.27 

.2813200 

4.2454000 

3.28 

.2808100 

4.2626000 

3.29 

.2803000 

4.2798000 

3.30 

.2798000 

4.2969000 

3.31 

.2793000 

4.3141000 

3.32 

.2788100 

4.3313000 

3.33 

.2783200 

4.3484000 

3.34 

.2778300 

4.3656000 

3.35 

.2773400 

4.3828000 

3.36 

.2768600 

4.3999000 

3.37 

.2763800 

4.4170000 

3.38 

.2759100 

4.4342000 

3.39 

.2754300 

4.4513000 

3.40 

.2749600 

4.4684000 

3.41 

.2745000 

4.4855000 

3.42 

.2740300 

4.5026000 

3.43 

.2735700 

4.5197000 

3.44 

.2731100 

4.5367000 

3.45 

.2726600 

4.5538000 

3.46 

.2722100 

4.5708000 

3.47 

.2717600 

4.5879000 

3.48 

.2713100 

4.604900C 

3.49 

.2708700 

4.6219000 

3.50 

.2704300 

4.6389000 

3.51 

.2699900 

4.6559000 

3.52 

.2695500 

4.6729000 

3.53 

.2691200 

4.6899000 

3.54 

.2686900 

4.7068000 

3.55 

.2682600 

4.7238000 

3.56 

•2678300 

4. 7407000 

3.57 

.2674100 

4.7576000 

3.58 

•2669900 

4.7745000 

3.59 

.2665700 

4.7914000 

3.60 

.2661600 

4.8082000 
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TABLE  11>4  (continued) 


9 

^  (er) 

<  (»•) 

9.61 

.2657400 

4.8251000 

3.62 

.2653300 

4.8420000 

3.63 

.2649200 

4.8588000 

3.64 

.2645200 

4. 8756000 

3.65 

.2641100 

4.8924000 

3.66 

.2637100 

4.9092000 

3.67 

.2633100 

4.9260000 

3.68 

.2629200 

4.9427000 

3.69 

.2625200 

4.9595000 

3.70 

.2621300 

4.9762000 

3.71 

.2617400 

4.9929000 

3.72 

.2613500 

5.0096000 

3.73 

.2609600 

5.0263000 

3.74 

.2605800 

5.0429000 

3.75 

.2602000 

5.0596000 

3.76 

.2598200 

5.0762000 

3.77 

.2594400 

5.0928000 

3.78 

.2590600 

5.1094000 

3.79 

.2586900 

5.1260000 

3.80 

.2583200 

5.1425000 

3.81 

.2579500 

5.1591000 

3.82 

.2575800 

5.1756000 

3.83 

.2572100 

5.1921000 

3.84 

.2568500 

5.2086000 

3.85 

.2564800 

5.2251000 

3.86 

.2561200 

5.2416000 

3.87 

.2557600 

5.2580000 

3.88 

.2554100 

5.2745000 

3.89 

.2550500 

5.2909000 

3.90 

.2547000 

5.3073000 

3.91 

.2543500 

5.3236000 

3.92 

.2540000 

5.3400000 

3.93 

.2536500 

5.3563000 

3.94 

.2533000 

5.3727000 

3.95 

.2529600 

5.3890000 

3.96 

.2526100 

5.4053000 

3.97 

.2522700 

5.4216000 

3.98 

.2519300 

5.4378000 

3.99 

.2515900 

5.4541000 

4.00 

.2512600 

5.4703000 

These  results  are  displayed  in  Figs.  11-3,  II-4  and  II-5. 

In  the  range  34.67  <  t  <  138.01,  the  values  of  the  dis¬ 
placement  of  the  free  surface  (r,  t)  have  also  been  computed 
from  Eq.  (II- 1)  using  the  Bendix  G-15  computer.  However,  be¬ 
cause  of  their  bulk  and  the  fact  that  they  do  not  have  the  general 

interest  of  Table  11-4,  these  results  have  not  been  presented  in 

e 

this  report.  Also,  only  the  wave  envelopes  B  (r,  t)  and  wave 
period  distribution  are  used  for  practical  purposes.  The  values 
of  7)  (r,  t)  have  been  plotted  only  for  Case  I  in  Fig.  II- 6  as  a 
sample  of  the  obtained  results.  Similar  graphs  for  Cases  II  and 
III  may  eventually  be  obtained  from  the  performed  computations. 
Fig.  II- 6  is  very  difficult  to  plot  because  of  the  large  number  of 
crests  and  troughs  (56  crests  and  55  troughs).  A  more  accurate 
plot  of  Fig.  II- 6  could  be  obtained  by  increasing  the  number  of  en¬ 
tries  in  the  performed  computations.  This  would  require  about 
four  times  as  many  computations,  i.e.  o*  would  have  to  be 
incremented  by  0.0025  in  the  intervzd  from  0  to  4. 

The  values  of  7f  (r,  t)  for  t  >  138.01  have  not  been  com¬ 
puted.  However,  the  times  of  arrival  of  the  crests  and  troughs 
may  be  computed  using  Eq.  (11-19)  for  positive  integers  m  >113. 
While  the  times  of  arrival  of  crests  and  troughs  are  the  same  for 
all  three  cases,  it  should  be  noted  that  a  crest  in  Case  I  may  cor¬ 
respond  to  a  trough  in  Case  U,  and  vice-versa.  For  example,  the 
time  of  arrival  =  262.  07  minutes  corresponds  to  a  trough  in 

Case  I,  a  crest  in  Case  II,  and  a  crest  in  Case  III.  The  question 
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WAVE  ENVELOPE  AND  WAVE  PERIOD  VERSUS  TIME  FOR  CASE  I 
(r  a  300  NAUTICAL  MILES) 


FIGURE  n-4 

WAVE  ENVELOPE  AND  WAVE  PERIOD  VERSUS  TIME 
FOR  CASE  II  (r  =  300  NAUTICAL  MILES) 
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FIGURE  U-5 

WAVE  ENVELOPE  AND  WAVE  PERIOD  VERSUS  TIME  FOR 
CASE  III  (r  =  300  NAUTICAL  MILES) 
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CASE  I  (r  =  300  NAUTICAL  MILES) 


of  whefJier  a  given  time  of  arrival  correaponda  to  a 
a  trough  ia  aettled  by  meana  of  the  following  rule: 

For  poaitive  integera  n,  the  timea  of  arrival  t^ 

o 

1  trougha  for  n  even  and  t^  <  t^  <  t^  ; 

2°  create  for  n  odd  and  t  <  t  <  t,  ; 

o  e  e 

3  trougha  for  n  even  and  t ^  <  t^  <  t^ : 

o  e  e 

4  create  for  n  even  and  t,  <  t  <  t^  , 

1  —  n  —  2 


creat  or 


repreaent 


U*27 


REFERENCES 


Kaplan,  Kenneth,  Norman  R.  Wallace  and  Thomas  C.  Goodale, 
0962),  "Studies  and  Analyses  of  Explosion-Generated 
Surface  Water  Waves,  "  United  Research  Services,  Bur¬ 
lingame,  California  (Preliminary  Draft). 


Kranaer,  Herbert  C.  and  Joseph  B.  Keller  (1955),  "Water  Waves 
Produced  by  Explosions,  "  N.  Y.  U,  Inst,  of  Math.  Sci. 
AFSWP  -  713. 


UST  OF  SYMBOLS 
APPENDIX  II 


X 

T 

h 

1 

R 

Vo 

t 

<r 

■E(s) 

A 

r 

0* 

^(o-) 


Distance  from  the  center  of  disturbance 
Wave  length 
Wave  period 
Depth 

Free  surface  elevation  above  still  water  level 

Effective  radius  of  disturbance 

Depth  of  disturbance  at  r  =  0 

Time  variable 
2  IT  h 

X 

Hankel  transform  of  function  E(r) 

Vo 


R  ys" 


2  TT  (y  -  ^) 


( 


a 


^(<t)  1  2 

♦'(O’)  J 
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APPENDIX  III 


SURFACE  WAVES  GENERATED  BY  DISTURBANCE  ON  SEA 
BED  IN  CONSTANT- DEPTH  OPEN  SEA 


by 


J.  A.  Hendrickson 


1. 


INTRODUCTION 


Even  though  the  sea  bed  disturbance  may,  in  actuality,  be 
of  a  complex  nature,  it  is  reasonable  to  assurhe  that  the  surface 
waves  generated  by  the  disturbance,  at  distances  far  removed  from 
the  source,  will  approximate  those  created  by  an  equivalent  source. 
We  will  assume  equivalence  on  the  basis  of  input  energy  and  geo¬ 
metrical  similitude.  Since  the  actual  source  (nuclear  explosion 
at  sea  bed)  is  three-dimensional  in  nature,  we  will  assume  an  equi¬ 
valent  source  in  the  form  of  a  cylindrical  upthrust  of  the  sea  bed. 
The  time  dependency  of  this  upthrust  will  be  assumed  to  be  ex¬ 
ponential  in  nature,  with  an  extremely  rapid  time  decay.  This 
appears  to  approximate  the  energy  dissipation  of  actual  nuclear 
blasts.  Finally,  in  order  to  complete  the  equivalence  between  our 
assumed  model  and  the  actual  disturbance  source,  one  can  match 
the  appropriate  fraction  of  the  nuclear  blast  energy  to  the  energy 
input  of  our  assumed  source. 

2.  DESCRIPTION  AND  SOLUTION  OF  PROBLEM 

Fig,  III- 1  shows  the  coordinate  system  and  nomenclature 
appropriate  to  the  assumed  disturbance. 

We  assume  the  source  is  described  as  follows: 


FIGURE  m-1 
NOTATION 


I 


m-2 


where 


-  8t 


g (t )  ■  I  -  • 
B 
0 


f(r) 


1  r  <  ro 
J  r  >  To 


and 


B 

fir)  =  ~2  ’  ’’  * 


The  field  equation  of  the  potential  function  ^  is  assumed 
to  be  of  the  time  independent  form 


4> 


I  b  <i> 


- 2  + - 

b  r  r  b  r 


-I-  - 


b*  4> 


=  0 


(III- 2) 


b  z 


while  the  boundary  conditions  appropriate  to  our  problem  may  be 
written  as  follows: 


and 


77  ‘  4>  Ot  Z  «  0 

't  ^Z 


g  A  4-  A  »  0 


(III- 3) 


ot  z  =  d 

z  ' 

Since  we  are  dealing  with  a  semi-infinite  media,  we  must  include 
all  possible  frequencies  and  wave  lengths  in  the  spectra  of 
generated  surface  waves.  Hence,  it  may  be  shown  that  tlie  appro¬ 
priate  bounded  solution  to  Eq.  (Ill- 2)  may  be  written  as 

(ni-4) 


r*  r*  r 

=  Jo  (ur)  A  (u, 


O')  sinhuz  +  B(u,(r)  cosh 


u  z 


cos  O'  t  du  do- 


III- 3 


It  is  easily  shown  that  the  boundary  condition  at  z  =  d  is 


satisfied  provided 


<T  sinh  ud  -  ug  cosh  ud 

B(u,a)  =  A(u,o-)  - 5 - 

ug  sinh  ud  -  cr  cosh  ud 


(III-5) 


From  Eq.  (III-l) 


*  . 


Bse 


-St 


r  <  r, 
r  >  r, 


It  may  be  easily  demonstrated  from  the  Fourier- Bessel 

theory  that  »  ^  may  equivalently  be  expressed 
'  t 


as 


*  2s"Br 


/-eo 


COS  (rt  der 
«r  T  s 


J,  (uO  Jjur)  du  (111-6) 


However,  from  Eq.  (III-4),  it  is  easily  seen  that 

no 


z  =  0 


COS  <rt  do- 


u  A  (u,(7- )  Jo  (ur)  du 


(III- 7) 


Hence,  using  the  first  of  Eqs.  (Ill- 3)  and  Eqs.  (Ill- 6)  and  (III- 7), 
it  is  seen  that  the  coefficient  A{u,  a  )  given  by  the  relation 


A  (u,o-)  = 


2s  B  To  J|  (u  To) 


w  u 


(<r^+  s^) 


(III- 8) 


III- 4 


Thus,  using  the  last  of  Eqs.  (Ill- 3)  and  Eqs.  (III-4),  (III-B) 
and  (III-8),  we  obtain,  after  some  simplification  and  rearranging 


2s*  B 


J|  (urj  J,  (ur) 
cosh  ud 


(111-9) 


du 


<r  cos  <rt  der 


2  2  1 

(cr  4-s  )  ug  tonh  ud 


If  we  let  X  =  ud,  and  integrate  Eq.  (III-9)  with  respect  to 
time,  we  obtain 


(lU-lO) 


2s‘B  r. 


i7« 


IT  d 


Jilx  ^  )  Jo(x^)dx 


cosh  X 


<T  sin  (rt  do" 


((r*+  s*) 


It  may  easily  be  shown  that  the  second  integral  in  Eq.  (Ill- 10) 
has  the  solution 

(in-ii) 


^co 

<r 

sin  <rt  dcr 

v 

_  ■m  _  ^ 

cos 

^  X  tonh  X 

\n 

—  exp  (-St) 

(<r  +  s^) 

O 

<r^  -  X  tonh  x 

2 

w 

s*  +  X  tonh  X  J 

We  are  interested  in  the  solution  to  Eq.  (Ill- 10)  for  the  case 
t  ^  0,  S  >  ^  I  Since  the  integral  of  Eq.  (Ill- 10)  is  extremely 
convergent  with  x,  we  may  substitute  Eq.  (Ill- 11)  into  Eq.  (Ill- 10) 
and  use  the  fact  that  S  >  >  I  with  the  approximate  result 
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» 


B(^) 


r  f  riit^  *ni/2 

Jilx-J)  CO*  ^xtonhx  dx 


cosh  X 


(III- 12) 


If  we  write 


W  *  IT  To  B  g  d 


(III- 13) 


where  -  density  of  fluid,  we  may  rewrite  Eq.  (111-12)  as 

follows: 


2K  (^) 


‘*1  ®  yx  tonh  X  )  dx 

cosh  X 


(III- 14) 


where  W  and  ^  *  - _ . 

K  •  - J  e  •  -|r  JTi' 

Z  tr  g  d 

Since  we  are  ultimately  interested  in  the  asymptotic  expan¬ 
sion  of  rj  for  >  >  1.  we  note  from  Eq,  (111-14)  that  the 
integrand  vanishes  as  x  — *  0.  Hence  the  only  contributions  occur 
when  X  >  0.  Thus  we  assiune  that  x  is  large  enough  for  the 
following  asymptotic  expansion  to  be  valid: 


(•Jx) 


y 


-TT 

WT  % 


COS 


_  JL 

4 


(111-15) 


Using  the  relation 


COS  0  COS  b  = 


2 


j^cos(o+b)  +  cos  (o  — b) 
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we  may  write 


j^cos  ^  j^cos  tonh  x  J 

»  ^  "l^cos (x +d>/TFanF~x  -  ^)  (III- 16) 

+  cos  {x-9y/x  tonh  X  — 

Using  Eqe.  (Ill- 15)  and  (III- 16)  with  Eq.  (Ill- 17),  we  obtain 
the  following  approximate  integral  form  for  the  surface  profile 
valid  for  >  I: 

(m-17) 
r* 


J,(x^){cos^(x+d>/x  tonhx  ~  cos-j(x~0>/x  tonh  x  -  ^)}dx 


-/%  cosh  X 


Since  1,  it  is  clear  that  the  contribution  of  the  term 

COS  -^(x  +  ^V^X  tonh  X  -  )  in  Eq.  (Ill- 17)  may  be  ignored  since 

a  slight  change  in  x  causes  the  integrand  to  rapidly  oscillate  from 
positive  to  negative  values  in  such  a  fashion  that  the  net  integral 
is  approximately  zero.  The  remaining  term,  however,  possesses 
stationary  points  (dependent  on  6)  and  will  hence  contribute  a 
finite  value  to  the  total  integral. 

Hence  we  may  write  Eq.  (Ill- 17)  as  follows: 


<//(x)  exp  (i-j)  ^  (x  )  dx 


(III- 18) 


in- 7 


where 


\p  {%) 


J|  {-^  X )  sech  X 


and 


</>(x)  =x-flyx  tanh  x  “ 


Eq.  (111-18)  is  of  such  a  form  that  the  method  of  "stationary 
phase"  may  be  used  to  obtain  a  satisfactory  solution.  Although  the 
method  will  not  be  described  herein,  the  results  will  be  used  in 
solving  Eq.  (111-18)  and  the  details  of  the  method  may  be  found  in 
Stoker  ( 1957). 

Given  an  integral  of  the  form  I(K)* 

-  0 

where  K  >  >  1  while  ^  (C.k)  is  not  a  rapidly  oscillating 
function  of  the  variable  ^  Then  the  value  of  I  (k)  may  be 
written  as  follows: 


I(k) 


=  Z  \/'(a^,k)  (- 


2  TT 


1/2 


4>  (a^) 


exp  < 


i(k<^{a^) 


4-0(^) 


(111-19) 


where  the  sum  is  taken  over  the  zeros  of  (^'  (^)  in  the 

region  a.^  ^  ^  b  at  which  ^"  (  a^.)  /  0  The  sign  of  the 

*  Stoker,  J.  J,  (1957),  "Water  Waves,  "  Pure  and  Applied  Mathe¬ 
matics,  Vol.  IV,  p,  181. 
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quantity  ir  /4  in  Eq.  (Ill- 19)  ia  taken  to  be  the  same  as  the  sign 
of  iOf)  •  Hence,  using  Eq.  (HI- 18)  we  desii'e  the  solution  to 
the  relation 


<i>  (x)  «  I 


6  lionh 
■  2  ^ 


X  +  X  sech 


X  tanh  X 


J. 


(Ill- 20) 


Let  be  the  r  roots  of  Eq.  (Ill- 20).  Hence 


e 


2  yg,  tonn  g^' 
tonh  4-  sech* 


(III-21) 


It  may  be  shown  that 


(m-22) 


/ (o  )  -  - 
r  r  2 


■  2  12  2 

tonh  0-4- a.  sech  a-  2sech  a- 

^  . .  m.Lm 


2  (a^  tonh  a,.) 


(a^  tonh  Of ) 


(l-a^  tonh  a,) 


and  it  is  to  be  noted  that  )  ^  0  >  and  that  there  is  but 

one  solution  to  Eq.  (III-21)  for  a  given  value  of  6^1. 

Note  that  the  value  9=1  corresponds  to  the  arrival  of  a  surface 
disturbance  propagating  with  a  phase  velocity  of  value  gd  . 

Finally,  using  Eq.  (IH-M)  we  find,  for  >>  1,  the 
asymptotic  behavior  of  to  be  the  following: 


III- 9 


(Ill- 23) 


V 


K  ( 


(Oi 


JtL 

A. 


g(a,)  cos( -y  f  (a^)) -t-Ol-j) 


3/2 


where 


g(a,)  = 


2  sech  a 


and 


f  (a,)  »  9  tonh 


It  may  easily  be  shown  that  Eq.  (Ill- 23)  has  its  maximum  at 
0^1  and  that  the  peak  amplitude  (M.  W.  L.  to  crest  elevation)  is 
given  by  the  expression 


V 


9-1 


(III- 24) 


The  interesting  thing  to  note  from  Eq.  (111-24),  remembering 

w 

that  K  =  —r  X  cst,  is  that  the  peak  surface  disturbance  is  only 
d^ 

dependent  on  the  source  energy  (W),  the  water  depth  (d),  and 
the  ratio  of  radial  position  to  water  depth  (~^)>  Also  the  amplitude 
is  inversely  proportional  to  the  radial  position  (in  contrast  to  the 
inverse  cube  root  of  position  for  a  two-dimensional  source). 


Ill- 10 


Because  of  the  rapid  convergence  of  the  integral  (see  Eq. 

(Ill- 15)  )  it  is  to  be  expected  that  the  surface  disturbance  is  rela- 

r 

lively  insensitive  to  the  ratio  (radius  of  surface  disturbance 

to  water  depth)  for  values  of  <  1.  (This  is  so  due  to  the 

To  J,  {  ^  K  ) 

relative  independence  on  of  the  function  — L_lJl -  ), 

IT'S-  TIT  -»  _i - *1 —  — 1-*.: —  -  yg  0  gg  obtained  from 


Fig.  111-2  shows  the  relation 


I* 

Eq.  (111-23)  for  =  10  and  -jp  <  1.  Also  shown  in  Fig^  III- 2 

is  the  value  of  the  integral  of  Eq.  (111-14)  vs.  6  for  -j-  =  10  and 
*■0  1 

obtained  by  numerical  integration. 

It  may  be  seen  that  the  agreement  between  the  asymptotic 
approximate  solution  and  the  numerical  solution  is  not  too  good  for 
the  first  positive  peak.  However,  the  agreement  on  the  first 
negative  peak  is  quite  good  (10%  error).  Also,  the  agreement  on 
the  absolute  maximum  value  of  the  integral  is  quite  good. 

It  is  felt  that  better  agreement  would  be  obtained  for  larger 
values  of  since  such  a  case  would  improve  the  asymptotic 

approximations . 

Finally,  Figs.  Ill- 3  and  III- 4  show  the  results  of  the  asymptotic 
*^01  r 

integration  for  -g-  =  and  =  100  and  300  respectively. 


UI-ll 
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r 

z 


t 

8 


r 

o 

B 

4> 

d 

u 

cr 

g 

w 

K 

e 


Horizontal  distance  from  the  origin 
Vertical  coordinates  (z  =  0  at  the  bottom) 

Height  of  the  bottom  disturbance  as  function  of  time 
Time 

Coefficient  of  decay  of  the  bottom  disturbance  with 
respect  to  time 

Radius  of  the  original  disturbance 

Final  height  of  the  bottom  disturbance  when  time  tends 
to  infinity 

Velocity  potential  function 
Depth  (mid- water  level:  z  =  d) 

Space  frequency.  Dimensionally  equal  to  -y—  where 
L  is  length 

Time  frequency.  Dimensionally  equal  to  where 

T  is  an  interval  of  time 


Density  of  fluid 
Gravity  acceleration 


Energy  of  the  original  disturbance 


Coefficient  = 


r  ^B 
o 


%  gd“ 


Dimensionless  time  — 

r 
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A  PPENDIX  IV 


THE  PRINCIPLE  OF  SUPERPOSITION  AND 
THE  THEORY  OF  CAUCHY- POISSON 


By 

B.  Le  Mehaute 


1. 


FOREWORD 


The  author  hae  received  eome  comments  on  the  validity  of  the 
following  calculations.  Some  consider  this  new  approach  as  exact  and 
powerful  as  the  theories  based  on  the  principle  of  the  stationary  phase. 

At  the  opposite,  some  others  consider  this  approach  as  invalid.  What¬ 
ever  its  exactness,  this  work  has  been  included  in  this  report.  The 
author  already  indicates  some  limits  of  validity  of  his  calculations.  It 
will  be  the  subject  of  further  investigation  and  up  to  other  investigators 
to  give  the  pros  and  cons  of  this  new  theory.  It  is  pointed  out  that  if 
this  theory  is  found  to  be  valid  in  some  respects,  it  may  become  a  power¬ 
ful  way  of  solving  many  related  problems.  It  has  been  seen  in  Appendix  I 
that  the  wave  motion  due  to  a  disturbance  of  infinitesimal  dimensions  has 
been  investigated  by  Cauchy-Polsson.  In  this  appendix  the  principle  of 
superposition  is  applied  in  order  to  calculate  the  wave  motion  due  to  a 
disturbance  of  finite  dimension.  It  is  demonstrated  that: 

a.  The  elevation  becomes  instantaneously  infinite  at  the 
edge  of  the  initial  disturbance,  whatever  the  intensity,  the  area  and  the 
shape  of  this  disturbance.  So,  the  splash  phenomenon  is  demonstrated 
and  explained. 

b.  Far  from  the  impulse  the  wave  pattern  does  not  de¬ 
pend  upon  the  area  and  shape  of  the  impulse,  but  as  a  first  approxima¬ 
tion  its  amplitude  is  proportional  to  the  total  value  of  its  intensity. 

The  wave  motion  caused  by  some  particular  simple  shape  of  an 
initial  disturbance  of  finite  dimensions  is  calculated. 


IV.  1 


Then,  a  diacueslon  give*  the  limits  of  validity  and  the  physical 
meaning  of  the  presented  mathematical  construction  based  on  the  prin¬ 
ciple  of  superposition  and  the  theory  of  Cauchy-Poisson. 


2.  INTRODUCTION 


The  theory  of  waves  produced  by  an  impulse  or  an  initial  local 
disturbance  on  a  free  surface  has  been  investigated,  particularly  by 

i 

Cauchy-Poisson  in  the  case  of  a  three-dimensional  motion,  and  by  Lamb 
in  the  case  of  a  two-dimensional  motion.  (In  the  following,  the  notations 
of  Lamb  are  used).  In  both  cases,  the  solutions  have  been  obtained  by 
neglecting  the  convective  inertia  and  friction  forces,  thus  the  theory  is 
linear  and  valid  for  slow  motion. 

The  calculations  can  be  carried  out  by  application  of  the  Fourier 
double  integral: 


f(d)  cos  k(X-a)  da 


for  a  two-dimensional  motion,  or 

tw)  H  dk  /  f(a)Jo(ko)  ada 


Jm  = j  Jo(kw)  k  6k  J 


for  a  three-dimensional  motion. 

It  has  been  assumed  that  the  initial  disturbance  or  initial  impulse 
is  confined  to  the  immediate  neighborhood  of  the  origin  so  that  f(  a  ) 
vanishes  for  all  but  infinitesimal  values  for  a  (See  Fig.  IV- 1) 

That  is: 


f(a)da=  I 


*  Lamb,  H.  (1932  Edn.)  Hydrodynamics ,  Cambridge  Univ.  Press,  Eng, 


1V.2 


in  the  case  of  a  two-dimenaional  motion,  and 

Ha)2ra  da  f(d)2wada  ■  I 

in  the  case  of  a  three-dimensional  motion  where  a 
infinitesimal. 


is  always  an 


FIGURE  IV- 1 

DEFINITION  OF  AN  ELEMENTARY 
FREE  SURFACE  DISTURBANCE 


From  this  asstimption,  it  has  to  be  noted  that  since  the  integral 
has  a  finite  value  equal  to  unity  and  since  a  j  is  an  infinitesimal, 
f  (  a  )  must  be  infinity.  This  assumption  is  in  complete  contradiction 
with  the  assumption  of  linear  theory  which  requires  the  motion  to  be  very 
slow  and  even  infinitesimal  in  order  to  be  exact. 

Later  on,  other  studies  were  made  to  investigate  the  wave  motion 
due  to  a  disturbance  of  finite  dimensions.  These  were  also  based  on  the 
assvunption  of  linear  theory  and  on  the  use  of  the  Fourier  integral.  In 
this  case  the  initial  disturbance  was  not  assumed  to  be  confined  to  the 
immediate  neighborhood  of  the  origin,  but  instead,  it  was  assumed  to 
be  distributed  on  a  given  range. 


IV- 3 


Hence,  the  later  theories  were  more  realistic  and  more  in 
accordance  with  physical  phenomena.  Unfortunately,  the  difficulty 
in  the  evaluation  of  the  Fourier  integral  limited  the  initial  disturbance 
f(  a  )  to  relatively  simple  form. 

In  this  appendix  an  attempt  has  been  made  to  find  a  more  gen¬ 
eral  method  for  application  to  more  complicated  forms  of  initial  dis¬ 
turbance  including  those  varying  with  time.  It  is  the  results  of  these 
investigations  which  are  presented  here. 

This  method  is  based  on  the  principle  of  superposition  and  the 
theory  of  Cauchy- Poisson.  Indeed,  since  this  theory  is  linear,  the 
principle  of  superposition  is  theoretically  valid  and  the  wave  motion 
caused  at  a  given  time  and  at  a  given  location  by  a  definite  disturbance 
varying  with  time  should  be  the  sum  of  wave  motions  due  to  infinitesi¬ 
mally  small  disturbances  occurring  at  various  initial  times  t^. 

Here  a  quotation  of  Lamb  is  worth  recalling: 

"In  any  practical  case,  however,  the  initial  elevation 
is  distributed  over  a  band  of  finite  breadth;  we  will  denote 
this  breadth  by  1.  The  disturbance  at  any  point  P  is  made 
up  of  parts  due  to  the  various  elements,  da  ,  say,  of  the 
breadth  1;  these  are  to  be  calculated  by  the  preceding  for¬ 
mulae,  and  integrated  over  the  breadth  of  the  band.  In  the 
result, the  mathematical  infinity  and  other  perplexing  pe¬ 
culiarities,  which  we  meet  with  in  the  case  of  a  concentra¬ 
ted  line- source,  disappear.  It  would  be  easy  to  write  down 
the  requisite  formulae,  but,  as  they  are  not  very  tractable, 
and  contain  nothing  not  implied  in  the  preceding  statement, 
they  may  be  passed  over.  " 
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GENERAL  PROCESS  OF  CALCULATION 


It  it  recalled  that  the  aolutions  for  the  elevation  ae  a  function 
of  time  t  and  distance  from  the  origin(  X  or  0  )  are,  in  the  case 
of  a  two-dimensional  motion 


(iv-i) 


and  in  the  case  of  a  three-dimensional  motion 


. ] 


(IV-2) 


From  the  assumption  that  the  motion  is  linear,  it  is  evident 
that  if  the  integral  of  initial  elevation  (or  impulse)  has  a  value  A  in¬ 
stead  of  unity: 


f  (a)  da  ■  A 


f(a) 2irada  ■  A 


tf  and  (  are  given  by  the  above  formulas  (IV- 1)  and  (IV-2)  re¬ 
spectively  multiplied  by  A. 

These  two  e]q>ressions  may  be  written  as 


A 


l.3.5.7...(4n*l)  I2x) 


(IV- 3) 


IV. 5 


and 


Now,  consider  the  cases  where  the  initial  disturbance  is 
located  at  a  point  x  =  X,  in  the  case  of  a  two-dimensional  motion 
(See  Fig.  IV-2),  and  at  a  point  Z  =  in  the  case  of  a  three- 

dimensional  motion  (See  Fig.  IV-3).  Z  is  a  complex  number. 


FIGURE  IV-2 

NOTATION  FOR  TWO-DIMENSIONAL  DISTURBANCE 


FIGURE  IV-3 

NOTATION  FOR  THREE-DIMENSIONAL  DISTURBANCE 
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The  above  formulae  (IV- 1)  and  (IV-2)  are  valid  provided  x 
and  V  are  replaced  by  X-X  •  v  -  pe*^  respectively.  The 
argument  for  v  is  taken  to  be  sero  by  reason  of  symmetry. 

A  similar  change  of  origin  of  time  may  also  be  inserted  by  re¬ 
placing  t  by(t-tj)  involving  an  initial  perturbance  at  time  t^  In¬ 
serting  these  values  and  assuming  that  A  depends  upon  the  location 
X,  or  (  P$  6  )  and  time  t^  one  obtains: 


n*o 


From  the  assumption  that  the  theory  is  linear,  the  principle 
of  superposition  may  be  applied,  and  it  is  deduced  that  the  free  surface 
is  obtained  by  a  'iuperposition  of  all  the  perturbances  created  at  differ¬ 
ent  places  and  at  different  times.  In  the  following  t^  vdll  be  taken 
equal  to  zero. 

In  particular,  one  can  imagine  that  there  is  a  perturbance 
A(X)  at  any  place  -R  $  R 

at  any  place  <  p  ^  R 

lO 
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such  that 


A(X) 


f(X)  dX 


np,B)  de  6p 


In  that  case  i(  X  )  and  f(  p,B)  do  not  need  to  be  infinite  to 
produce  finite  wave  and  consequently  A  may  be  infinitesimally  small, 
since  it  is  now  the  product  of  a  finite  value  and  an  infinitesimally 
small  area.  However,  the  total  values 


f(X)  dX 


and 


2.ir 

Hpfi)  AB  dp 


in  which  R  has  a  finite  value,  have  a  finite  value  and  involve  a  finite 
amount  of  energy,  which  is  the  total  energy  of  the  initial  perturbance. 

Now  at  a  given  time  t  and  at  a  given  place  x  or  (  v,o  )  the 
wave  motion  will  be  according  to  the  principle  of  superposition,  the 
sum  of  all  the  elementary  wave  motions  caused  by  the  various  places 
X  or  (  p,B  ). 

Hence,  17  and  ^  are  given  by  the  following  integrals  vhere 
X  and  t  are  considered  as  constant  and  X  or  (  p,  9  )  as  variable: 


rR  ^ 

/  _f(X)_V 

/  ir(X-X) 

-'-R  n«c 


(-1) 


1.3  ..(4n  ♦!) 


gt^ 

2(x-X)  . 


2n»  I 


dX 


7  (-1)" 

.1.3  ...(2n*l). 

2 

1  2  (W.p9'¥l 

(4n  ♦  2)1 

_  w-pe'*^  J 

2n«l 


dB  dp 


n»o 
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or  rearranging: 


•  .  XV  (-0*  I  gtMg">«r^  f(X)dX 

^  »/_.  I.3....(4n*l)  I'T)  J  (x.X)®"*2 

n«o  -'-R 

j_^  (-l)"[l  3  .(2n  J)]^  [  qfg  f(pg)dgdp 

"  2»/_^  (4n*2)l  1^1  J  J 

n»o  -'o  -'o 


Without  any  more  calculation,  it  may  be  seen  that  when  x  or 
or  ie  large  with  respect  to  the  maximum  values  for  X  or  p  that  is 
R,  such  that  X  and  p  may  be  neglected  in  the  expressions  (x-  X) 
and  (  er.pe'^  )f  the  obtained  formulas  for  if  and  {  are  the 
same  as  formulas  (IV- 1)  and  (IV- 2)  above  provided  the  second  term 
is  multiplied  by 


f(X)dX 


and 


2ir 

HpjSidQ  dp 


respectively;  that  is,  by  the  total  value  of  the  initial  disturbance. 

This  means  that  far  from  the  initial  disturbance  the  wave  pattern  does 
not  depend  upon  the  shape  and  stretch  of  the  disturbance,  and  the  ele¬ 
vation  is  linearly  correlated  with  the  total  impulse  or  disturbance. 
This  result  is  consistent  with  the  above  remark  on  the  multiplication 
factor  A.  But  in  that  case,  this  multiplication  factor  is  the  sum  of 
all  the  infinitesimal  local  values  A  (X)  or  A  (  p,6  ). 


IV.  9 


Moreover,  it  may  be  noticed  that  when  x  and  «  approach 
R,  then  t;  and  (  tend  to  infinity  (except  when  t  «  0),  since  the 
numerator  has  a  finite  value,  whatever  the  shape  for  f(x)  and 
f(  p,9  ).  This  is  in  accordance  with  the  well  known  fact  that  when 

a  circular  plate  falls  into  water,  the  free  surface  suddenly  rises  very 
high  in  the  shape  of  a  dvike's  crown  at  the  edge  of  the  plate. 

Later,  the  validity  of  such  a  theory  will  be  discussed.  Now 
some  specific  simple  cases  are  calculated. 


4.  THE  ANALYSIS  OF  SOME  SPECIFIC  CASES 

First  the  case  of  a  two-dimensional  motion  caused  by  a  uniform 
disturbance  is  analyzed.  In  this  case  f(X)  is  a  constant:  say  f(X)  -  b 
(See  Fig.  IV. 4). 


Then: 


FIGURE  IV-4 

MATHEMATICAL  MODEL  DISTURBANCE 
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The  following  general  integral  permits  the  direct  establishment 
of  the  complete  series: 


11  *0 

or  developing 
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It  is  verified  that  when  X  -—  +  R  ■»?  — ♦  »  except  when  t  =  0. 

2  2 
On  the  other  hand,  when  x  tends  to  be  large  with  respect  to  R  , 

it  is  found  after  a  number  of  simplifications: 


which  is  formula  (IV-3)  in  which  2Rb  has  taken  the  place  of  A; 
only  the  elevation  depends  upon  the  value  of  the  disturbance,  but  the 
wave  pattern  is  Identical. 

Now  considering  the  base  of  an  impulse  in  the  form  of  parabola; 


1 


f(X)  -  b  -  oX® 
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Then 
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The  general  integral  value  In  ie  found  to  be 
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finity.  Similarly,  after  a  ntimber  of  computations  (neglecting  R 
before  x^),  it  is  found  that  when  X*  »R* 


2(bR-^) 


Inserting  this  value  ,  the  Eq.  (IV-3)  for  rj  is  recognized  pro¬ 

vided  A  is  replaced  by  2(bR-  ¥> 

It  may  be  verified  that  this  is  the  total  value  of  the  initial 
disturbance  given  by  the  integral: 


r”  r" 

/  f(X)dX  •  /  (b 
•J-R  -^-R 


aX*)dR  *  2(bR  -  -—) 
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The  case  of  an  initial  diaturbance  in  the  form  of  a  cylinder 
may  also  be  easily  analyzed  by  calculating  the  integral 


I 


rZir 

.  /  /  dg  d/> 

Jo  k 


The  result  is  obtained  in  the  form  of  a  recurrence  formula,  the  number 
of  terms  increasing  with  the  value  for  . 

The  length  of  the  obtained  formulas  does  not  justify  writing  them 
down  in  this  appendix  since  they  contain  nothing  not  already  implied  in 
the  previous  formulas. 


5.  DISCUSSION  OF  THE  VALIDITY  OF  THIS  PROCESS  OF 

CALCULATION 

Now  it  is  time  to  question  the  validity  of  such  a  theory  as  It  has 
been  developed  for  the  Cauchy-Poisson  theory.  Does  this  mathematical 
model  have  a  physical  significance? 

First  of  all,  it  is  recalled  that  the  Cauchy-Poisson  wave  motion 
involves  an  infinite  amount  of  energy.  As  previously  mentioned,  this 
paradox  is  explained  by  saying  that  the  initial  disturbance,  being  in¬ 
finitesimally  small  in  width,  must  have  an  infinite  amplitude  for  its 
total  area  to  be  finite.  It  has  been  noted  previously  that  this  statement 
is  in  complete  disagreement  with  the  fact  that  the  theory  is  based  on 
the  assumption  that  the  motion  is  infinitesimal  everywhere,  including 
the  origin.  Moreover,  as  also  previously  mentioned,  the  fact  that  the 
initial  disturbance  is  exerted  instantaneously  involves  an  infinite  power 
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during  an  infinitesimal  time  interval.  This  power  causes  an  infinite 
velocity  at  origin  which  is  again  incompatible  with  the  assumption  of 
linear  theory. 

The  case  of  an  initial  disturbance  of  finite  dimension,  as  it  has 
been  analyzed  in  this  paper,  also  contains  a  paradox:  as  it  has  been 
e3q>lained,  the  value  for  f(X)  does  not  need  to  be  infinite  to  create  a 
wave  pattern  of  finite  dimensions.  Even  this  statement  permits  the 
principle  of  superposition  to  be  theoretically  valid  as  a  first  approxima¬ 
tion.  It  would  not  be  valid  if  f(X)  was  infinite.  However,  it  is  found 
that  the  amount  of  wave  energy  produced  by  an  initial  disturbance  of 
finite  energy  is  infinite,  and  in  particular  that  the  elevation  at  the  edge 
of  the  disturbance  tends  to  infinity. 

In  the  case  where  the  initial  disturbance  is  distributed  over  a 
finite  breadth,  we  again  quote  Lambs  "the  mathematical  infinity  (for 
f(a)  )  and  other  perplexing  peculiarities,  which  we  meet  with  in  the 
case  of  a  concentrated  line- source,  disappear."  Indeed  they  should 
disappear;  however,  they  do  not  disappear.  Hence  these  considerations 
have  led  the  author  to  maintain  that  the  Cauchy-Poisson  theory  is  a 
pure  mathematical  construction  which  can  represent  only  a  limited 
amo\mt  of  physical  facts.  It  must  be  recalled  here  that  frictional 
effects  are  neglected.  The  infinite  amplitude  of  the  initial  disturbance 
and  its  infinite  power  are  not  compatible  with  the  assumption  that  an 
infinitesimal  motion  is  linear,  thus,  problem  presents  a  singularity 
at  the  origin.  This  limits  the  validity  of  the  theory  and  does  not  per¬ 
mit  strict  application  of  the  principle  of  superposition. 
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However,  a  number  of  results  do  have  a  physical  significance. 
These  are:  (1)  far  from  the  origin  the  wave  motion  does  not  depend 
upon  the  width  and  shape  of  the  initial  disturbance.  Its  amplitude, 
under  the  assumption  that  the  wave  motion  is  small,  is  proportional 
to  the  total  value  of  the  impulse;  (2)  the  value  for  the  elevation  at  the 
edge  of  the  initial  impulse,  which  is  zero  at  time  zero,  tends  at  once 
to  infinity  at  time  t  =  c  .  despite  the  fact  that  the  initial  impulse 
has  a  limited  energy.  Indeed,  it  is  necessary  to  keep  in  mind  the  fact 
that  this  limited  energy  is  theoretically  exerted  instantaneously,  that 
it  involves  an  infinite  power.  But  since  the  amount  of  available  energy 
is  finite,  this  instantaneous  rise  of  water  must  last  a  very  short  time. 

In  fact,  this  is  observed  when  a  body  is  thrown  into  water.  Water 
rises  very  high  at  the  edge  of  the  body,  then  falls  down  and  is  followed 
by  small  oscillations  about  the  still>water  level.  It  is  evident  that  the 
friction  forces  also  have  a  definite  influence  in  causing  this  damping. 

6.  CONCLUSION 

To  conclude  this  discussion,  the  present  theory  based  on  the 
principle  of  superposition  permits  an  explanation  of  the  splash  of  water 
caused  by  a  disturbance  of  finiH  dimensions. 

Incompatibility  between  the  linearization  and  the  infinite  power 
at  the  origin,  and  friction  forces,  limit  its  validity  for  describing  the 
true  physical  situation.  Moreover,  it  has  always  been  assumed  that 
the  initial  disturbance  acts  instantaneously.  In  the  case  of  a  disturbance 
of  finite  dimension,  the  inertia  of  the  involved  mass  also  has  an  effect. 
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It  has  been  seen  in  Appendix  I  that  other  theories  exist  which  seem  to 
give  fairly  good  agreement  with  physical  facts  but  their  validity  is 
limited  to  the  large  values  for  x  and  w  respectively.  These  do  not 
e3q>lain  the  splash  phenomenon  as  has  been  attempted  in  this  paper. 

These  other  theories,  based  on  the  use  of  the  Fourier  integral 
in  which  f(  a  )  is  not  confined  to  the  neighborhood  of  the  origin, 
could  probably  be  generalized  as  indicated  in  this  paper  for  a  disturb' 
ance  of  any  shape  which  could  also  be  a  function  of  space  and  time. 
Their  investigation  may  also  be  of  particular  interest  for  calculating 
the  long  waves  due  to  a  traveling  atmospheric  perturbance  varying  with 
time. 
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LIST  OF  SYMBOLS 


APPENDIX  IV 

X  Horizontal  coordinate  in  a  two-dimensional  motion 

07  Horizontal  coordinates  in  a  three-dimensional  motion 

(complex  number) 

a  Horizontal  spread  of  an  original  disturbance 

1{o)  Function  characterizing  the  original  disturbance 
Tf  Free  surface  elevation  in  a  two-dimensional  motion 

C  Free  surface  elevation  in  a  three-dimensional  motion 

A  Value  for  the  initial  total  impulse 

t  Time 

n  Integer  index 

X  Coordinate  for  an  elementary  initial  elevation  into  a  two- 
dimensional  motion 

(  p,d  )  Cylindrical  coordinates  for  an  elementary  initial  elevation 

into  a  three-dimensional  motion  (  p  :  modulus,  9  :  argument) 

Z  (complex  number) 

R  Horizontal  spread  for  the  original  disturbance 

b  Height  of  the  original  disturbance 
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APPENDIX  V 


THE  SHOALING,  DAMPING,  BREAKING  AND  RUN-UP 
OF  LONG  WAVES  OVER  THE  CONTINENTAL  SHELF 


ON  SATURATED  AND  NONSATURATED  BREAKERS 


by 

B.  Le  Mehaut^ 


1. 


INTRODUCTION 


In  this  appendix  the  problems  associated  with  a  wave  traveling 
over  the  continental  shelf  are  analyzed. 

First  a  choice  must  be  made  from  among  all  existing  theories 
of  waves  traveling  over  gentle  slopes  in  order  to  select  the  most  con¬ 
venient  theory  for  application  to  this  problem.  It  is  shou'n  that  in  most 
cases  the  waves  appear  as  a  succession  of  solitary  waves  and  from  ex¬ 
isting  eiq>erimental  data  it  is  demonstrated  that  no  reliable  theory  exists. 

Wave  damping  due  to  bottom  friction  is  analyzed  and  a  formula 
is  proposed  which  takes  into  account  bottom  friction  and  shoaling  effects. 

The  very  important  distinction  between  saturated  and  nonsatu- 
rated  breakers  is  introduced.  This  study  is  of  major  importance  in 
this  report  because  the  demonstration  of  the  natural  protection  afforded 
by  the  continental  shelf  relies  on  this  theory. 

Finally  a  brief  literature  survey  of  previous  work  done  on  the 
wave  run-up  due  to  long  waves  is  given. 

2.  THE  CHOICE  OF  THEORIES 

It  is  most  important  to  be  aware  of  the  most  convenient  existing 
theories  for  treating  the  problem  of  wave  motion  on  the  continental  shelf. 
These  wave  motions  of  interest  are  characterized  by: 

a.  Their  period  T  ranges  from  50  to  200  seconds. 

b.  Their  wave  height  H  ranges  from  0  to  150  feet. 

c.  The  depth  d  is  smaller  than  600,  and  most  often 
smaller  than  100  feet. 
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A  first  indication  of  the  essential  wave  characteristics  is  obtained 


from  the  graph  presented  in  Fig.  1-4  of  Appendix  I.  It  is  seen  that  the 

d  ^  H 


most  significant  factors  are 


and 


where  X,  is  the  wave 
7)  X^d"^  i.e. 


length  in  deep  water,  or  alternatively,  simply 

T  g  d”  with  X  =  T  ^^d  where  is  the  elevation  of 

the  wave  crest  above  the  still  water  level.  Knowing  H  ,  can  be 

calculated  from  Fig.  1-3  of  Appendix  I.  Then  it  is  easily  seen  that  the 

value  to  be  considered  for  is  almost  the  same  as  the  value  for 

H.  When  - y  is  very  small,  the  smallest  value  which  occurs  in  this 

2  -2 

study  for  T  g  d  is  17H  over  the  continental  shelf.  Hence  it 

is  seen  that  the  waves  over  the  continental  shelf  can  be  considered  as  a 
succession  of  solitary  waves  and  the  solitary  wave  theory  should  be  a 
sound  basis  for  preliminary  analysis  of  the  phenomena  involved.  Later 
it  will  be  seen  that  this  statement  may  require  some  reservations  and 
that  many  refinements  will  be  necessary.  The  cnoidal  wave  theory  has 
been  shown  to  often  be  a  better  representation  for  the  wave  motion  in 
V6iy  shallow  water.  (See  Wiegel  (I960),  ) 


3.  ON  THE  SHOALING  EFFECT 
a.  Hydrodynamic  Solutions 

A  number  of  theoretical  studies  have  been  carried  out  in  recent 
years  on  the  problems  of  water  waves  on  a  slope.  In  particular,  signif¬ 
icant  contributions  have  been  made  byMiche  (1944),  Stoker  (1947), 
Biesel  (1952),  Peters  (1952),  Roseau  (1952),  Keller  (1958),  Carrier  and 
Greenspan  (1958)  and  Williams  (1959). 
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Solutions  have  been  obtained  by  direct  integration  of  the  momen¬ 
tum  equation  and  continuity  relationships  for  various  boundary  condi¬ 
tions,  including  that  of  an  inclined  sea  bottom.  Some  of  these  theories 
ars  solutions  of  linear  equations  valid  for  short  waves.  Hence  they 
cannot  be  applied  to  the  present  problem  because  of  the  relative  impor¬ 
tance  of  nonlinear  terms  (convective  inertia).  An  interesting  nonlinear 
solution  for  long  waves  has  been  developed  by  Carrier  and  Greenspan. 
But  the  vertical  acceleration,  which  is  important  for  near-breaking 
waves  (see  Appendix  VI),  and  the  damping  due  to  bottom  friction  have 
been  neglected.  Also,  the  solutions  would  not  hold  true  for  breaking 
inception. 

A  theoretical  solution  for  the  present  problem  has  also  been 
obtained  by  Kishi  (1962).  The  most  important  results  of  this  study 
are  reproduced  in  this  report.  Kishi  starts  from  the  basic  long  wave 


equations  in  the  usual  symbols: 

Momentum:  ^  ^  + 

b* 


=  -  g 


Continuity: 


^  [u  (d  -i  ^  )3  ^  Q 


with  c  s 


assuming  that 


u  « 


l]  ,  (V-l) 


it  is  found  that  for  two  locations  along  the  wave  path  (1  and  2)  and 
with  S  H: 
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which  permits  the  calculation  of  H2  as  a  function  of  from  the  know* 

resting  to  note  that  wfa 

■■  ^  ■  (4-) 


ledge  of  H,  and  d,  .  It  is  interesting  to  note  that  when  H/d  is  small. 

‘  ‘  H,  /  d,  \  1/4 

the  classical  Green's  law  is  found: 


Although  this  theory  may  be  quite  useful,  it  must  be  pointed  out 
that  the  derivation  of  Eq.  (V-Z)  is  based  on  the  assumption  (V-1).  Some 
doubts  can  be  e]q>ressed  as  to  the  validity  of  (V-1).  Also,  in  practice  the 
use  of  (V-Z)  would  involve  trial  and  error  computation  of  H2  in  terms  of 
Hj  ,  dj  and  d2  •  Because  of  these  two  objections,  Eq.  (V-Z)  is  not 
further  used  but  has  been  mentioned  as  a  possible  approach  for  further 
investigation. 

Also,  it  will  be  seen  in  the  following  appendix  that  the  vertical 
acceleration,  neglected  by  Kishi,  takes  on  great  importance  over  a  gentle 
slope  such  as  encountered  on  the  continental  shelf.  In  fact,  it  is  the  key 
for  explaining  the  paradox  of  long  waves  indicated  by  many  authors.  This 
paradox  --  the  Earnshaw  paradox  (1845)  --is  that  the  long  wave  theory 
gives  rise  to  a  bore,  whatever  the  wave  height  even  over  a  horizontal 
bottom.  In  fact,  it  is  known  that  some  tsunami  waves  over  a  steep  slope, 
as  on  the  Pacific  coast,  never  break.  The  sea  level  varies  gently  with 
time  only. 


b.  The  Energy  Method 

The  second  method  consists  of  assuming  that  the  wave  motion  on 
a  sloping  bottom  is  .the  same  as  on  a  horizontal  bottom.  Then  when  the 
wave  motion  has  been  so  determined,  it  is  assumed  that  the  rate  of 
transmitted  energy  is  constant. 
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The  great  advantage  of  this  method  is  its  simplicity,  permitting 
an  estimation  of  the  damping  due  to  friction. 

This  method  works  particularly  well  on  a  nearly -horizontal 
plane  as  is  encountered,  for  example,  on  the  continental  shelf  of  the 
Atlantic  coast  where  the  slope  is  between  2/10,  000  and  5/10,000.  Near 
(he  shoreline  the  slope  becomes  steeper,  ending  with  a  1/15  slope  on 
the  beach.  In  this  area  the  wave  motion  can  no  longer  be  considered 
as  similar  to  the  motion  on  a  horizontal  plane.  The  wave  profile  loses 
its  symmetrical  shape.  Also  a  part  of  the  wave  energy  is  reflected  in 
such  a  way  that  the  assumption  on  the  transmission  of  wave  energy  is 
not  fully  satisfied.  It  is  then  necessary  to  apply  an  analytical  or 
numerical  method.  Despite  these  limitations,  the  energy  method  is  a 
reasonable  preliminary  guide  for  studying  wave  motion  over  the  con¬ 
tinental  shelf. 

When  applied  to  the  solitary  wave  theory,  the  principle  of  con¬ 
servation  of  rate  of  transmitted  energy  gives 


E  V  b  =  — 5 —  d^^^  V  b  =  constant 

3^3 

where  H  is  the  wave  height 

V  is  the  wave  celerity  =  .^g  (d  +  H) 

b  is  the  distance  between  orthogonals,  proportional  to  the  radial 
distance  R  in  the  case  of  circular  waves  over  concentric 
bottom  contours. 

Then  between  two  points  (1)  and  (2)  at  radial  distances  and  R^ 


h; 


d,+ 


S')* 


d>H- 

2  2 


di  Ri  1  di  ^  Ra  ^ 

®2  *'2  ®2  "l 


(V-3) 


V-5 


c.  Experimental  Facts 

Experiment  shows  that  due  to  partial  reflection  over  a  slope, 
the  exponent  4/3  in  (V-3)  is  generally  too  high  (Ippen  and  Kulin  (1954)  ). 
It  can  only  be  considered  as  a  limit  when  the  slope  tends  to  zero.  The 
following  graphs,  based  on  experimental  results,  give  the  value  for 
n  as  a  function  of  the  slope  (Fig.  V-  1).  It  is  then  seen  that  over  the 
continental  shelf,  the  following  law,  where  n  =  1,  is  more  realistic; 


^  '  (4) 
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It  must  be  noted  that  such  a  law  will  also  be  in  accordance  with 

the  assumption  that  the  transmitted  energy  is  a  constant;  E  b  =  constant. 

This  assumption  is  far  more  reasonable  that  £  V  b  =  constant  because 

this  latter  law  implies  that  the  energy  of  the  wave  increases  when  the 

celerity  decreases,  which  is  not  physically  reasonable. 

Also,  according  to  some  experimental  results  presented  in  Fig. 

V-2,  i,  seems  that  the  law  of  variation  of  wave  height  with  the  depth 

follows  the  solitary  wave  theory  as  presented  in  Eq.  (V-3)  only  in  a 

narrow  range  near  the  breaking  depth  (Mur.k  (1949)  ).  According 

to  these  experiments,  when  -t—  is  larger  than  1.4,  the  Airy  law 

“b 

applies  despite  the  fact  that  d~^  is  much  larger  than  the 

required  value  presented  in  Appendix  I,  i.e.  0.01.  The  Airy  law  is 
based  on  the  assumption  that  the  rate  of  transmitted  energy  of  a  periodic 
gravity  wave,  given  by  the  linear  theory,  is  a  constant,  i.e.  V  b  = 
constant,  where  V  is  the  group  velocity . 

Applied  to  long  waves  (V  »  ^gd  ),  the  well-known  Green's  law 


is  obtained: 


FIGURE  V-1 

WAVE  HEIGHT  VARIATIONS  VERSUS  BOTTOM  SLOPE 


E  IN  DEPTH 


A  modification  of  Green's  law  is  suggested  to  take  account  of  the 


nonlinear  effect  as  follows; 


'  d2+  Hj' 


I  I 

^  b,  7 
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As  a  general  conclusion,  despite  the  fact  that  the  wave  motion 

over  the  continental  shelf  is  represented  fairly  well  by  a  solitary  wave, 

experiments  show  that  the  shoaling  effect  given  from  the  solitary  wave 

theory  cannot  be  applied  blindly.  The  shoaling  effect  is  certainly  very 

dependent  upon  the  slope  and  the  Airy  theory  seems  to  be  much  more  in 

accordance  with  fact  when  —t — ^  >  1.4. 

% 

It  must  be  noted  that  it  is  the  essential  characteristic  of  a  solitary 
wave  to  travel  without  deformation  on  a  horizontal  bottom.  Hence  these 
discrepancies  are  not  surprising.  The  real  law  can  only  be  given  by  a 
general  application  of  the  method  of  characteristics  with  very  small  inter¬ 
vals  with  the  use  of  a  digital  computer.  This  problem  will  be  considered 
in  another  appendix. 


4.  THE  DAMPING  EFFECT  DUE  TO  BOTTOM  FRICTION 

Most  of  the  theoretical  studies  on  wave  damping  have  been  carried 
out  for  periodic  waves:  Boussinesq  (1877),  Hough  (1877),  Basset  (1888), 
Lamb  (1932),  Keulegan  (1948),  Biesel  (1949),  Putnam  and  Johnson  (1949), 
Miche  (1954),  Reid  and  Bretschneider  (1954),  among  others. 

As  in  the  case  of  the  analysis  of  a  wave  on  a  sloped  bottom,  two 
methods  exist  to  attack  the  problem.  The  first  method  --  the  analytical 
method  --  consists  of  solving  directly  the  basic  differential  equations  -- 
momentum,  continuity  —  taking  a  friction  term  into  account.  Theoretically, 
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this  method  presents  the  advantage  of  giving  not  only  the  damping,  but 
also  the  deformation  of  wave  motion  due  to  friction  forces  (and  convective 
inertia).  This  is  important  for  very  long  waves  in  shallow  water  such 
as  tidal  waves  in  an  estuary.  This  method  is  essential  for  studying  the 
transformation  of  the  wave  into  a  tidal  bore. 

For  the  problem  under  study  the  second  method,  the  energy  method, 
is  simpler  and  will  be  as  accurate  for  the  wave  over  the  continental  shelf 
before  breaking.  It  consists  also  cf  the  determination  of  the  wave  motion 
independently  of  the  friction  forces  over  a  horizontal  bottom.  As  a  con¬ 
sequence,  the  wave  profile  is  symmetrical.  This  approximation  is  valid 
because  the  decay  of  wave  amplitude  over  a  wave  length  is  small  and  the 
slope  is  very  gentle. 

The  damping  effect  is  then  simply  defined  as  a  decay  in  wave  height 
calculated  from  energy  considerations:  the  loss  of  energy  over  a  given 
wave  length  is  equal  to  the  variation  of  wave  energy.  This  cotild  be  ex¬ 
pressed  in  two  ways:  the  first  is  particularly  convenient  for  periodic 


waves: 


(V-6) 


where  P  is  the  power  per  unit  width  or  energy  propagated  per  unit  time 
through  a  vertical  area  of  depth  d  and  unit  width: 

»  r  47rd  1 


P  =  ^  ^  "2  '  - - 4 


sin  h 


4  Trd 


where  H  is  the  wave  height,  V  the  wave  celerity  (  d  the  depth, 

L  the  wave  length  (T  x  is  the  distance  measured  along  the  wave 

ray  in  the  direction  of  propagation  of  the  wave,  b  is  the  distance  between 


wave  orthogonals  of  two  wave  rays.  Hence,  for  long  waves: 
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(V-7) 


I  H  ^  i 

p  *  ?  />  g  (y)  (gd) 


In  the  present  case  b  is  proportional  to  R,  i.e.  the  distance  from  the 
original  disturbance;  hence  b  can  be  replaced  by  R  in  Eq.  (V-5). 

is  the  average  amount  of  energy  dissipated  per  unit  area  at  the  bottom 
per  unit  time.  It  can  be  calculated  theoretically  when  the  flow  is  laminar, 
but  it  has  been  shown  that  the  flow  is  turbulent  when  (Collins  (1961)  ) 


>  1 . 08  sinh 


2  IT  d 


H  >  1.2  in  ft. 


sec.  units  for  long  waves.  This  condition  is 


always  verified  in  the  present  case.  Hence  will  instead  be: 


*  zL^  t/ 
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where  u„  is  the  bottom  velocity  and  T  the  shearing  stress.  Here 


for  a  coefficient  of  friction  f, 


=  P  1%  I 


and  according  to  the  Airy  theory  applied  to  long  waves: 


Ug  =  U  (x,t)  s  ^  \/^  (Itt-ni*) 


where  k  =  and  m  = 


^  "  T  yid 

Then,  introducing  these  values  it  is  found  that 

and  inserting  (V-7)  and  (V-11)  into  (V-6): 
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(V-10) 
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^(h'rVT).-^  Ih’  Rd' 
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Dividing  by  H^R’\ld,  integrating  over  a  small  interval  X2  -  s  Ax  =  Ar 


and  replacing  e 


>ax 


by  1  >  ax,  it  is  found: 


3ir 
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which  permits  step-by-step  calculation  of  H.  When  f  =  0,  the  classical 
Green  law  is  easily  recognized. 

This  formula  will  have  to  be  inverted  if  one  wants  to  evaluate  the 
deep  water  wave  height  as  a  function  of  the  shallow  water  wave  heights.  In 
particular  the  negative  sign  in  the  expression  between  brackets  becomes 
positive. 

Fur  the  solitary  wave  theory  it  is  more  exact  to  apply  the  same  type 
of  calculation  to  the  transmitted  energy  rather  than  to  the  power.  Then, 


d(EVb) 


-  b 


dE 
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dx  ''  dt 

The  energy  of  a  solitary  wave  is 

and  V  =  "^g  (d  +  H),  where  H  is  the  maximum  height  above  the  still  water 


3 

•2  d 
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level  of  the  solitary  wave,  and  d  the  depth.  (In  this  case  =  H) 

The  rate  of  loss  of  energy  due  to  bottom  friction  is: 


- 


dx 


(V-16) 


77 

By  inserting  the  classical  relationships  =  u  =  V  -g —  and  17  = 

/I  /  Hv 

(-3-) 


H 


where 


X  -  Vt 

— a — 


and 


cosh  a 

the  free  surface  elevation  above  the  still  water  level,  it  is  found  that 
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By  use  of  the  Cfude rmannian  of  a  it  is  found  that  the  Integral 


has  a  value  16/15.  Hence,  finally  . 

dL,  ^  ..  h!_NL 

dt  laTf  P  ^  1 

d  ^ 
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After  insertion  of  these  values  and  since  b  is  proportional  to 
R,  after  integration  it  is  found  that 


f  H,  AR 


22  “2  "2 
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This  formula  also  permits  step-by-step  calculation  of  H. 

It  is  seen  that  the  loss  of  energy  in  a  solitary  wave  height  H  is 
more  important  than  in  a  periodic  wave  of  the  same  wave  height. 

When  f  =  0  and  the  motion  is  two-dimensional,  the  classical  law 


d.+  H  ' 


(V-20) 


is  easily  recognized.  Now  the  problem  is  how  to  determine  the  friction 
factor  f  in  Eq.  (V-13)  and  (V-20). 

In  the  case  of  viscous  flow,  the  periodic  motion  in  the  boundary 
layer  on  a  smooth  plane  is  relatively  well  known.  The  thickness  of  the 
boundary  layer  8  is  proportional  to  -/v  T,  that  is,  it  increases  with 
the  period. 

In  the  case  of  turbulent  flow,  it  is  known  that  8  also  increases 
with  T.  The  exact  increase  of  8  1>  \mknown.  In  the  case  of  long 
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waves  over  a  shallow  bottom  it  is  reasonable  to  admit  that  all  velocity 
distribution  along  a  vertical  at  a  given  time  is  very  similar  to  that  of 
a  steady  flow.  In  a  word,  unsteady  motion  appears  as  a  succession  of 
steady  flows  insofar  as  the  friction  effect  is  concerned.  Then  f  can 
be  expressed  as  a  function  of  the  Chezy  coefficient  ,  which  can  be 

expressed  as  a  function  of  the  Manning  coefficient  n: 


2 

n 


(1.486* 


14.6 


(ft- sec.  units) 


(V-21) 


The  choice  of  n  is  delicate;  however  it  can  tentatively  be  con¬ 
sidered  as  equal  to  0.02,  which  is  usually  given  for  a  bottom  composed 
of  gravel.  But,  it  may  be  expected  to  range  between  0.015  and  0.025  due 
to  the  dual  uncertainty  of  the  bottom  roughness  and  the  application  of  the 
Manning  formula  to  unsteady  motion.  In  relatively  deep  water  the  depth 
d  in  formula  (V-21)  should  rather  be  replaced  by  the  thickness  of  the 
turbulent  boundary  layer  8  .  Unfortunately,  it  has  been  seen  that  8  also 
is  unknown. 

However,  it  must  be  mentioned  that  the  tidal  problem  in  an  estuary 
has  been  studied  satisfactorily  by  use  of  the  Chezy  (or  Manning)  coefficient. 
All  other  factors  being  equal,  only  a  slightly  different  value  for  has 
been  found  for  a  decelerated  flow  than  for  an  accelerated  flow. 


5.  ON  THE  BREAKING  INCEPTION  AND  NONSATURATED  BREAKERS 
It  is  commonly  admitted  that  breakers  on  a  beach  can  be  separated 
into  spilling  breakers  on  a  very  flat  slope  and  plunging  breakers  on  a 
steeper  slope.  Plunging  breakers  are  sometimes  called  surging  breakers 
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on  a  very  eteep  elope.  The  separation  of  breakers  into  categories 
is  based  on  visual  observations  rather  than  on  any  hydrodynamical 
criterion.  However,  the  essential  hydrodynamical  characteristics  of 
these  breakers  will  be  reviewed;  then  a  theory  for  nonsaturated 
breakers  developed  (Le  Mehaute  (1962)  ). 

The  profile  of  a  spilling  breaker  remains,  for  the  most  part, 
almost  symmetrical  and  the  wave  breaks  by  curling  over  slightly  at 
the  crest  (Fig.  V-3).  As  long  as  the  foam  of  the  breaker  is  small  by 
comparison  with  the  "bulk"  water,  which  happens  on  a  very  gentle 
slope,  the  wave  presents  roughly  the  main  characteristics  of  a  solitary 
wave,  even  after  breaking  inception.  But,  due  to  the  spilling  breaker, 
a  given  amount  of  energy  is  dissipated  in  such  a  way  that  the  wave 
crest  follows  the  breaking  index  curve  defined  by  H  =  0,  78  d.  Then 
the  spilling  breaker  is  transformed  into  a  bore  when  the  slope  be¬ 
comes  steeper.  When  the  slope  is  steep  before  breaking  inception, 
the  wave  profile  first  loses  its  symmetrical  shape,  then  a  plunging 
breaking  wave  generates  a  bore  directly. 

In  the  following  an  attempt  is  made  to  investigate  analytically 
these  described  phenomena.  As  before,  two  methods  exist.  The  first 
method,  the  energy  method,  is  only  approximate  but  gives  a  great 
amount  of  information  from  relatively  simple  calculations.  The 
second  method,  the  analytical  method,  is  more  accurate  but  requires 
tedious  computations  for  each  particular  case. 

It  will  be  seen  that  the  second  method  is  an  application  of  the 
method  of  characteristics  and  requires  some  refinements  for  analyzing 
the  surf  zone  and  calculating  the  wave  run-up.  This  will  be  presented 
in  Appendix  VI. 
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In  tlw  energy  method  the  rate  at  which  the  energy  Is  lost  due  to 

fie  dE 

bottom  friction  is  I  and  to  spilling  breakers  is  •  Hence: 


where  h^  and  are  the  depths  before  and  after  the  front  of  the  bore> 
respectively,  and  Q  is  the  discharge  due  to  the  moving  bore.  It  must  be 
remembered  that  the  above  formula  is  baaed  on  the  assumption  that  the 
vertical  distribution  of  the  horizontal  velocity  u  is  uniform. 

In  the  case  under  study  the  spilling  breaker  is  due  to  the  fact  that 
the  horizontal  velocity  at  the  crest  becomes  greater  than  the  wave  celerity 
V  .  By  analogy  (see  Fig.  V-3)  h^  =  d  +  H  and  h^  =  d  where  fi 

is  always  smaller  than  unity  and  can  be  zero  at  the  limit.  The  vertical 
velocity  distribution,  and  consequently  the  discharge,  is  directly  related 
to  the  average  horizontal  velocity.  Hence  the  discharge  could  be  written: 

0«  (d  +  H)  Uj  -(d  +  /9H)u,  •  V  [d  +  H  -  (d  + /3  H)] 

Inserting  these  values  into  Eq.  (V-23)  and  defining  B  as  follows: 
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(V-24) 

B  will  be  called  the  "breaking  coefficient."  The  breaking  coefficient  B 
is  the  ratio  of  the  rate  of  energy  dissipated  by  the  spilling  breaker  to  the 
rate  of  energy  which  could  be  dissipated  by  a  bore  of  front  height  equal 
to  the  height  of  the  solitary  wave  which  generated  it.  B  =  0  corresponds 
to  no  breaking  (  /9  =  1).  A  small  value  for  B  corresponds  to  a  little 
spilling  breaking  near  the  crest  (  close  to  unity):  i.e.  a  partial  breaking 
or  a  nonsaturated  breaker.  It  is  difficult  to  ascertain  the  maximum  value 
for  B  by  the  energy  method.  However,  it  is  certain  that  B  cannot  be 
larger  than  unity  {  0  -  0).  Then  total  breaking  occurs  and  the  breaker 

becomes  a  saturated  breaker  or  fully  developed  bore.  Further  consid¬ 
eration  will  be  given  to  the  physical  meaning  of  B  later. 

Now,  by  introducing  equalities  {V-15),  (V-18A)  and  (7-24)  into 
Eq.  (V-22),  it  is  found  that 


dx 


4  f 


H 


3 

d? 


3/2  3/2 

which  gives  after  division  by  H  d  V, 

A*£tX 

X  =  •  ^*^4  since  e  = 


integration  over  a  small  in- 
1  -  ax: 


5 

Ji  BH,"  Axi 

-TT - 1 - J 

d^ 
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When  all  friction  effects  are  neglected  (f  =  0)  and  there  is  no  breaking 
(B  =  0),  the  classical  law 
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is  easily  recognized.  It  has  already  been  pointed  out  that  such  a  law 
is  not  too  well  verified  experimentally  (see  Fig.  V-1  and  V-2). 

Despite  these  limitations  the  physical  interpretation  of  this 
study  will  be  based  on  Eqs.  (V-25)  and  (V-26)  because  the  spilling 
breaker  effect  tends  to  replace  the  variation  of  wave  height  by  a  simple 
law  expressed  as  H  -  0.78  d,  even  though  Eq.  (V-27)  is  not  verified 
experimentally.  When  H  <  0.  78  d,  there  is  no  breaking  and  the 

breaking  coefficient  B  =  0.  Then  V  =  ^  g  (d  +  H)  j  More¬ 

over,  assuming  H  is  small  by  comparison  with  d, 

d,  1  r  8  ^  1 

(V-28) 

°l-2 

This  equation  has  already  been  found  (c.f.  V-19). 

When  =  0.78  d^^,  there  is  inception  of  breaking  and  the  break¬ 
ing  coefficient  B  becomes  greater  than  zero.  In  the  case  of  a  small 
spilling  breaker,  V  retains  its  value  V  =  ^  g  (d  +  H)j  Then, 

replacing  these  values  for  H  and  V  in  Eq.  (V-25), 


d 

■3^ 


i.  e.  the  slope  S  = 


1.1  i  d’/"  4  0.07  Bd^'^ 
g 
d 


(d)  =  0.  01  f  +  0.  02  B  or  within  the  known  limits 


0  ^  B  =  SOS  -  O.Sf  < 


(V-29) 


It  is  seen  that  the  breaking  coefficient  B  increases  with  the  slope: 
the  steeper  the  slope,  the  greater  the  rate  at  which  the  energy  is  dis¬ 
sipated  by  the  spilling  breaker. 
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It  may  occur  that  due  to  bottom  friction  B  always  retains  a  aero 
value  despite  the  shoaling  when  S  <  0.01  f  as  is  easily  seen  from  Eq. 

(V-29)«  This  result  can  also  be  found  directly  from  Eq.  (V-25)  by  re-  “ 


placing  H  by0.78d  and  equating  B  to  zero. 

Inserting  the  value  (V-21)  for  f  ,  a  criterion  for  damping  without 


breaking  is  proposed: 


I4.6n 


100  d® 

i.e.  with  the  Manning  coefficient  n  =  0.02:  (d  in  feet) 

S<^ 

d3 


(V-30) 


On  the  other  hand,  it  has  lieen  seen  that  B  cannot  exceed  unity. 

This  happens  when  S  =  0.02-f0.01f  =  0.02.  When  S  =  0.02,  then 

the  breaker  is  "saturated."  Fig.  V-4  illustrates  these  considerations. 

Now  a  complete  physical  interpretation  can  be  drawn  from  the  pre- 

-5  1/3 

vious  considerations.  If  the  slope  is  always  smaller  than  6.10  /  d  , 

then  the  wave  height  is  completely  damped  by  bottom  friction.  There  is 
no  breaking  and  no  run-up.  This  occurrence  is  very  rare. 

On  a  steeper  slope  there  is  a  maximum  amount  of  wave  energy 
that  a  solitary  wave  can  transmit  towards  the  shoreline  over  a  given  depth. 
This  maximum  energy  is  reached  when  H  =0.78  d.  If  the  amount  of 
energy  passing  through  a  given  plane  tends  to  be  larger  than  this  maximum 
value,  a  spilling  breaker  will  dissipate  the  difference.  This  occurs  on 
a  relatively  gentle  slope  and  such  a  condition  represents  a  nonsaturated 


breaker,  in  which  case  the  wave  height  is  directly  related  only  to  the  depth. 
The  run-up  is  negligible. 
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FIGURE  V-4 

BREAKING  COEFFICIENT  VERSUS  SLOPE 


It  is  seen,  also,  that  there  is  a  limiting  amount  of  energy  which 
coulH  be  dissipated  by  a  breaker  over  a  given  length.  Hence,  when  the 
slope  becomes  steeper  and  steeper,  the  regulating  effect  of  the  spilling 
bread(er  reaches  its  limit  when  B  =  1.  Then  the' breaking  index  curve 
is  surpassed  by  the  height  of  the  bore  front.  There  is  run-up.  The  words 
"saturated"  and  "nonsaturated"  breakers  are  now  defined,  explained  and 
justified. 

A  very  important  conclusion  is  also  drawn:  on  a  beach  having 
its  curvature  upwards,  the  maximum  possible  wave  run-up  is  given  by 
the  wave  which  breaks  at  a  depth  where  the  slope  is  equal  to  0.02.  It  is 
known  that  if  d^  is  the  depth  over  that  slope,  the  corresponding  wave 
heigh,,  is  =  0.78  d^.  Any  wave  having  a  greater  height  breaks  sooner, 

dissipating  its  energy  following  the  breaking  index  curve  up  to  the  plane 
where  the  slope  becomes  larger  than  0.02. 

In  fact  the  theoretical  value  0.02  for  the  critical  slope  (corre¬ 
sponding  to  B  =  1)  may  be  replaced  by  a  more  factual  and  conservative 
value  0.01  .  The  exact  determination  of  this  value  requires  further  in¬ 
vestigation  by  the  method  of  characteristics  as  given  in  Appendix  VII. 

The  results  of  this  section  are  summarized  in  Fig.  V-5  by  three 
typical  cases.  It  must  be  noted  that  the  run-up  in  cases  I  and  II  is  the 
same  despite  the  different  deep  water  wave  heights. 

A  very  important  conclusion  can  be  drawn  from  this  study.  Since 
the  continental  shelf  has  a  very  gentle  slope  averaging  4/10,000,  it  acts 
as  a  natural  protection.  The  breaking  coefficient  B  keeps  a  very  small 
value  most  of  the  time.  The  breaker  is  far  from  being  saturated.  The 
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FIGURE  V-5 

ILLUSTRATION  OF  SATURATED  (BORE)  AND  NONSATURATED  (SPILLING)  B 


maximum  wave  destruction  depends  upon  the  depth  where  the  slope  is  near 
0.01  and  the  depth  at  this  slope  is  less  than  30  feet.  Hence  the  maximum 
possible  destruction  is  due  to  a  wave  height  based  on  H  =  0.78  d  or  0.78 
X  30  =  23  feet.  Also,  the  maximum  possible  wave  entering  the  Chesapeake 
Bay  is  also  directly  given  by  the  very  simple  rule  H  =  0.78  d.  This  wave 
is  a  gentle  spilling  breaker  in  the  shape  of  a  limit  solitary  wave. 

6.  A  LITERATURE  SURVEY  ON  THE  WAVE  RUN-UP  OF  LONG  WAVES 
No  theories  exist  for  calculating  the  wave  run-up  on  a  slope  prior 
to  this  study.  Appendix  VI  of  this  report  establishes  a  theory  for  this 
purpose.  Some  experimental  information  on  this  phenomenon  does  exist. 
Unfortunately,  these  experiments  were  carried  out  on  slopes  which  were 
too  steep  for  application  of  the  results  to  the  problem  under  study.  However, 
the  results  do  give  some  information  on  the  general  trend.  This  information 
will  be  applied  to  the  problem  under  consideration.  It  must  be  noted,  how¬ 
ever,  that  this  application  requires  some  extrapolation  and  is,  therefore, 
subject  to  possible  error. 

It  was  found  by  Kaplan  (1955)  that  on  a  slope  of  1/30 

0.381  (^)  (V-31) 

where  R  is  the  vertical  run-up  above  the  still  water  level,  and  H  and  L 
are  the  wave  height  and  wave  length  at  the  toe  of  the  slope.  (See  Fig.  V-6) 

On  a  slope  of  1/60,  he  found  that 

•g  =  0.206  (^) 


(V-32) 


dn-NflM  3AllV13d  ;H/a 
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The  interesting  fact  emerging  from  these  tests  is  that  the  value 

D 

of  does  not  depend  upon  the  depth  at  the  toe  of  the  slope.  These 

H  -2 

laws  were  verified  for  values  of  -j-  ranging  from  8  x  lO"  to 
Z  X  10*^.  They  were  also  proposed  for  extrapolation  for  smaller 
values  of  . 

Some  recent  experiments  on  a  slope  of  1/30  by  Kishi  (1962) 

•  3  H  -  3  R 

for  7  X  10*  <  -j-  <  2  x  10  gave  smaller  values  for  -pj-, 

R 

namely  =  2.  But  this  could  be  due  to  the  fact  that  Kishi  worked 

to  a  very  small  scale  for  which  scale  effects  may  not  be  negligible. 

Some  experiments  have  also  been  carried  out  by  J.  V.  Hall  and 
G.  M.  Watts  (1953)  on  the  run-up  of  a  solitary  wave.  They  give: 


R 

d 


0.67  u  |.9(S) 

IKS)  (^) 


0.35 
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for  any  slope  between  5°  and  12^.  However,  on  reworking  the  data 
provided  by  Hall  and  Watts,  it  is  found  that  -yj-  is  close  to  3  on  a 
slope  5°  and  10°. 

From  this  scattered  information,  combined  with  the  previous 
theory  on  nonsaturated  breakers  and  the  results  obtained  in  the 
following  appendix,  it  can  be  tentatively  concluded  that  for  an  average 
slope  as  usually  encountered. 


Hj^  being  the  wave  height  (0.  78  d)  at  a  depth  d  where  the  slope  is 

0,02. 


V-26 


REFERENCES 


Basset  (1888).  Hydrodynamique 

Bieseli  F.  (1949).  "Notule  Hydraulique. "  La  Houille  Blanche ,  Sep-Oct 

Biesel,  F.  (1952).  "Study  of  Wave  Propagation  in  Water  of  Gradxially 
Varying  Depth."  Gravity  Waves,  U.S.  National  Bureau  of 
Standards  Circular  !^21 

Boussinesq  (1877).  Essais  sur  la  theorie  de  eaux  courantes 

Bretschneider,  C.  L.  and  R.  O.  Reid  (1954).  "Changes  in  Wave  Height 
due  to  Bottom  Friction,  Percolation  and  Refraction."  Tech. 
Memo  No.  45,  Beach  Erosion  Board,  Dept,  of  the  Army  Corps 
of  Engineers,  36pp. 

Carrier,  G.  F.  and  H.  P.  Greenspan  (1958).  "Water  Waves  of  Finite 
Amplitude  on  a  Sloping  Beach."  Journal  of  Fluid  Mechanics, 

Vol.  4,  Part  1,  p,  97,  May 

Collins,  J.  I.  (1961).  "Effect  on  Mass  Transport  of  the  Onset  of 

Turbulence  at  the  Bed  Under  Periodic  Gravity  Waves."  ASME- 
EIC  Hydraulics  Conference,  Montreal,  Paper  number  61,  EIC-8 

Earnshaw,  S.  (1845),  Transactions,  Cambridge  Philosophical  Society, 
Vol.  VIII 

Hall,  J.  V.  and  G.  M.  Watts  (1953).  "Laboratory  Investigation  of  the 
Vertical  Rise  of  Solitary  Waves  on  Impermeable  Slope."  Tech 
Memo  No.  33,  Beach  Erosion  Board,  U.S.  Army  Corps  of 
Engineers . 

Hough  (1877).  Proceedings ,  London  Math.  Soc.  ,  Vol.  XXVIU,  p.  264 

Ippen,  A.  T.  and  G.  Kulin  (1954).  "The  Shoaling  and  Breaking  of  the 

Solitary  Wave."  Proceedings ,  Vth  Conference  on  Coastal  Engi¬ 
neering,  Grenoble,  France,  September,  p.  27 

Kaplan,  K.  (1955).  "Generalised  Laboratory  Study  of  Tsunami  Run-up. " 
Tech.  Memo  No.  60,  Beach  Erosion  Board,  Dept,  of  the  Army 
Corps  of  Engineers 

Keller,  J.  B.  (1958).  "Surface  Waves  on  Water  of  Nonuniform  Depth." 
Journal  of  Fluid  Mechanics,  V.  4,  Pt.  6,  Nov. 

Keulegan,  G.  H.  (1948).  "Gradual  Damping  of  Solitary  Waves."  Jour. 
Res. ,  Nat.  Bur.  of  Standards,  Vol.  40,  pp.  487-498 


V-27 


Kiahi.  Tautomu  (1962).  "Tranaformation.  Breaking  and  Run-up  of  Long 
Wavea  of  Finite  Height."  Preaented  at  Vlllth  Conference  on 
Coaatal  Engineering,  Mexico  City,  November,  1962. 

Lamb,  H.  (1932).  Hydrodynamlca .  Dover  Publlcatlona,  Sixth  Edition, 
p.  625 

Le  M^haute",  B.  (1962).  "On  Nonaaturated  Breakera  and  the  Wave  Run-up." 
Vlllth  Conference  on  Coaatal  Engineering,  Mexico  City,  Nov.  1962. 

Miche,  A.  (1944).  "Mouvementa  ondulatiorea  de  lamer  en  profondeur 
conatante  on  decroiaaante . "  Ann.  dea  ponta  et  chauaaeea, 
pp.  131-164,  270-292.  369-406 

Miche,  R.  (1954).  "Train  d'ondes  oceaniquea,"  COEC,  Miniate  re  defenae 
nationale.  No.  135 

Munk,  W.  H.  (1949).  "The  Solitary  Wave  Theory  and  Ita  Application  to 

Surf  Problema."  Annala ,  New  York  Academy  of  Sciencea,  Vol.  51, 
Art.  3,  p.  392 

Petera,  A.  S.  (1952).  "Water  Wavea  over  a  Sloping  Beach,"  Communi- 
cationa  on  Pure  and  Applied  Mathematica ,  Vol.  5,  pp.  6'^-l68 

Putnam,  J.  A.  and  J.  W,  Johnaon  (1949).  "The  Diaaipation  of  Wave  Energy 
by  Bottom  Friction."  Tranaactiona ,  Amer.  Geophyaical  Union, 

Vol.  30,  No.  1,  pp.  67^71 

Roseau,  M.  (1952).  "Contribution  a  la  theorie  dea  ondes  liquides  de  gravite 
en  profondeur  variables."  Publications  Scientifiquea  et  Techniques 
du  Miniate  re  de  I'Air,  No.  275,  Paris 

Stoker,  J.  J.  (1947).  "Surface  Water  in  Water  of  Variable  Depth." 
Quarterly  of  Applied  Mathematics ,  Vol.  5,  pp.  1-54 

Stoker,  J.  J.  (1957).  Water  Waves.  Interacience  Publishers,  N.  Y.  p.  319 

Wiegel,  R.  L.  (I960).  "A  Presentation  of  Cnoidal  Wave  for  Practical 
Applications."  Journal  of  Fluid  Mechanics,  pp.  273-286 


LIST  OF  SYMBOLS 


APPENDIX  V 


d 

X. 

X  or  L 
T 

H 

u 

X 

E 

V 

b 

R 

Sub  ^ 
n 


S 

P 


Water  depth 

Wave  length  in  deep  water 
Wave  length  in  shallow  water 
Wave  period 

Elevation  of  wave  surface  above  still  water  level 

Maximum  value  of  ‘q  at  the  wave  crest 

Wave  height 

Horizontal  velocity 

Horizontal  coordinates 

Wave  Energy 

Wave  celerity  for  a  long  wave  in  shallow  water 

Distance  between  two  orthogonals 

Radial  distance  from  the  original  disturbance 
Also  Vertical  wave  run>up 

Related  to  breaking  characteristics 

Exponent  for  variation  of  wave  height  with  depth  vs. 
bottom  slope 

Bottom  slope 

Wave  power  per  unit  width 

Amovint  of  energy  dissipated  by  bottom  friction  per  unit  area 
Shearing  stress  per  unit  area 
Bottom  velocity 
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Depth  on  the  high  side  of  the  bore 
Depth  on  the  low  aide  of  the  bore 
Discharge  due  to  the  moving  bore 
Friction  coefficient  = 

Gravity  acceleration 

Chezy  coefficient 

Manning  coefficient 

Coefficient  characterizing  the  height  of  the  front  of 
a  spilling  breaker 

Breaking  coefficient:  ratio  of  the  energy  dissipated  by 
a  spilling  breaker  and  the  energy  which  can  be  dissipated 
by  a  bore  of  same  height 
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APPENDIX  VI 


THE  WAVE  RUN-UP 

BY  THE  METHOD  OF  CHARACTERISTICS 


By 

B.  L>e  Mehaute 


1.  INTRODUCTION 


This  appendix  presents  the  results  of  theoretical  investigations 
on  the  wave  run-up  which  have  been  carried  out  by  the  method  of 
characteristics. 

This  powerful  method  has  been  of  great  use  in  many  scientific 
fields.  It  has  been  in  common  use  in  hydraulics  for  studying  flood 
routing  and  tidal  waves  in  estuaries.  It  is  also  used  for  studying 
water  hammer  in  pipe  lines,  and  has  been  introduced  in  meteorology 
by  Freeman  (1951). 

It  seems  that  Stoker  (1957)  was  the  first  to  propose  the  applica¬ 
tion  of  the  method  of  characteristics  to  the  problem  of  a  wave  breaking 
on  a  beach.  However,  he  did  not  study  the  problem  of  the  run-up  on 
a  dry  bed.  Hence  this  particular  problem  has  been  solved  in  this  report. 

The  various  topics  of  investigation  are: 

a.  Establishment  of  the  basic  equations  of  the  motion  prior 
to  wave  breaking.  It  is  shown  that  the  vertical  acceleration,  usually 
neglected,  has  an  important  effect.  The  corresponding  correction  term 
is  given. 

b.  The  basic  principle  of  the  method  of  characteristics  is 
recalled.  Dimensionless  parameters  are  introduced, 

c.  The  input  and  limit  of  the  waves  are  defined  by  a  limit 
solitary  wave  where  the  slope  becomes  steeper  than  a  given  value 
such  as  0.  02.  It  is  assumed  that  such  a  wave  on  a  more  gentle  slope 
follows  the  breaking  index  curve  as  demonstrated  in  the  previous  appendix. 


d.  The  problem  of  nonsaturated  breaker  inception  ie 
diicuaeed. 

e.  The  method  for  studying  nonsaturated  breaker  propagation, 
and  bore  propagation,  is  given. 

f .  The  particular  problem  of  the  bore  reaching  the  shore¬ 
line  is  theoretically  analyzed  when  the  bottom  friction  is  neglected. 

g.  The  run-up  on  a  dry  bed  with  bottom  friction  is  also 
theoretically  analyzed  in  various  cases  corresponding  to  various  simpli¬ 
fying  assumptions.  Then  the  links  between  the  motion  over  a  dry  bed 
and  the  bore  reaching  the  shoreline  with  bottom  friction  are  established. 

h.  Finally,  an  application  of  the  method  is  carried  out  and 
suggestions  foi^  further  investigations  are  given. 

2.  THE  BASIC  EQUATIONS  AND  ASSUMPTIONS 

The  classical  equation  for  long  waves  in  shallow  water  neglects 
the  vertical  component  of  velocity,  vertical  acceleration  and  the 
bottom  friction  forces.  They  are:  (see  Fig.  VI- 1  for  notation) 

Momentum:  +  U  ~  9  (VI- 1) 

Continuity:  ^  U  ^ «  q  (VI-2) 

o  I  o  x 
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FIGURE  VI- 1 
NOTATION 


It  is  recalled  that  these  equations  are  based  on  the  assumption  that  the 
distribution  of  horizontal  velocity  u  is  uniform.  Hence  u  is  the 
average  velocity  in  a  cross  section. 

It  has  already  been  pointed  out  by  Stoker  and  others  that  the  method 
of  characteristics,  as  previously  applied,  gives  a  wave  profile  much 
more  unsymmetrical  than  usually  observed  at  sea  or  in  a  wave  flume. 
This  limitation  on  the  validity  of  the  theory  is  of  small  importance  on  a 
comparatively  steep  slope,  but  becomes  very  important  on  a  gentle  slope 
such  as  is  usually  encountered  over  the  Atlantic  continental  shelf. 

The  theory  of  nonsaturated  breakers  (LeM4haut4  1962)  presented 
in  the  previous  appendix,  which  is  based  on  the  assumption  that  the  wa\’’e 
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profile  !•  symmetricali  certainly  gives  a  better  approximation  on  a 
very  gentle  slope.  For  example,  at  the  limit,  if  the  study  of  the 
propagation  of  a  solitary  wave  over  a  horizontal  bottom  is  treated 
by  the  method  of  characteristics,  the  wave  profile  becomes  rapidly 
unsymmetrical  (Fig.  VI-2).  A  bore  appears  even  if  the  starting  wave 
is  not  a  limit  solitary  wave.  (H/d  <  0.78),  However,  it  is  well 
known  by  direct  integration  and  experiments  that  a  solitary  wave  over 
a  horizontal  bottom  travels  without  deformation,  the  wave  profile  always 
being  symmetrical  (see  McGowan  (1891)  and  Munk  (1949)  ).  This  fact 
proves  some  deficiencies  in  the  method  of  characteristics  as  previously 
used.  Ur  sell  and  Birkhoff  (1949)  present  this  as  a  paradox,  the  Earnshaw 
paradox  (1845)*,  also  encountered  in  gas  dynani^s.  Birkhoff  even  states 
that  nobody  knows  the  explanation.  Stoker  (1957),  Laitone  (1961)  and 
others  adso  discuss  the  problems  of  steady  state  in  long  waves.  Stoker 
states  that  the  steady  state  can  be  reached  provided  the  theories  are 
carried  out  at  a  high  order  of  approximation  (p.  342).  Briefly,  the  ex¬ 
planation  of  this  paradox  and  this  deficiency  of  the  method  of  character¬ 
istics  applied  to  long  waves  are  due  to  the  fact  that  the  vertical  accelera¬ 
tion  and  path  curvature  effects  (particularly  important  at  the  crest  of  a 
near-breaking  wave)  are  neglected  and  consequently  the  pressure  distri¬ 
bution  is  assumed  to  be  hydrostatic.  The  solitary  wave  theory  takes 
account  of  this  path  curvature  effect  by  assuming  that  the  vertical  velocity 
is  linearly  distributed  from  the  bottom  to  the  free  surface: 


w 


d+z  dq  ,»,d+z  b9 

TTT  "Ht  *  -mi 
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SOLITARY  WAVE  OVER  A  HORIZONTAL  BOTTOM  BY  THE  METHOD 


Then  it  is  found  that  the  correcting  term  to  be  inserted  on  the  right 
aide  of  the  momentum  Eq.  (VI- 1)  ie: 

-  ji±JL  — (VI-3) 

3  bt'bx 

A  similar  term  can  also  be  added  to  take  into  account  the  influence 
of  the  slope  of  the  sea  bottom  on  the  path  curvatures.  The  length 
of  the  equation  which  is  thus  obtained  is  such  that  its  insertion  for 
practical  computation  is  not  justified.  Moreover,  the  corresponding 
error  is  very  small  for  two  reasons: 

1.  On  a  steep  slope,  the  slope  has  an  effect  on  the  path 
curvature,  but  it  will  be  seen  that  the  total  path  curvature  effect 
is  small  by  comparison  with  the  term  -  gS  which  appears  in  the 
momentum  equation  (S  is  the  slope). 

2.  On  a  gentle  slope  the  total  path  curvature  effect  is  of 
great  importance,  but  the  correction  to  the  path  curvature  effect 
caused  by  the  bottom  slope  becomes  negligible  as  long  as  the  slope 
of  the  free  surface  is  large  by  comparison  with  the  bottom  slope. 

This  is  always  the  case  for  a  near-breaking  wave  close  to  the  crest 
and  at  some  distance  from  the  crest,  the  path  curvature  is  unimportant. 

To  conclude  this  discussion,  the  proposed  correction  term  Eq.  (VI-3) 
can  be  considered  as  a  sufficient  approximation  --  but  also  as  a  neces¬ 
sary  correction  --  for  a  wave  traveling  on  a  very  gentle  slope.  It  is 
pointed  out  that  when  the  bottom  slope  tends  to  zero,  the  solitary  wave 
theory  is  more  exact.  Consequently,  the  theory  previously  presented 
in  Appendix  V  can  be  considered  as  the  limit  case,  obtained  by  a 
direct  integration  where  the  path  curvature  effect  is  taken  into  account. 
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It  haa  bean  pointed  out  that  the  bottom  friction  forcea  have  been 
neglected  in  the  previoua  momentum  Eq.  (VI- 1).  If,  aa  already 
aaaumed  in  Af^endix  V.  the  ahearing  atreaa  at  the  bottom  la 
quadratic  and  the  wave  motion  appeara  aa  a  aucceaaion  of  ateady 
flowa  then  the  Chesy  formtila  givea  for  the  friction  term  F: 

_  g  u|u| 

"cj^  (d  +  V) 

(d  -f  q  )  being  the  hydraulic  radiua  in  thin  cane. 

It  can  be  eeen  that  by  taking  "reaeonable"  valuea  for  the 
bottom  friction  effect  ia  alao  relatively  negligible  when  the  alope  ia 
larger  than  1/10,  but  becomea  very  important  on  a  vary  gentle  alope, 
aa  doea  the  path  curvature  effect.  However,  it  will  be  aeen  that 
the  friction  term  ia  particularly  important  when  the  bore  reachea 
the  ehoreline  and  for  the  run-up  computationa  on  a  dry  bed,  whatever 
the  alope. 

Thene  two  neceaaary  correction  terma  nhow  the  importance  of 
the  notion  of  nonaaturated  breakera  introduced  in  Appendix  V  for 
calculating  the  maximum  posaible  wave  run-up.  Since  nonaaturated 
breakera  follow  the  breaking  index  curve,  it  ia  aufficient  to  atart 
the  computetiona  by  the  method  of  characteriatica  from  a  limiting 
aolitary  wave  on  a  alope  where  the  breaker  tenda  to  be  aaturated: 

S  >  0.02  or  S  s  0.01  for  aafety.  Thia  conaideration  ahould 
permit  elimination  of  a  tremendoua  amount  of  calculation  in  atudying 
the  wave  traveling  over  the  continental  ahelf  becauae  it  would  require 
-taking  account  of  the  path  curvature  term  and  bottom  friction  terme. 
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It  will  be  seen  that  the  method  of  characterietice  apcumuiatea 
errors.  Hence  the  necessary  simplifying  hypothesis  which  has  been 
assumed  to  perform  the  calculation  in  the  energy  method  in  Appendix 
V  does  not  give  a  greater  error  than  does  the  more  theoretically 
exact  method  of  characteristics  when  the  calculations  have  to  be 
performed  over  a  long  distance. 

Finally,  it  is  pointed  out  that  the  basic  momentum  and  continuity 
Eq.  (VI- 1)  and  (VI-2)  assumed  that  the  vertical  velocity  distribution 
was  uniform.  Actually  the  velocity  distribution  of  a  near -breaking 
wave  does  not  satisfy  this  assumption.  As  already  pointed  out,  for 
a  spilling  breaker  the  velocity  near  the  crest  becomes  larger  than 
the  wave  celerity.  Because  of  this  phenomenon,  a  loss  of  energy 
is  not  taken  into  account  by  the  method  of  characteristics;  The 
solution  may  consist  of  imposing,  for  example,  a  maximum  value  to 
the  wave  height:  H  =  0.78d.  This  correction  may  have  some 
importance  if  the  method  of  characteristics  is  begun  on  a  very  gentle 
slope.  But  since  we  already  know  that  the  breaking  index  curve  is 
followed,  this  correction  is  unnecessary  when  we  start  the  method 
of  characteristics  on  a  slope  steeper  than  the  critical  slope  0.02. 
Then  the  spilling  breaker  generates  a  bore  very  quickly  and  the 
correction  due  to  non-uniform  velocity  distribution  is  negligible. 

The  basic  starting  equations  will  finally  become: 


bu 

6t 


+ 


u  ^  »  -gAl 

bx  bx 
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3.  THE  BASIC  PRINCIPLE  OF  CALCULATION  OF  THE  METHOD 
or  CHARACTERISTICS 


Now  soma  elemantary  tranaformationi  muat  ba  parformad. 
Dimanaionlaaa  tarma  ara  alao  introduced. 

Firat  it  ia  aeen  that  the  wave  motion  ia  completely  defined  aa 
a  function  of  time  and  apace  by  the  elevation  nf  (x,  t)  and  horiaontal 
velocity  u(x.  t).  It  will  be  more  convenient  to  define  the  wave 
motion  by  two  terma  u  (x,  t)  and  c  (x,  t).  which  are  both  dimen- 
aionally  equivalent  to  velocitiea,  the  definition  of  c  being: 

C  «  [g(d  +  »>)]i  (VI-6) 


The  definition  for  c  ia  arbitrary  and  can  be  modified  in  an  attempt 
to  eliminate  the  path  curvature  term  in  the  calculationa.  For  example, 
c  can  be  taken  to  be 


C 


(,.  *^)]* 


But  auch  a  definition  for  c  alao  requirea  the  neglect  of  a  few  terma. 

1/2 

In  the  preaent  phaae  of  thia  atudy  the  value  c  =  [g(di-‘7)] 

'frill  be  uaed.  However,  it  ia  intereating  to  note  thia  poaaibility  in 
view  of  further  inveatigationa.  Another  way  of  defining  c  will  alao 
be  obtained  by  aaauming  that  a  aolitary  wave  muat  travel  without 
deformation  over  a  horizontal  bottom. 

From  Eq.  (VI- 6)  it  ia  deduced  that 

It  ■ » -H-  ♦  « ^ 
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Hence 


g 


"  3  Tx*  ’  2c-|7  +  gS 


Also 


g  u  |u  I 
Ch*  (d  +  V) 


and  from  the  continuity  equation: 


bX  b  t 


d-*-  b^d  4  V)  ,  _ltl  [u(d4^)]  _ 
3  bt^bx  3  bx*bt 


(VI-8) 

b^(uc^) 

^g*  bx*bt 


Inserting  these  relationships  into  the  momentun  Eq.(VI-4)  gives: 


^  +  U-^  +c-^  -  G 

bt  bx  b  X  * 


biiy£!i 

bx*bt 


(VI-9) 


In  the  following  the  right-hand  side  of  Eq,(VI-9)  will  be  called 
even  if  friction  and  path  curvature  terms  are  neglected.  The 
continuity  Eq.  (VI-8),  in  terms  of  u  and  c,  also  becomes: 
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b2c 


+  c 


bu 

TT 


(VI- 10) 


Adding  and  subtracting  •uccessively  Eq.(VI-9)  and  (VI-10)  gives 


(u+ 2c)  +  (U +  c)  (u+ 2c)  ■  G*  (VI-11) 

(u  -  2c)  +  (u  -  c)  (u  -  2c)  ■  Q*  (VI-12) 


It  is  seen  that  these  eiqpressions  are  the  total  differential  with  respect 

to  time  of  (u  4-  2c)  and  (u  -  2c)  provided  ^  ~  u  f  c  and 

dx 

^  s  u  -  c  respectively.  Hence, 

^  (u  +  2c)  =  along  a  line  *  u  +  c 

d  dx 

^  (u  +  2c)  =  along  a  line  ^  =  u  -  c 

Before  eiqilaining  the  mathematical  process  based  on  this 
equation,  it  is  particvilarly  convenient  when  performing  the  calculations 
and  also  in  order  to  give  more  generality  to  the  obtained  results  --to 
use  dimensionless  terms.  For  this  purpose,  let  dj  where 

dj  is  an  arbitrary  length.  It  is  most  convenient  to  let  d^  equal  the 
depth  at  the  point  vdiere  the  calculation  of  the  method  of  characteristics 
will  begin.  This  depth  co\ild  be  chosen  as  the  depth  where  the  slope 
(s  0.01,  i.  e.  where  the  breaker  tends  to  become  saturated  and  where 
the  input  can  be  taken  to  be  a  limit  solitary  wave.  Then  inserting  the 
relationships: 
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X  = 


X 


X  d 


1 


t 


X 


T  = 


t 


u  =  C*U 


c  =  C 


the  above  Eq.(VI-ll)  and  (VI- 12)  (where  is  expressed  as  given 
in  Eq.  VI-9)  )  become,  after  some  elementary  calculations: 

_  ^  b^C^) 

■a  »  01^ 

3  bT*bX 

^  C*  b’(UC*)  (VI- 13) 

*  ■J"  SFTr 

adong  ^ijr  -  U  C.  Now  the  principle  of  calculation  based  on  this 

set  of  equations  can  be  explained. 

It  is  seen  that  =  g  G.  Knowing  U,  C,  and  G  at  two 

points  (1)  and  (2)  in  a  T-X  diagram  (see  Fig.  VI-3),  namely  Up 

Cp  Gp  U^,  C^.  G^,  U  and  C  can  be  calculated  at  a  third  point 

(3),  and  by  the  following  processes.  A  line  of  slope 

^ —  is  drawn  from  point  (1)  and  a  line  of  slope  ^ ^ — 
11  2  2 
from  point  (2).  Point  (3)  is  defined  by  the  intersection  of  these  two 

lines.  Then  from  Eq.  (VI-13)  we  know  that: 


A(u 

20  =  G,  ATj_3 

along  the  line 

1 

dT 

Tj-rc7 

di!^ 

A(u  - 

20  =  Oj  ATj.j 

along  the  line 

1 

dT 

dX 

^(UT2C)^G.-S-^(^) 
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Hence. 


and 


U3  ^ 

2C3  =  Uj  +  2Cj 

*  °i 

U3  - 

2C3  =  U2  -  2C2 

+  °2 

fT  — 

.  C  - 

C,  +  G 

.  T. 

”3  ■ 

2  ^1 

^2 

1  1- 

(V 


U,  -  U,  C,  +  C, 

C3=-V-^  ^  T,.2-G2AT2_3 


which  permita  calculation  of  Uj  and  Cj.  Repeating  this  process  of 
calculation  for  each  point  of  the  T  -  X  diagram  permits  calculation 
of  the  complete  wave  evolution  as  a  function  of  space  and  time. 


FIGURE  VI-3 

BASIC  PRINCIPLE  OF  THE  METHOD  OF  CHARACTERISTICS 


3 


The  calculation  of  G  la  tedious  but  not  difficult.  The  friction 


term 


i«  often  small  and  can  usually  be  neglected 
near  breaking  inception.  But,  it  becomes  important  when  the  wave 
reaches  the  shoreline  and  is  essential  for  the  run-up  on  a  dry  bed. 

On  the  other  hand,  the  curvature  term,  while  being  important  near 
breaking  inception,  becomes  negligible  for  the  run-up  on  a  dry  bed. 
Both  of  these  terms  must  be  compared  quantitatively  with  (-S). 


FIGURE  VI- 4 

CALCULATION  OF  CURVATURE  TERM 


The  calculation  of  the  curvature  term  is  complex.  Moreover,  it  is 
effectively  negligible  when  the  slope  is  greater  than  1/10.  It  can  be 
evaluated  from  the  knowledge  of  U  and  C  at  six  points  on  the  T-X 
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diagram  (tee  Fig.  VI-3A)  by  the  formula; 


bX*bT  AX*  AT  ‘■ 

+  2  UjCg*  -  U,C,*  ] 

The  evaluation  of  these  points  generally  requires  a  great  amount 
of  interpolation:  hence,  despite  the  fact  that  it  is  theoretically 
possible  and  despite  its  importance,  it  may  be  theoretically  unrealis¬ 
tic  to  take  this  term  into  account  even  on  a  gentle  slope,  unless  the 
calculations  are  done  by  a  computer. 

These  considerations  again  demonstrate  the  importance  of 
the  study  presented  in  the  previous  appendix  on  saturated  and  non- 
saturated  breakers.  Since  we  know  from  these  computations  that  a 
nonsaturated  breaker  remains  roughly  symmetrical  in  shape  and 
follows  the  breaking  index  curve,  a  great  number  of  computations  are 
saved  by  taking  the  input  on  a  slope  of  0.01  as  has  already  been 
pointed  out. 

4.  THE  INPUT  DEFINITION  AND  WAVE  "UMITS" 

As  already  pointed  out,  the  maximum  possible  wave  run-up  can 
^e  determined  when  the  input  is  defined  by  a  limit  solitary  wave 
.^here  the  slope  tends  to  become  steeper.  According  to  the  classical 
solitary  wave  theory  (see  Munk  (1949))  C  and  U  can  be  obtained 
respectively  from  the  following  equations 

a 
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(VI- 15) 

U  >  V  and  V  ■  [g(d4-H)]? 

0  I  j 

from  which  U  and  C  at  time  T  =  0  can  easily  be  determined. 

However,  it  is  recalled  that  u  as  given  by  this  expression  is  an 

approximate  velocity.  It  is  more  exact  to  define  an  input  by  the 

average  velocity  u(x)  *  — -  /  u(x,z)  dz 

d  +  V 

where 

2  3Q 

I .+  cos  M  -r  cosh  M  -r 

u(x,z)  »  N  _ 2 _ 2.  and  N  and  M  are  given  by 

(cos  M  +  cosh  M 

N  •  -I  .in'  [m  (i  +f  ■7'  [“  ('  +  I)] 

Possible  input  definitions  expressed  as  a  function  of  X  at 
T  =  0  are  represented  in  Fig.  VI-5.  They  are  limit  solitary  waves 
over  a  horizontal  bottom.  They  are  supposed  to  represent 
spilling  breakers  over  a  very  gentle  slope,  reaching  a  steeper  slope. 
They  must  be  expressed  by  dimensionless  parameters  U  and  C 
as  functions  of  X  and  T. 


H 


cosh® 


where 
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T  DEFINITIONS 


It  ic  known  that  the  solitary  wave  has  no  theoretical  horizontal 
limits.  However,  in  practice  U  S  0  and  C  is  very  close  to 
unity  over  a  length  of  X  >  10.  The  friction  and  curvature  terms  also 
tend  to  zero.  Hence  this  can  be  considered  as  defining  the  limit  of 
the  wave.  The  location  of  the  complete  limits  of  the  wave  in  the 
T-X  diagram  and  the  corresponding  value  for  C  can  be  calculated 
by  exact  integration  in  the  case  where  S  is  a  constant. 


At  the  limit 


U - T-C  " 


Hence, 


dT 

dX 


because  U  =  0. 


1/2 

because:  C 


i.  e. 


dT  _ 

17-  ■ 


»  0 


For  example,  in  the  first  case  presented  in  Fig.  VI-4,  the  shoreward 
limit  is  defined  by  the  equation: 


since  X  =  0  for  T  ■  0  .  The  seaward  limit  is  defined  by  the  following 
two  functions:  (1)  on  the  horizontal  bed  (X  <  0) 

dT  I  / <i|  V''* 

dX  •  o’  "UV  “  ■’■'L  +  X  (Vl-18) 

L  being  the  "wave  length"  (L  =  10) 
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(2)  on  the  elope,  it  is  easily  seen  that 
T.L+|.  [l  -d  (VI.19) 

Along  these  limits,  U  =  0  and  C  »  ^  *  (1  -  (VI-20) 

A  number  of  calculations  can  be  eliminated  by  considering  that 
the  slope  does  not  have  any  influence  on  the  wave  motion  as  long 
as  the  characteristic  U  -  C  (of  slope  ^ ^  — =  -1) 

coming  from  the  point  X  ^  0,  T  =  0  does  not  cross  a  characteristic 
U  -f  C.  Hence,  the  input  can  be  taken  along  this  characteristic 
(U  -  C)  (line  AB  on  Fig.  VI-5))  by  a  simple  projection  of  the  values 
of  U  and  C  from  AD  to  AB.  Even  projecting  them  to  AE  will 
introduce  only  a  very  small  error  because  the  effect  of  the  slope  on 
the  wave  motion  remains  small  as  long  as  the  wave  elements  do  not 
encroach  on  the  slope. 


FIGURE  VI- 6 

WAVE  LIMITS  AND  INPUT  DEFINITION 
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Now  the  very  important  phenomenon  of  the  bore  inception  is 
analyzed. 

When  the  two  characteristics  (U  +  C)  cross  each  other  (point  4  on 
Fig.  VI-7)),  then  two  sets  of  values  for  U  and  C  are  obtained  at 
point  4,  namely  and  from  the  equations: 


U  +  2C  =  U,  +  2C,  +  G  A  T,  . 
u  u  I  I  1-4 

U  -  2C  =  U,  -  2C-  +  G  A  . 
i,  u  u  3  3  3-4 

and 


(VI-21) 


(VI-22) 
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Ffiytleally  this  uMtat  IlMt  th«  way*  profile  la  vortical  at  tfaat 
point.  THoro  la  bore  Ineeptlea  and  characterlalng  the  va^e 
0<  U  and  C  just  In  front  of  the  bore  and  and  juat  at  ^e  top 
of  the  front  of  the  bore. 

Then  there  la  a  diacontlnulty  In  the  T-X  diagram.  The  claaalcal 

dX 

method  no  longer  appllea  along  a  line  which  croaeea  the  line  V  » 

(W  s  w  /  gdj,  w  being  the  apeed  of  the  bore). 

Before  atudylng  the  bore  equation  a,  aome  further  conaiderationa 
on  the  bore  inception  muet  be  given. 

It  can  eaaily  be  aeen  diat  all  the  (U  4-  C)  characteriatica  on  the 
front  aide  of  the  wave  converge,  while  all  characterietica  on  the  back 
•  aide  of  the  wave  diverge.  Hence  the  bore  inception  on  the  T-X  diagram 

dependa  upon  the  choaen  interval  A  X  for  the  inpvt  definition. 

» 

The  amaller  the  interval,  the  aooner  the  bore  appeara.  In  fact,  aince 
input  ia  a  limit  aolltary  wave,  it  ia  normal  that  the  bore  begina  aa  ao<m 
aa  the  effect  of  the  elope  ia  felt,  1.  e.  at  the  interaection  of  the 
characteriatica  (U  -  C)  coming  from  T  *  0,  X  »  0.  But  the  loaa 
of  energy  due  to  a  email  apilling  breaker  near  the  creat  ie  negligible 
aa  long  aa  the  initial  interval  ^  X  for  two  characteriatica  (U  4-  C) 
ia  email. 

Along  a  line  croaaing  the  W  line,  the  momentum  and  continuity 

» 

equationa  are  thoae  of  a  ahock  wave,  i.  e.  from  elementary  hydraulica. 

^  (See  Fig.  Vl-8). 
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FIGURE  VI- 8 
BORE  NOTATION 


(d+^)d 


i.  e.  after  inserting  the  dimensionless  notation  U^  ,  ,  W 

-  =d‘  '  '  ("'  •  "d)  ("„  ■  “d)  («-25) 

«u  -  «d  ■  W  (C/  -  C/)  (VI-26) 

from  which  it  is  found  that  the  dimensionless  velocity  of  the  bore  W  is: 


VI.22 


w  .  u^  + 


1/2 


(VI-27) 


Then  a  W  line  from  fhe  bore  inception  muet  be  drawn  in  the  T-X 

die 

diagram  with  a  elope  such  as  — —  =  W.  It  must  be  noted  that 

this  line  is  between  the  line  of  slope  and 

(see  Fig.  Vl-9).  But  W  also  varies  as  a  function  of  X  .  Hence  it 
must  be  computed  step  by  step 


4 


/ 


FIGURE  VI-9  FIGURE  VI- 10 

BORE  CHARACTERISTICS  ON  THE  LOW  SIDE  OF  THE  BORE 

AFTER  INCEPTION 


FIGURE  VI- 11 
INTERPOLATION 


FIGURE  VI- 12 
EXTRAPOLATION 


V 


First  it  must  bs  noted  that  U  and  C  can  always  be  calculated 
by  the  classical  process  on  the  low  side  of  the  bore.  i.  e.  in  front 
of  die  bore  because  W  is  always  greater  than  (U  i-  C)  on  that  side  (see 
Fig.  Vl-10).  But  point  (3)  of  Fig.  VI-10,  calculated  from  points 
(1)  and  (2).  does  not  necessarily  coincide  with  any  point  of  the  W  line. 
Hence  and  Just  in  front  of  the  bore  must  be  determined  either 
by  interpolation  between  (2)  and  (3)  along  the  line  (U^,  -  (see  Fig. 
VI- 11)  or  by  extrapolation  from  (2)  to  (3)  along  a  line  (U^  -  C^) 

(see  Fig.  VI- 12). 

This  work  is  greatly  simplified  when  the  W  line  crosses  the 
"limit"  of  the  wave  (defined  by  Eq.  (VI- 17)  because  then  the  bore 
travels  on  still  water  with 


1/2 

(l-SX) 


(VI-28) 


Now  three  other  unknowns  remain:  namely  .  and  the 

W  values  for  the  following  step.  Hence  three  equations  are  necessary. 
These  three  equations  are  given  by  the  continuity  and  momentum 
for  the  bore  and  have  been  already  expressed  in  (VI-25)  and  (VI-26). 
Another  equation  is  given  from  the  classical  relationship 


d  (U  +  2C)  »  G  dT 

applied  along  the  (Ug  4-  Cg)  line  on  the  high  side  of  the  bore  and 
which  crosses  the  W  line  at  the  point  to  be  determined.  (See  Fig. 
VI- 13). 


VI-25 


X 


FIGURE  VI- 13 

ON  THE  HIGH  SIDE  OF  THE  BORE 

For  this  purpose  the  determination  of  the  values  Uq  and  Cq  of 
point  (0)  can  be  done  theoretically  on  the  (U^  -  C^)  line  by  inter¬ 
polating  the  U  and  C  values  between  (1)  and  (2).  Then 


■  “o*“o*°*  Vu 

In  practice  points  (0)  and  (1)  are  so  close  to  each  other  that  the  values 

for  Un  and  C.  can  sometimes  be  taken  as  U  and  C 

0  0  "l  “l 

Then  Eq.  (VI-29)  is  replaced  by: 


VI.26 


I 


* 


Now  *04  W  can  ba  daterminad  from  Eqa.  (Vl-S5)«  (Vl-26), 

and  {VI-29)  or  (Vl-JO). 

Tha  ■olutioa  of  tiiia  ayatam  ia  givan  by  tha  following  aat  of  a^piattona: 


where 


y^.i 


2Y  [l  +  Y^] 


+  2Y  »  K 


C 

Vs  ** 

y  a  ^ 


..  (v^^o)  -  v« 
^  - 

“0+“0  “  (“«♦“„) 


(VI-31) 


(VI-32) 


(VI-33) 


Tha  function  K  »  f(Y)  haa  bean  drawn  for  a  range  of  poaaible  values 

for  Y.  Then  Y  is  determined  graphically  from  the  enclosed  curve 

Fig.  VI- 14.  If  Eq.  (VI-31)  is  mathematically  solved,  then  six  values 

for  Y  can  be  found,  but  only  one  has  a  physical  significance.  This 

C„ 

value  can  easily  be  known  since  Y  -  -  always  varies  slowly  along 

^d 

tha  W  line. 

When  is  calculated,  then  is  easily  obtained  from 
Eq.  (VI- 30).  Then  W  can  be  calculated  for  the  following  interval  by 
formula  (VI-27). 


VI-27 


100 


FIGURE  Vl-14 

GRAPHICAL  REPRESENTATION  OF  THE  FUNCTION 
C 

•  iiK)  FOR  BORE  CALCULATIONS 
VI-28 


ts  aMlyMkU  Thi*  yM»bUm  atvAMUlir  riM  to  «  aiwnlwr 
o£  tntorttstiim  ttedlM  of  a  matoainaatical  natara:  Wliitoam  {l^Hh  KaUar* 
at.  al.  (I960),  Hto  and  Mayaz  (1962).  Tha  followiag  ealctitoHotta  aza  . 
tevalopad  from  a  dlffarmt  approach  but  give  aimilar  raaalta.  A 
caratol  eompariaoa  of  thaaa  varioua  atudiaa  hat  not  haaa  emiaidartoi  aa 
baiag  witoin  tha  acopa  of  toia  stadp  bacauaa  of  toa  ra^ar  acadamic 
nature  of  the  reaearch.  It  ariU  be  aeen  in  toe  foUowiag  aeetion  that 
for  practical  parpoaea  tha  botUmt  friction  foreaa  comfdately  ehaaga 
tha  reaulta  preaanted  in  diia  aaction. 

Conaider  the  bore  formula: 


W 


(VI-34) 


It  ia  aeen  that  —  0  near  the  ahoreline,  and  at  the  aame  time 

^ .  d  [  ^ri  -  sxl 

”331 -  *  - - *  -  -  00 

Hence  the  variationa  of  W  ^  and  alao  and  muat  be  expaetad 

to  be  very  large  near  the  ahoreline.  The  atep-by-atep  proceaa 

calculation  baaed  on  the  aaaumption  of  alow  variationa  ia  no  Im^ar 

valid  as  the  ahoreline  ia  approached. 

From  formula  (VI-34),  it  ia  seen  that  vrhwi  •>*  0,  W  cm 

Itova  any  value  dapanAng  uptm  the  value  for  .  In  fiact,  it  hM  ftoiai 


Vl-29 


{bund  by  Trmvnmn  that  W  can  never  eMaaid  ^  .  Indeed  the 

wave  elemente  which  yi  eld  their  energy  to  the  bore  arrive  at  a 
epeed  U  -t-  C.  The  diesipation  of  energy  due  to  the  bore  cMtcee 
to  decrease  in  such  a  way  that  the  speed  of  .the  bore  W  also 
decreases  to  a  value  -t-  after  which  more  energy  is  provided 
from  behind  the  bore.  Hence  the  Freeman  limit  W  < 
shows  that  W  is  finite  if  and  are  finite.  It  can  easily  be 
seen  from  characteristic  equation  along  any  characteristic 


2C J  =  (U  +  2C)  +  GA  T  =  Cst 


(VI.35) 


that  and  can  never  be  infinity.  (They  are  both  positive). 
Moreover,  consider  the  relationship  (VI>34)  in  which  W 

is  finite  and  C  .  0.  It  is  seen  that  this  can  be  achieved  only  if 

^  C 
^11 

both  C  -*0  and  -  — •  oo  .  Hence  C  .  must  be  an  infinitesimal 

u  t^d  ® 

of  higher  order  than  C  .  The  Eq.  (VI-35)  becomes,  at  the  shoreline: 


U  =  Uc  =  U  +  2C  +  GA  T  (VI-36) 

u  o 

This  means  that  the  bore  (or  shock  wave)  disappears  at  the 
shoreline  and  that  the  potential  energy  (function  of  )  is  suddenly 
completely  transformed  into  kinetic  energy  (function  of  ).  More¬ 
over,  by  continuity  it  must  be  expected  that 


W  --  (VI-37) 

which  is  in  agreement  with  the  Freeman  limit.  It  is  seen  also  diat  the 
Freeman  limit  is  reached  only  at  the  shoreline  where 
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W  —  II4^C  -•  UoU 

tt  »  tt  ,  * 


Hence  et  the  ehoreUne  tfie  chnrecterietic  eqoetion  aui  be  written: 


d;(u^+2C^)  ■  -  S  dT 


S  ^ 


or 


Letting  e  ^ 


and  <VI-37)  give: 

W 


U 


and 


$ 


[-^1' 


/2 


(VI-38) 

(VI-39) 
Eq«.  (VI- 34) 


(VI.40) 


Dividing  (VI-38)  by  ,  inserting  o  and  /9  ,  and  the  value  for  a 
given  by  (Vl-40)i  then  expressing  a  as  a  function  of  /9  and 
as  a  fxmction  of  S  and  X  ,  it  is  found  that: 


$ 

$ 


14/9^ 

“r“ 


77^ 

— r“ 


1/2 


ITT 


1/2 


+ 


TTVT 


+  i 


SdX  .  _ 

■nrgjr-  +  ® 


(VI-41) 


-2 


0  is  always  larger  than  unity  (and  even  tends  to  infinity)  hence,  the 
denominator  of  the  left  hand  term  is  positive.  Hence  Eq.  (VI-41) 


implies  that 


0 


-•  00  I  since  0 


(■^) 


already  tends 


to  infinity  for  W  to  be  finite,  d  0  must  also  tend  to  mfinity  but  at 

c  c 

a  higher  rate.  Consequently,  since  ^  -  00 ,  then 

also  -  00  but  at  a  much  greater  rate. 

Moreover,  from  the  characteristic  Eq.  (VI-38),  it  is  seen  that 
since  -  -  00  .  —  +  00  •  Accordingly  —  4 


00  • 
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To  sum  up  theao  rooult*.  near  the  choreline  If  Ae  bottcnn 


friction  forces  are  neglected: 

Ud  -  0 

Cd^o 

C  —  0 
u 


W  —  U  :  finite 


d(Ud) 

~T3r“ 


C  0 


d(Cd) 


00 


dW 


d  (U  ) 
u' 

ssr 


00 


The  bore  disappears  at  the  shoreline  and  is  replaced  by  a 
"rarefaction  wave"  which  appears  as  an  edge  of  water  climbing  over 
a  dry  bed. 

at  the  shoreline  (called  in  the  following)  can  then  be 

approximated  very  simply.  Since  it  is  known  that  at  the  shore- 

dX 

line  tends  to  zero  and  that  “yp"  =  W  -•  +  C^,  can  be  given  by: 


U.  =  (U .  +  2  C  )  +  G  ^  T  (VI-42) 

In  practice  the  values  and  can  be  taken  from  the  step  prior 

to  the  quick  variations  of  W,  and  C^.  For  example,  the  point 

dW 

could  be  chosen  where  “gjf-  =  0.  In  practice,  also,  the  rate  of 
variation  which  follows  is  so  quick  that  if  G  A  T  is  taken  to  be 
-S  AT,  it  is  found  to  be  negligible  in  Eq.  (VI- 42)  and  this  equation 
may  be  taken  quite  simply  as  +  2C^.  This  equation 


VI-32 


•xplalna  th«  audden  Ineraaaa  of  water  velocity  near  the  ahoreline 
when  the  ahock  wave  or  bore  diaappeara. 

7.  THE  WAVE  RUN-UP  CALCULATION  ON  A  DRY  BBD 

At  the  extreme  edge  of  water  C  »  0  ,  becauae  H  »  -d. 
Hence,  the  characteriatic  equationa  ia  simply: 


dU  >  GdT  along  ^ 


(VI-43) 


Hence  d 


(^) 


-S  dX  and  neglecting  the  friction  forcea: 


T*  ”  "T"  “  S(X-X^) 


(VI -44) 


where  X  is  the  shoreline  coordinate. 

s 

The  vertical  run-up  can  be  calculated  directly  in  the  case 
where  it  is  the  first  wave  element  which  has  the  most  energy 


(VI-45) 


However,  the  following  wave  elements  usually  overtake  the  first 
one,  so  increasing  the  run-up  and  generating  a  "roll-wave."  The 
exact  calculation  of  the  run-up  is  then  difficult  because  at  the  front 
of  the  wave:  U  +  C  =  U  -  C  .  This  difficulty  is  automatically 
solved  if  the  friction  forces  are  taken  into  account  as  follows; 


It  is  seen  that  if  at  the  front  of  the  wave  C  -•  0  ,  the  friction 

I  U  V 

term  f  t  p'  )  tends  to  infinity.  Hence  the  edge  of  water  is 


cut  short  and  the  leading  front  appears  as  an  almost  vertical  wall 
of  water  which  has  the  physical  appearance  of  a  bore.  Then,  the 
basic  equation  for  the  leading  wave  element  is: 
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d(U-t-  2C)  > 


-[••■(s) 


dX 

u  +■  c 


(VI-46) 


For  reason  of  similitude,  there  exists  between  U  and  C  a 
linear  relationship:  U  =  A  C.  The  value  of  A  depends  upon  the 
friction  coefficient.  This  value  can  be  determined  by  anology  with 
a  boundary  layer  problem.  If.  more  simply,  one  assumes  that  at 
the  front  of  the  wave  all  particles  have  the  same  velocity  u  ,  the 
convective  inertia  is  zero.  Also,  the  local  inertia  p  ^  is 
negligible  by  comparison  with  the  pressure  gradient  and  friction 
term.  Hence  by  equating  pressure,  gravity,  to  bottom  friction 
forces,  it  gives 'F  ig,.  VI-15. 


(  d  +  17 )  [  -  -  S  ]  -  u^ 


(VI-47) 


FIGURE  VI- 15 

WAVE  FRONT  ON  A  DRY  BED 
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AlsOt  the  elope  S  Is  smell.  Hence,  by  integrating  Bq.  (Vl-47) 

«>+’)*  ■  -T  (S-*) 

being  the  coordinate  at  the  extreme  front  of  water.  It  is  seen 
that  the  front  of  water  is  in  the  form  of  a  parabola. 

Also,  for  reason  of  similitude,  the  length  of  the  front  of  water 

2  2 
(  )i  proportional  to  u  such  as  (  )  =  Bu  , 

Hence  =  2  f  B  U**  and  A  *  (2fB)^^^ 

An  adaptation  from  a  work  of  Keulegan  (1949)  gives  A  ^  1/2. 

Whatever  the  value  for  A  .  Bq.  (VI-45)  can  be  written 

(1  +  A)(1+2A)  ^  +  s)  -0  (VI.48) 


which  gives  for  U; 


(VI-49) 


The  leading  front  characteristic  is  determined  graphically  step  by 

dX 

step,  by  the  intersection  of  the  curve  -jpjr-  =  U  where  U  is 
given  by  the  Eq.  (Vl-49)  and  the  (U  +  C)  characteristics  catching 
up  the  leading  front.  The  value  of  U  at  the  intersection  is 
determined  along  this  characteristic  ( U  i-  C  )  from  the  equation: 


(Uo+2Co)  -  (l  +  2A)U  +  GAT  (Vl-50) 

The  wave  front  becomes  increasingly  insensitive  to  the  following 

dX 

wave  elements.  The  maximum  run-up  is  found  vdien  jy  -»  0, 

for  all  characteristics  (  U  -t-  C  ). 
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Now  the  problem  of  the  bore  near  the  ehoreline  with  btfttom 

friction  can  be  solved.  It  has  been  seen  that  when  C  .  0, 

dC  ** 

“•  0,  W  -•  U^,  and  — jy —  “  00  •  1“  f*ct  such  limits 

are  unrealistic  because  they  are  obtained  by  neglecting  the  bottom 
friction,  which  is  never  negligible  near  the  shoreline.  Hence,  in 
practice  the  theoretical  curves  shown  in  Fig.  VI- 16  should  be 
replaced  by  more  factual  curves.  These  factual  curves  are  obtained 
by  considering  the  relationship  which  exists  between  U  and  C 
at  the  front  of  the  run-up  on  a  dry  bed.  This  relationship  has  been 
previously  established  and  depends  upon  the  friction  coefficient  f. 

Hence,  as  soon  as  the  relationship  AC^  is  verified,  the  motion 

must  be  considered  as  a  motion  on  a  dry  bed.  The  bore  theory  is 
no  longer  valid. 


- WITH  BOTTOM  FRICTION 

FIGURE  VI- 16 

BORE  CHARACTERISTICS 
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8.  APPLICATION 


Th«  method  described  in  the  previous  section  has  been  partly 
applied  by  R.  8.  Grewal  over  a  beach  of  constant  1/10  slope,  corres¬ 
ponding  to  one  of  the  highest  beach  slopes.  Since  it  was  the  purpose  of 
this  study  to  find  a  method  of  computation  for  the  run-up,  such  a 
steep  slope  was  chosen  because  it  permitted  investigation  of  all  the 
various  theoretical  aspects  of  the  problem  with  less  calculation  than 
a  more  gentle  slope.  Also,  since  the  path  curvature  is  very  tedious 
to  calculate  and  does  not  present  any  theoretical  difficulties,  it  can 
justifiably  be  neglected  for  a  slope  of  1/10.  The  bottom  friction  term 
is  also  negligible  except  near  the  shoreline. 

For  sufficient  accuracy,  at  least  four  characteristic  numbers 
are  necessary  to  perform  the  calculations  because  errors  accumulate 
rapidly. 

The  input  definition  has  been  calculated  as  indicated  in  section  4. 
The  results  are  presented  in  Table  VI- 1  for  0  >  X  >  5  .  They  are 

symmetrical  around  X  »  -5  for  -5  >  X  >  -10  . 


TABLE  VI-  1 
INPUT  DEFINITIONS 
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The  corresponding  results  ere  presented  in  Fig.  VI- 17.  Fig.  VI- 18 
is  an  enlarged  portion  of  Fig.  VI- 17  near  the  shoreline.  Fig.  VI- 19 
represents  the  successive  wave  profiles  obtained  from  Fig.  VI- 17 
and  VI- 18. 

Because  such  results  have  a  definite  value,  the  values  of  W  , 

U  and  C  for  all  the  points  are  given  in  Tables  VI-2  and  VI-3  and 
Figures  VI-20  and  VI-21  give  the  values  of  , 

W  and  ‘>1  as  fvinctions  of  X  and  T  respectively. 

Now  the  r\m-up  can  be  calculated  from  the  following  set  of 
eqviations.  Since  G  Zi  T  is  small,  U  and  C  at  the  shoreline  are 
obtained  from 

U  +2C  =U  +2C 

8  8  U  U 

Taking  the  last  obtained  values  for  and  at  X  -  9.  52 
(the  shoreline  is  at  X  =  10  ),  one  obtains: 

Ug  +  2  Cg  =  1.223  +  2  (.635)  =  2.493 
And  from  the  relationship  =  A  where  A  is  arbitrarily  chosen 
as  0.  34  for  example,  it  is  found  that  U  si.  568,  and  C  =  0.  466. 

The  application  of  the  method  presented  in  section  7  has  not  been 
performed  quantitatively  because  it  was  beyond  the  scope  of  the  present 
project.  The  Fig.  Vl-22  gives  a  qualitative  aspect  of  the  characteristic 
method  over  a  dry  bed.  Some  quick  calculations  seem  to  show  that 

D 

■n  ~  Z.  82.  In  fact  the  final  result  depends  upon  the  value  for  A. 
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FlGURi:  Vl-17 

APPUCATION  OF  THE  METHOD  OF  CHARACTERISTICS  OVER  A  1/10  BOTTOM  SLOPE 


9.4 
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.n* 


FIGURE  VI- 20 

BORE  CHARACTERISTICS  VERSUS  DISTANCE 
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BORC  CHARACTERISTICS  VERSUS  DISTANCE 
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BORE  CHARACTERISTICS 
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FIGURE  Vl-22 

QUALITATIVE  ASPECT  OF  THE  METHOD  OF  CHAR¬ 
ACTERISTICS  FOR  STUDYING  SURGE  OVER  A  DRY  BED 
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It  Bamia  that  th«  viiua  for  A  is  botwssn  ths  valuo  proposed 
by  Ksttlogan,  namoly  A  «  0. 5.  and  the  value  chosen  in  this  example. 
A  «  0.34 

9.  CONCLUSION 

A  theoretical  method  has  been  developed  for  computing  the  wave 
run>up.  The  various  steps  of  computation  have  been  described, 
namely: 

a. .  Input  definition  for  possible  maximum  run-up 

b. .  Wave  deformation  on  a  slope 

c.  Wave  breaking  inception 

d.  Bore  traveling  within  the  "limits"  of  the  wave  (U^  0). 

e.  Bore  traveling  outside  the  "limits"  of  the  wave  (U^  =  0). 

f .  Bore  reaching  the  shoreline  without  bottom  friction 

g.  Bore  reaching  the  shoreline  with  bottom  friction 

h.  The  run-up  on  a  dry  bed  with  bottom  friction,  including 
the  determination  of  the  shape  and  height  of  the  wave  front. 

It  can  therefore  be  said  that  a  great  step  forward  has  been  made 
in  afield  which,  despite  its  importance,  has  had  a  relatively  small 
number  of  theoretical  studies.  Most  previous  studies  on  wave  run-up 
have  been  mainly  limited  to  experimental  studies  on  rather  steep  slopes 
in  laboratories.  Previous  theoretical  studies  on  the  climb  of  a  bore 
on  a  beach  mainly  covered  topics  e  and  f  above.  Practical  considera¬ 
tions  have  shown  that  the  bottom  friction  modified  the  theoretical 
results  obtained  in  f  . 
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Deapite  the  preeent  echievemente,  the  problem  deserves  further 
investigation  and  refinement.  A  computer  program  has  to  be  developed 
for  analyzing  rapidly  and  accurately  the  cases  of  waves  traveling  on 
various  slopes.  The  friction  coefficients  can  be  determined  by 
correlating  the  present  theories  with  some  experimental  results. 

Then  other  cases  can  be  investigated  rapidly. 

Finally  the  theory  of  run-up  on  a  dry  bed  can  be  refined  by 
further  analytical  and  experimental  investigation. 

Among  other  topies  to  be  investigated  prior  to  writing  the 
computer  program  is  the  choice  of  interval  and  characteristic  numbers 
as  a  function  of  the  error.  Also  a  numerical  method  can  be  work  out 
to  calculate  U  ‘and  C  for  given  values  of  X  and  T  defining  a  fixed 
rectangular  net  as  it  has  been  proposed  by  Stoker  for  flood  routing. 

This  process  will  permit  the  path  curvature  term  to  be  taken  into 
account  more  easily. 
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UST  OF  SYMBOLS 
APPENDIX  VI 


X 

z 

t 

d 

V 

u 

w 

e 

g 

Ch 

£ 


X 

T 

U  = 

C  = 


Horizontal  coordinate  positive  towards  the  shoreline 

Vertical  coordinate  positive  upwards  (Mid-water  level 
z  s  0) 

Time 

Depth  (bottom:  z  =  -d) 

Vertical  distance  from  the  free  surface  to  the  mid-water 
level. 

Horizontal  water  velocity 
Vertical  water  velocity 
.y  g  (d+  1?) 

Gravity  acceleration 
Chezy  coefficient 


Friction  coefficient 


2,2  2 

Bottom  slope 


2  . 

b  (  u  c  ) 

b  X  ^ bt 


>/g  ^  1 

Depth  at  the  origin 


~T~ 

u 

viar 


Dimensionless  horizontal  distance 

Dimensionless  time 
Dimensionless  horizontal  velocity 
Dimensionless  velocity  for  wave  elements 


VI- 51 


G  * 
H 

V 

Sub^ 

Subd 

Sub^ 

vr 

W 

Y 

K  = 


a  = 

0  = 


g  :  Dimenaionlea*  coirecting  term  for  character ifltlcs 
Wave  height  of  a  solitary  wave 
Solitary  wave  velocity  (  (d  +  ti)  ) 

Related  to  the  high  aide  of  the  bore 

Related  to  the  low  aide  of  the  bore 

Related  to  the  shoreline 
Related  to  the  front  of  the  wave 

Bore  velocity 

:  Dimensionless  bore  velocity 


(U  ^  +  2C)  -  Ud  +  G  At 


C 

u 

-C- 

d 


R 

R*  = 

A 

B 


Vertical  Run-up 
R 

U 


Dimensionless  vertical  run-up 
Coefficient  for  the  wave  front 

Coefficient  characterizing  the  length  of  the  wave  front. 
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APPENDIX  Vn 


TWO-DIMENSIONAL  NONUNEAR 
WAVE  MOTION  IN  AN  ESTUARY 


By 

John  C.  Freeman,  Jr. 
Larry  Armijo 
Mary  Ann  Noeer 


1.  INTRODUCTION 

The  anelyeie  of  two-dimeneional  waves  in  ah  estuary  has  re¬ 
sulted  in  two  studies: 

a.  The  nonlinear  motion  of  a  single  finite  wave.  This  has 
yielded  some  new  developments  in  wave  analysis. 

b.  The  presentation  of  a  finite  difference  scheme,  suit¬ 
able  for  numerical  machine  computation,  that  allows 
for  evaluation  of  two-dimensional  wave  motion  by  the 
method  of  wave  derivatives.  Although  other  approaches 
m'ay  exist,  this  was  selected  as  the  best  and  most 
economical  computation  method  for  the  wave  problem 

at  hand. 

2.  NONLINEAR  MOTION  OF  A  THIN  WAVE  THROUGH  A  BAY 

A  moving  wave  is  not  a  general  two-dimensional  disturbance 

over  a  40  mile  square.  Indeed,  even  a  wave  with  a  period  of  a 
hundred  seconds  influences  only  about  a  mile  and  most  indications 
are  that  the  wave  passes  over  an  area  in  a  finite  time  and  leaves 
the  given  area  relatively  undisturbed  after  passage.  This  implies 
that  the  wave  is  a  relatively  thin  long  disturbance  moving  through 
a  bay  that  is  certainly  undisturbed  ahead  of  it  and  probably  undis¬ 
turbed  behind  it.  Thus,  a  solitary  wave  may  be  regarded  as  a 
disturbance  influencing  only  a  narrow  strip  at  any  given  time. 
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We  write  the  equations  of  motion  in  shallow  water  in  the  following 


form: 


(vn-1) 


)y  /  ~  vox  oy 

In  order  to  discuss  the  problem  of  leaving  the  fluid  undisturbed,  we 
consider  a  one -dimensional  problem. 


This  give  us  — ‘s  2  Co 


-E, 


(VII-2) 


along  =  U  +  C 
d  t 


which  is  the  way  we  follow  the  wave 


d  (u  -  2^  s  +  2  C  •“— **  "  E  along  =  u  - C 
— g-j  bx  *  dt 


The  wave  moves  along  the  x  axis  with  a  speed  of  the  order  of  c  so 
<lx 

that  ■  =  u  -  c  moves  through  the  wave  with  a  speed  of  the  order 

of  -2c.  Thi?  means  that  the  line  — ^ —  =  u  -  c  is  almost  parallel 

,  jj  d(u  -  2c)  _  dt  d(  u  -  2c) 

to  the  X  axis  and  we  can  consider  - — '  =  -g^'  — i— - . 


If  we  neglect  friction,  _^i^££L 


Co 


hCo 


.  Since 


-Co--.  «- 


and 

u  -  c 

parallel  to  t  =  zero, 
d  (  u-  2c) 

— asj - 


dx 

“ar 


U  -  C  i)  X  .  ' 

dx  =  (u  -  c  )  is  almost 


-St 


6(u  -2c) 
S"x 


b(-2  Co) 
STx 
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This  tells  us  that  a  good  approximation  is 


(u  -  2c)  (x,  t)  =  -2  Co  (x). 


This  approximation  includes  as  a  result  that  the  water  is 

dx  , 

—  =  u  +  c  passes 

u  __  Co 


undisturbed  after  the  last  characteristic 
a  point.  Note  that  this  says  u  -  c  =  Co 


so  that 


-Co  - 


u 

T 


1  *f 

rs^  *  ' 


U 


The  thin  wave  and  the  observed  fact  that  most  of  the  motion  and 
wave  height  changes  occur  along  a  normal  to  the  wave  augg«sts  the 
desirability  of  a  system  in  which  a  normal  to  the  wave  is  used  as  a 
coordinate  system  and  in  which  the  largest  fluid  motions  are,  parallel 
fo  the  wave  so  that  v  ~  can  be  considered 

to  the  first  order  only,  u  and  c  will  be  considered  to  higher  order. 

Since  changes  in  the  y-direction  in  the  ordinary  flow  would 
move  through  the  wave  much  faster  than  along  it  (See  .‘tl-  1) 

it  is  assumed  that  the  chief  contribution  to  v  would  be  the  turning 
of  the  wave  front  as  you  move  parallel  to  it. 
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FIGURE  VU-1 

BASIC  PRINCIPLE  OF  CALCULATION 


If  the  wave  travels  from  A  to  B,  a  disturbance  at  A  could 
have  only  affected  the  wave  in  the  shaded  area  and  that  only  while  the 
wave  travels  to  C  .  The  shaded  area  is  determined  completely 
by  condition  between  D  and  B  in  the  time  it  takes  to  travel  from 
O  to  B  .  These  considerations  lead  to  the  approximation  that 

V  (y  -  yj  =  U(y  -  y^)  B{y  .  y^) 

or  (u  810(9**^)  in  extreme  cases). 

Q 

The  term  C  in  Eq.  VII-1  can  then  be  written  C  ■^(u^)  . 

We  will  ultimately  find  that  W  (the  speed  of  a  forward  advancing 
bore)  will  control  our  wave  and  that 


which  we  will  show  later. 
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This  leaves  us  with  the  following  equation  fot  the  motion  of  a 
wave  with  a  leading  bore  moving  through  a  shallow  bay  in  a  coordinate 
system  with  x  normal  to  the  wave. 

c  * 

-gf ;  V  •  f  K.c,.cJ  ;  c^«  C,(x^.y^, 

Xb  .yb  1^1,  are  coordinates  of  a  leadiif[  bore. 

This  give  us 

along  =  u  +  c 


We  are  also  given  as  an  approximation  that 


u  -  2  C  *  -2 
or 

u  =  2  (c  -Co) 


So  we  have 


.1 


^(4c-2c,)=2Co^-c^2(c-Co)0  *4  | 


hg*  (c-Co)  |c-c,j 


along 

This  can  be  written 


Jf-3c-2c. 


f  c  .40 ^  ^  ,c-c,e- 
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along 


Note  that: 


(3c~2co)-Co 

3 


w‘  *'■*‘0 

(3c-2co)  +  2co 

g  .  i -  ggd  c  -  Co  ■ 

are  handy  equations  to  use  with  this  differential  expression. 

We  can  rewrite  this  expression  setting  t>  3c~'2c,  so  that 

«•  -  TC  +4Co»  - -|-C  -  ^Co  +  4Co  •  -  •|'C  + 

2.  i-c*  T  '•-Co 

3.  (c-Co)  I  c  -Co|- 

4  ^2/3  ,  ^JLtZCo  ^2/3 

3 


The  resulting  Statement  of  the  problem  is  as  follows: 


<*«  .  I  s  ,  1  s  ,  (a  "‘"I  6  n  I,  r  \ 

dT  +  T'=o>  TT - 6?  ® 


^^“1/3  (C~Co)  *~Co 
3  ■  C  1 2Co  .2/3 

^  3  ' 

along 

dt 

It  is  obvious  from  the  structure  of  this  set  of  equations  that  if 
a  system  in  which  we  have  large  values  of  C  following  on  smaller 
values,  that  the  larger  values  will  overtake  the  smaller  and  a  jump 
(or  a  bore)  will  form.  Thus  we  are  most  interested  in  the  wave  with 
a  leading  bore. 

We  calculate  the  conditions  at  and  the  speed  of  this  bore  with 
the  equations  of  conservation  of  mass  and  of  momentum.  Consider 
the  bore  as  stationary  with  the  fluid  moving  through  it  at  speed  -V. 
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% 
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FIGURE  VU-2 
BORE  NOTATION 


Continuity  of  ma«a  tells  us  -V  *  (“V  +  u  )  {dj^+  t;^)  . 


From  our  previotisly  given  equations,  we  find 


V  Co«  [-V  +  2(C-C,)j  C*, 
V(c*-c5)  -2  C*  (c-Cp), 


or  the  speed  of  the  bore 


V  « 


2c^ 

C+Co 


We  express  this  in  terms  of  C 


C  » 


C  +  C, 


~r*T^ 


’  T  +  I 
2 


(«+2c.) 


(VII- 3) 


It  is  inherent  in  the  physics  of  the  problem  that  C  -  V  ^  0  .  If 
C  —  V  *  0  or  nearly  so  and  there  is  an  indication  c,  changes  so 
that  ~y  <  0  then  we  must  adjust  a  so  that  d(c-V)BO. 

(This  usually  occurs  with  decreasing  Cq  . )  From  Eq.  (VIl-3), 

>  W  except  for  C,  .  There  is  a  following  value 

of  C  <  which  we  can  say  c  >  V*  This  value  of  c  will 
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be  moving  at  the  same  speed  as  the  bore  and  the  values  between  it 
and  the  bore  will  be  trapped  unless  the  bore  slows  down. 

Over  a  large  part  of  the  course  of  a  wave,  you  might  esqiect 


C  **  .  We  write: 


2  r*^+5cCo  .  4Co -CC,  1 

V  •  T  I  c+Sco  ctSco  J 


4c|-CCe 


2  .  2  ^Cq-c 

T-*  +  T"e  1+5^: 


and  get  as  an  approximation: 


v«  I-. +^. 


We  are  now  in  a  position  to  compute  the  wave  and  the  bore. 


■r-dT'^T'o-br 


( .-Co)  Ic-cqI 
3  «-f2c.  .Z/i 


In  order  to  speed  up  computation,  we  now  make  two  assumptions. 
One  is  that  we  can  recognize  when  (  C  ~C,  )  is  negative  or  positive 

so  that  we  can  simplify  the  last  term  in  the  equation;  the  other  and 
less  justified  assumption  is  that  (  C  —  Co  )  has  a  weak  dependence  on  y 


so  that  the  second  term  can  be  written  . 


Co  ,  ,  b® 

6 


(  c-cq)^  be 

6  by 


We  now  set 


bt  ^  bx 


and  we  have 


_L  dCfl 

4  dt 


3  d'co 

T-ir 


Ce(c-Co)  be  (<--^Co)^rbe  .  1 

6  by  ^  IW  (c+2Co)2'5  J 
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be  moving  at  the  tame  epeed  as  the  bore  and  the  values  between  it 
and  the  bore  will  be  trapped  unless  the  bore  slows  down. 

Over  a  large  part  of  the  course  of  a  wave,  you  might  expect 


« 


We  write; 


V  = 


2 

T 


[- 


+  5cCo 

c  +5co 


+ 


4c|-CCo  1 

e+5co  J 


V  » 


+ 


3  Tfe; 


and  get  as  an  approximation: 


V« 


We  are  now  in  a  position  to  compute  the  wave  and  the  bore. 


dt  “  ”  4  dt 


+  -f-Co 


bCp 

6k 


by 


e  (c-Co) 


( t-cp)  Ic-CqI 
3  .  c  -t-  2Co  .2/3 

'  3  ' 


In  order  to  speed  up  computation,  we  now  make  two  assumptions. 
One  is  that  we  can  recognize  when  (  C  — C,  )  is  negative  or  positive 
so  that  we  can  simplify  the  last  term  in  the  equation;  the  other  and 
less  justified  assumption  is  that  (  C  —  Cp  )  has  a  weak  dependence  on  y 


so  that  the  second  term  can  be  written 

«  by  6  by 


We  now  set 


and  we  have 


Co(«-Co) 

6 


d' 

dt 

6 

6t 

+  (.-c)|j 

dc 

1 

dCo  3  d'co 

dt 

4 

dt  2  dt 

be  _  (<~Co)^r6g  2(3)^^\g~*^^  1 

by  ~S—  Iby  (e+Zc)2«  J 
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AIro  w«  recall: 


t»  _  2  .  .  2  ^®ob  *b 

'^b*  T‘b+T'«b 


Thus  we  can  see  that  the  total  computation  can  be  done  in  terms  of 

C|  ,  Co  and 

If- 

Thus  a  set  of  nomograms  is  indicated.  These  should  give: 

a. 

e-Co 

b. 

6 

c. 

(e-Co)* 

6 

d. 

i  .1 

2(3)’  kg’ 

- - J 

(e  +  2Co)* 

2  4Ca-e  2 

e. 

1  Co  ‘'V 

f. 

be 

57 

Some  of  these  quantities  should  be  expressed  as  numbers  and 
others  should  be  shown  as  the  distance  traveled  in  a  time  increment. 

The  nomogram  used  for  computation  is  included  in  Fig.  VII- 1. 
A  sample  approximate  computaliun  for  Chesapeake  Bay  is  included 
in  Fig.  Vri-2. 


3.  USE  OF  THE  NOMOGRAM 

This  nomogram  is  designed  to  compute  the  travel  of  a  wave  at 
speed  c  on  the  Coast  and  Geodetic  chart  No.  C  (c  GS  70-A  for 
Chesapeake  Bay  near  Norfolk.  The  travel  is  for  500  seconds.  The 


value  of  is  indicated  by  the  depdt  d  on  the  chart  and  it  is 
eesumed  that  values  of  c  are  written  on  the  chart.  The  wave 
front  traveling  at  speed  V  is  the  basis  of  Fig.  VII-2.  This  is 
very  much  like  a  wave  refraction  diagram. 

The  use  of  this  theory  and  nomogram  should  yield  interesting 
results  on  two  dimensional  wave  travel.  The  next  section  gives  a 
more  general  development. 


VllelO 


1 


1 


To  find  (-^<C)  enter  nomogram  with  (C  and  C^. 

1)  Follow  solid  line  to  intercept  on  axis  80/  dy  =  0. 

Z)  Follow  dashed  line  to  find  slope. 

3)  Draw  line  in  (-AC,  80/ By  )  plane. 

4)  Match  arcs  to  find  80/ 8y. 

5)  Read(-AC). 

FIGURE  VII-3 

NOMOGRAPH  FOR  TWO-DIMENSIONAL  NONLINEAR  LONG  WAVES 
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4. 


A  NUMERICAL  PROCEDURE  TOR  OBTAINING  AN  APPROXIMATE 
SOLUTION  FOR  THE  TWO-DIMENSIONAL  TIME  WATBR  WAVB 


PROBLEM 

The  following  ie  an  approximate  method  for  following  the  progreee 

of  a  solitary  wave  through  an  estuary  or  bay.  The  method  suggested 

is  a  slight  modification  of  "  The  Method  of  Wave  Derivatives,  "  proposed 

by  Freeman  and  Baer  (1957)  .  The  notation  used  here  ie  very  much 

the  same  as  that  used  in  the  above  reference  except  that  variable  depth 

and  friction  effects  are  taken  into  account  while  Coriolis  and  wind  effects 

are  ignored.  The  difference  equations  are  stated  in  terms  of  a  fixed 

space  grid,  but  only  those  portions  of  the  grid  in  which  the  effects 

of  a  solitary  wave  will  be  noticeable  are  used  in  the  computation. 

The  entire  procedure  is  motivated  by  physical  considerations  and 

die  need  to  obtain  reasonable  results  in  a  limited  amount  of  time. 

Consequently,  the  procedure  is  not  presented  with  an  appropriate 

analysis  of  stability.  We  apologize  in  advance  for  the  rather  glaring 

lack  of  mathematical  rigor  in  this  discussion. 

Let  D  be  a  bounded  domain  in  the  first  quadrant  of  the  x,  y 

plane.  Let  O*  be  a  domain  contained  in  D  such  that  the  distance 

e 

from  any  point  P  in  D  to  the  boundary  of  D  is  greater  than  some 
fixed  number  c  >  0  • 


% 


*  Freeman,  John  C. ,  Jr.  and  Ledolph  Baer.  "The  Method  of  Wave 
Derivatives."  Transactions,  American  Geophysical  Union, 

Vol.  38,  No.  4,  August  1957,  pp.  483-494. 
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FT.GURE  VU-5 

NOTATION  FOR  A  NUMERICAL  PROCEDURE 
FOR  NONLINEAR  LONG  WAVES 


Let  (x.,  y.)  be  a  square  grid  on  the  x,  y  plane  of  mesh  width  8  , 


We  are  concerned  with  the  system  of  partial  differential 


equations 


■(».>»  «0. '  >0.  (vn-4| 


vn-14 


4. 


where  u  and  v  are  the  x  and  y  components  of  the  fluid  velocity, 
d(x,  y)  is  the  depth  of  the  basin  represented  by  the  domain  O  from 
a  given  reference  level  (e.  g.  sea  level),  c  is  the  wave  velocity,  i.  e. 


c  * 


>/gh 


(VII.5) 


where 


h{  x,y,t)«  d(x,y)  +  »?(x,y,t) 


•tj  (x,  y,  t)  being  the  displacement  of  the  free  surface  relative  to 
the  reference  level.  The  quantities  and  are,  respectively, 

the  X  and  y  components  of  the  internal  friction  forces  per  unit 
mass.  These  are  assumed  to  be  of  the  forms 


Ex- 


S’ 


.2/3 


(VU.6) 


ri-l/3 

riTT 


vVu"  + V* 


The  values  u(x,y,  0),  v(x,  y,0),  d(x,  y)  and  (x,  y,  0)are 
assumed  to  be  known  in  D.  The  domain  D  or  equivalently,  the 
quantity  e  ,  is  to  be  determined  so  that  the  quantity  h  is  bounded 
away  from  zero.  For  example,  if  an  upper  bound 

M  »  8Up|  7;(x,y,t)| 

(x,y),D 
0<  t<  T 
e 

is  known,  then  the  domain  O  could  be  determined  by  the  requirement 
<■ 

that  if  (x,  y)  s  D  ,  then 


d(x,y)  —  M>  e>0, 
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where  •  ia  aome  pre-aaaigned  poaltive  number. 

By  the  methoda  of  1 1]  >  the  ayatem  VII- 1  may  be  replaced 

by  the  ayatem 


udiere 


"dT +  ■ 

« 

. .M  M 

at 

dt 

D4A4  . 
dt 

D,  A,  . 

O3  As 

dt  ' 

dt 

^2  4. 

dt  ^ 

D4  A4 

dt 

-0] 

b 

bd  5p 

/  bd  . 

'  bx  ^  by 

I  bd  bd 


dt  bf 


-4-(u+c)  +  (v+c) 


^.-|r+(u+c)|^+(v-c)|^ 


If-  -  +(u-c)  +(v-c) 


^  -  -|p  -H(u-c)  Ij-  +(v+c)^  . 


“Sy  ’ 

b 


(vn-7) 


(VU-8) 


A,  ■  u  +  V  +  2c  , 

A2  ■  u  -V  +  2c, 

►  (vi: 

As*  u+v-*2c, 

A4  •  u  -  v-2c, 

Di 

The  operatora  — gj —  are  cailled  "wave  derivatlvea"  and  the 
method  ia  a  generalization  of  the  familiar  method  of  characteriatica 


(VII-9) 
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A, 


( 


for  the  one  dimeneionel  time  dependent  problem. 

We  will  have  use  for  the  following  results  obtained  in  Freeman 
and  Baer 


^  (A|  -f  Ag)  -j-CAg-f  A4) 

U  «  ^ 

(A.  A^)-(A,+  A4) 

^  (A,+  A,)-(A3+  AJ 

2c  =  — ' - g— .---J - a-. 


(VII-IO) 


»  (A|  +  Aj)  ,  (A|  4-  Ag)  ”  (A3+  A4) 

A,*  '-V  - 4 - 

.  (A2+  A4)  .  (A,4-  A2)-(A3+  A4) 

Ag’  -  2  +  - 4 - ’ 

.  (A,  +  A,)  (A,+ A2)-(A3+ A4) 

A^»  - - - , 

(Ag+Aj  (A,+ A2)-(A3+ A4) 
^4*  2  ■  4 


We  will  replace  the  system  (VU-8)  by  a  system  of  difference 
equations  and  will  define  subregions  R(n  At),  n  =  0,  1, 

M  A  t  =  T,  of  the  domain  D*  in  which  our  computations  will  take 
place. 

The  region  R(0)  will  be  defined  by  our  initial  conditions  which 
will  describe  a  solitary  wave.  We  will  assume  that  for  points  (x,  y) 
outside  of  or  on  the  boundary  of  the  region  R(0), 

{u(x,y,0)  »  v(x,y,0)  »  0 

^U.yf0)«  0 
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FIGURE  VII- 6 

NOTATION  FOR  A  NUMERICAL  PROCEDURE 
FOR  NONLINEAR  LONG  WAVES 


We  define  R(0)  to  be  the  set  of  all  those  points  of  our  space  grid 
which  are  contained  in  R(0)  .  We  define  l[(  0)  to  be  all  those 
points  of  R(0)  which  are  at  a  vertex  of  a  square  contained  entirely 
in  R(0)  .  Analagous  definitions  will  be  assumed  for  other  regions 
in  the  x,  y  plane. 
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f* 


*• 


f 


Let  Py  a  (Xj,  y^)  and  let 

V,  (u,v,c)  «  (u  +  c,v  +c), 
Vg  (u,v,c)  «(u  +  c,v  -  c), 
Vj  (u,v,c)  »  (u-  C,  V  -  c), 
V4  (u,v,c)  ■  (u-  C,  V  +c), 


(VII -12) 


(VU-13) 

for  points  (Xj,  y^)  *]^'(0).  For  all  other  points  of  the  lattice  In  D*. 
define  all  of  these  quantities  to  be  equal  to  zero.  Extend  the  values 

to  all  points  of  the  piecewise  rectangular  region 
determined  by  1^(0)  by  linear  interpolation.  The  extended  functions 
u*(«,y),  v*(*,y),  will  then  be  continuous  on  the  piecewise 

rectangular  region  determined  by  iTlO). 

Define  R(At)  to  be  all  those  points  of  P  of  6*  which  satisfy 
the  requirements  that  the  points 


Also  let 

uj  •  U(x,  ,  yj.O), 
v;  »  v(x, ,  y,,0), 
“  ‘^<*1  .  yj  .0). 


pij  •  Pij . 
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are  in  the  piecewise  rectangular  region  determined  by  ^(0).  Define 
R  (At)  to  be  the  set  of  all  points  of  R  (At)  which  are  at  the 
vertex  of  a  square  contained  entirely  in  D 
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For  those  points  P,j«R*(Af)  we  can  solve  the  system  of 
difference  equations 

*2 ,1 + ^.1 = ^'  [  <  -,r  -  2  '  “j  ■'^;  •'^;  >  -  ^  s  ■' 

?;,)}]+r,(p,fi+'2;  <pr>’ 

bd(p„) 


^0 
V  * 

IJ 


wnere 


a  *  2g 
I) 


/9il  '  2« 


6k 

6d(P|j) 

br 


A„j  s  u,j -f-Vij  +  2c,j  t 


•y*  »N»*  »^4  1^9 

^21)*  +  » 


^3lj*  ^®l|  ’ 
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^9  ^9  ^4 

^41]  ■  “iJ  *■’'1)*”  ^-tj  ’ 

^9  0  ^9 

and  the  quantities  A/ {  P)  are  the  continuous  extensions  of  to 

H 

the  piecewise  rectangular  region  determined  by  R  (0)  .  The  quantities 

^I'l  '  "®i'|  '  ^n)'^'tl|**^'si|*^4IJ  respectively  determined  from 

Eqs.  (VII-10)  and(V’n-ll)  and  all  of  these  can  bo  extended  by  linear  inter¬ 
polation  to  functions  7»'(P),'v'(P)fc*(P),'Xj  (P),7(£(P),'Xj(P),7(J|(P) 


con 


tinuous  on  Ihe  piecewise  rectangular  region  determined  by  R  (At) 
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The  general  algorithm  for  continuing  our  computations  is  now 
obvious,  but  we  will  state  it  for  the  sake  of  completeness. 

Assume  that  the  quantities  'U||, . .  have  been 

H 

found  for  points  in  the  region  R  (kAt)  and  have  been  extended  by 
linear  interpolation  to  the  fimctions  . .  continuous 

n 

on  the  piecewise  rectangular  region  determined  by  R  (kAt)  ,  For 

points  P.j  of  D*  which  are  not  inR*(kAt)  ,  the  values  . •'^4ti 

are  to  be  obtained  by  setting  ‘ujj  *0  ►  Define  R((k+0At) 

to  be  all  those  points  P„  of  D*  which  satisfy  the  requirements  that 
the  points 


Vf  )  A  (vn-17) 

t 

are  in  the  piecewise  rectangular  region  determined  by  R(k  A  t  ) 

Define  r((K+I)  At)  to  be  the  set  of  all  those  points  of  R((k+I)  At) 
which  are  at  a  vertex  of  a  square  contained  entirely  in  D*. 

For  those  points  P^^  cR  (k  -t  l)At)  we  can  solve  the  system 
of 'difference  equations 


,wll+l  ,.^11+1 
Ai,j  +  Agjj 

At 

[«ll  -ZE. 

4wk  «N/k 

(  U,J  ,  V,j  ,  C|J( 

*^k 
+  Ag 

A  -1-  A  s 

3ir  ”4ij 

At 

[«„  -ZE. 

^k  ^k 
(u,j  ,  V,j  ,  C,J, 

^  k 
+ A4 

,  'J'S+I 

Aiij  +  A31J  » 

A  t 

k  ^k . 

V,j  ,  C,j) 

h 

II  • 

cr,)}]+K 

^k+l  ^N+l 

”2I|  <IJ 

At 

[(a,j-  ^,j) 

-z{e.  (<,. 

_  '>»k 

-E,  (u,j ,  v,j 

(pj)  +x;  (p*,V 
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The  quantities 


can  now  be  detennined 


from  Eq.  (Vn*10)and  Eq.  (VII- 11)  and  all  of  these  can  be  extended 
by  linear  interpolation  to  functions  u  (P)  , . .  ■ ,  (P)  continuous 

on  the  piecewise  rectangular  region  determined  by  R((k-f  I)  A  t) 

The  algorithm  is  now  completely  defined  and  the  computations  are 
iterated  M  times  until  the  final  time  M  At>  T  is  obtained. 


UST  OF  SYMBOLS 


X.  y 

u 

V 

t 

c 


g 

d 

V 

h 


n 


R 


horizontal  coordinates 

particle  velocity  along  the  x-aixis 

Particle  velocity  along  the  y-axis 

time  variable 

Vg  (d  +  i;  ) 

gravity  acceleration 

water  depth 

free  surface  elevation  around  the  still  water  level 
d+Tf 

friction  force  component  in  the  x  direction 

friction  force  component  in  the  y  direction 

angle  of  a  nearby  n  vector  with  the  u  vector 
with  which  we  are  working 

bore  velocity 

relative  to  the  bore  characteristics 
friction  coefficient 
3  c  -  2  Cq 

preassigned  positive  number 

wave  derivatives  operator 

u  +  v  +  2c 

integer 

domain 

a 

subregion  of  domain  D 
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